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DISCUSSION  AND  SUMMARY 


The  chairman  of  the  special  interest: -group on  control  theory  solicited 
from  those  attending  their  views  on  the  following  issues,  namely, 

a.  How  can  we  trigger  a  greater  DARCOM-wide  participation  to  this 
important  subject  to  the  Army. 

b.  Was  this  first  meeting  constructive  and  useful.  Was  there  a 
need  for  another  meeting.  If  so,  when,  where  and  who  should  attend. 

All  participants  considered  the  first  meeting  very  useful  addressing 
significant  problem  areas  in  fire  control  and  missile  guidance.  There 
was  unanimous  agreement  to  hold  a  second  meeting  in  early  FY81  (Nov/Dec  80) 
at  APG  and  general  agreement  to  invite  the  Navy  and  Air  Force  to  participate. 
Recommendation  was  made  and  favorably  indorsed  by  all  to  change  the  name  of 
the  group  to  "Coordinating:  Group  on  Modern  Control  Theory"  to  emphasize 
the  application  of  "modem  control  theory  and  techniques"  to  military 
problems  as  compared  to  classical  control  theory  techniques.  Suggestion  was 
made  that  the  chairman  explore  the  feasibility  of  using  the  JTCG/ME  by  which 
the  Navy  and  Air  Force  could  participate. 

It  was  strongly  felt  that  travel  fluids  severely  limited  number  of 
participants  and  that  DARCOM  HQ  support  should  be  solicited  to  encourage 
subordinate  commands  to  actively  participate  in  this  coordinating  group. 


RESEARCH  FOR  OS  ARMY  ADVANCED  G&C  SYSTEM 


Dr.  Harold  L.  Pastrick 
Guidance  &  Control  Directorate 
OS  Army  Missile  Laboratory 
Redstone  Arsenal,  Alabama  35809 


Abstract 

In  1973  the  US  Army  Missile  Commend  embarked 
upon  a  task  to  develop  an  advanced  guidance  and  con* 
trol  system  for  future  missiles.  It  was  Intended 
to  "leapfrog"  systems  currently  under  development  in 
order  to  meet  the  stringent  demands  and  constraints 
imposed  by  targets  with  predicted  characteristics  of 
the  1990s  and  beyond  and  by  the  predicated,  battle¬ 
field  environment  of  that  time. 

In  this  paper  the  problem  is  redefined  and  the 
latest  develocment  program  presented.  Results  that 
have  been  achieved  to  date  are  described,  particu¬ 
larly  in  the  areas  of  mathematical  models  of  the 
missiles  and  their  guidance  and  control  (G&C)  sys¬ 
tems  being  used  for  analyses  and  simulations;  aero¬ 
dynamics;  propulsion;  guidance  laws  being  developed 
ana  analyzed;  status  of  development  of  Disturbance 
Accommodating  Control;  signal  processing  to  locate 
and  track  the  target(s);  and  digital  design  tool 
development.  The  paper  is  concluded  with  a  section 
of  future  plans,  to  include  contractor  support. 

I.  Introduction 

The  US  Army  Missile  Command  (MICOM)  recently 
began  a  task  to  develop  an  advanced  GAC  system  for 
future  Army  Modular  Missiles.  The  intent  is  to 
"leapfrog"  systems  currently  under  development. 

The  purpose  of  this  paper  Is  to  describe  the  work 
that  has  been  completed  within  this  new  task  and  to 
provide  an  indication  of  future  efforts  that  are 
now  planned.- 

The  first  step  in  implementing  this  tesk  was 
to  conduct  a  literature  survey  to  estaollsh  a  tech¬ 
nology  base  starting  point.  Following  this  survey, 
guiaance  laws  ware-placed  in  five  categories  and 
defined  mathematically.  The  Impl wientation  and 
predicted  performance  of  each  category  was  then  in¬ 
vestigated  and  compared  in  light  of  current  and 
predicted  hardware  and  software  capabilities.  ' 

This  work  was  subsequently  updated  In  1979.  2 

The  program  objectives  are  three-fold: 

•  To  develop  and  prove  a  SAC  system  that  is 
capable  of  guiding  and  controlling  future  US  Army 
missiles  (generelly  defined  as  air  defense  and  sur- 
*ace-ta-surfaee  general  suoport)  to  destroy  pre¬ 
scribed  lines  of  future  targets.  This  must  be 
scconpl Ished  under  the  predicted  severe  battlefield 
jnvlrormenc  of  the  future. 

•  To  "leapfrog"  systems  under  current  develop¬ 
ment. 


•~i's  oaoer  Is  declared  a  work  of  the  US  Govern¬ 
ment  and  therefore  is  in  the  public  domain. 


’"Consulting  Engineer,  Associate  Fellow,  AIAA. 


e  To  broaden  and  deepen  the  existing  GSC  sys¬ 
tem  technology  and  design  base  within  the  GSC 
Directorate  of  MICCM. 

In  the  sequel  the  expected  threat  and  its 
characteristics  will  be  surma ri zed .  This  is  follow¬ 
ed  by  a  description  of  the  development  plan  for  the 
advanced  GSC  system.  The  progress  to  date  is  pre¬ 
sented,  and  the  paper  is  concludea  with  a  section 
on  future  plans. 

II.  Expected  Threat 

Theater  defense  typically  is  provided  by  a 
mixture  of  ground-based  and  airborne  defense  systems 
supported  by  radars,  command  and  control  systems, 
electronic  warfire  equipment,  and  passive  measures 
such  as  camouflage,  decoys,  and  equipment  disper¬ 
sion.  The  air  defense  objective  of  ground  based 
systems  is  to  limit  the  opponent's  effectiveness  by 
attacking  Ms  critical  assets  so  that  land  ‘orces 
may  maneuver  with  a  minima  of  Interference  from 
the  enemy  air  weaponry. 

For  many  years  now,  the  enemy  doctrine  has 
emphasized  large  mass  and  brute  force,  and  his  air 
attacks  will  provide  no  exception.  It  is  entirely 
feasible  to  assume  that  an  attack  in  the  Central 
Europe  area  will  ce  accompanied  by  several  thousand 
combat  aircraft,  a  rft  adoitlon,  his  doctrine  calls 
for  the  massing  of  large  quantities  of  artillery 
fire  on  a  section  selected  'or  a  tank-led  break¬ 
through.  It  is  unlikely  that  'IAT0  forces  either 
now  or  in  the  near  future  will  match  the  Warsaw 
Pact  forces  in. terms  of  numbers  of  weapons,  nor  is 
it  the  intent  to  aim  toward  that  end.  Rather,  it  is 
important  to  optimize  the  effectiveness  of  our3 
smaller  force  to  meet  the  anticipated  threat. 

The  Army  is  attempting  to  maximize  the  effec¬ 
tiveness  of  its  current  family  of  air  defense 
weapons  wnile  concurrently  developing  a  new  family 
to  meet  the  threat  of  the  1390's.  In  the  near 
term,  there  will  be  continued  modification  of  cur¬ 
rent  systems  as  necessary,  and  while  still  feasible, 
to  overcome  qualitative  and  quantitative  deficien¬ 
cies.  Longer-term  replacements  continue  in  devel¬ 
opment  or  procurement  for  all  the  major  field  army 
air  defense  systems.  Examples  of  this  strategy 
include  the  following:  high  to  medium  altitude 
missile  systems  -  PATRIOT  ‘or  MIKE  HERCULES  ana 
HAWK;  short  range  missile  systns  -  U.  S.  ROLAND  ‘or 
CHAPARRAL;  manoortable  mi ss i 1 es  STINGER  for 
REDEYE;  mobile  gun  systems  -  DIVAD  given  to  VULCAN. 
The  systems  will  provide  the  effective  aerial 
umbrella  needed  by  our  forces  to  not  only  survive, 
but  to  fight  effectively. 

For  security  reasons,  it  is  impossible  in  this 
forum  to  describe  specifically  the  air  threat  that 
will  be  encountered  in  the  scenario  described  above. 
However,  in  order  to  auantify  the  problem  somewnat, 
we  shall  attempt  to  attribute  vehicle  characteris¬ 
tics  to  the  enemy  based  on  our  current  technology 
in  the  field  of  air  defense  targets .  The  Targets 
Management  Office,  'JS  Army  Missile  Commend  since 
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1964  has  published  an  extensive  library  of  target 
program  reports.  In  particular,  they  have  classi¬ 
fied  a  variety  of  aerial  targets  as  test  and  evalua¬ 
tion  (TS£)  targets  used  for  air  defense  weapon 
systems.  Needless  to  say,  targets  used  for  this 
purpose  must  exercise  an  air  defense  weapon  system 
to  the  limits  of  its  capabilities.  5 

A  particularly  interesting  target  is  known  as 
HAHST,  an  acronym  for  High  Altitude  High  Speed 
Target.  It  is  designed  to  achieve  speeds  up  to 
Mach  4  at  altitudes  up  to  100,000  feet.  Additional 
performance  characteristics  for  HAHST  and  other 
existing  targets  are  given  in  References  2  and  S. 

The  effectiveness  of  any  missile  systm  is  con¬ 
ditioned  on  its  ability  to  function  in  an  Electronic 
Countermeasures  environment.  Stand  off  jammers, 
barrage  jammers  and  even  dispensed  chaff  will  be 
used  to  deny  the  air  defense  tracking  radars  the 
capabilities  needed  to  be  effective  in  their  sector. 
The  jammers  are  Intended  to  reduce  the  acquisition 
range  of  the  radars  and,  if  perfected,  will  elimi¬ 
nate  accurate  tracking  entirely.  Additionally,  the 
threat,  aircraft  and  missiles  will  be  fabricated  to 
present  the  smallest  possible  radar  cross  section. 
The  state-of-the-art  in  this  field  is  beyond  the 
scope  of  presentation  in  this  paper. 


utilize  different  modules  such  as  propulsion,  guid¬ 
ance,  warhead,  etc. 

A  program  plan  initiated  by  the  GSC  Director¬ 
ate,  MtCCM,  was  undertaken  approximately  two  years 
ago  with  the  program  objectives  enumerated  above. 
The  program  plan  contains  six  intermediate  objec¬ 
tives,  the  end  of  which  each  constitutes  a  program 
milestone.  They  are: 

e  Define  overall  program; 

e  Collect  elements  that  will  form  the  candi¬ 
date  GSC  systems; 

e  Define  candidate  GSC  systems; 

e  Evaluate  the  candidate  GSC  systems: 
select  best  system; 

e  Design  and  fabricate  the  selected  system; 

e  Demonstrate  "Prcof-of-Concept"  of  the 
selected  system. 

These  separate  program  elements  (or  intermediate) 
objectives  are  described  pelow  in  more  detail.  An 
accompanying  milestone  chart  is  provided  as  Fig.  1. 


From  a  defensive  viewpoint,  the  effect  of 
enemy  j  aiming  of  the  air  defense  radars  and  the 
minimization  of  enemy  attack  aircraft  radar  cross 
sections  have  a  profound  Impact  on  the  air  defense 
missile  system.  The  rationale  is  reasonable  and 
straightforward.  If  the  enemy  does  indeed  have 
aircraft  and  attack  missiles  and  remote  pilotless 
vehicles  (RPVs)  either  with  or  better  then  the 
characteristics  attributed  to  them  via  the  method 
atove,  and  if  the  enmity  minimizes  his  radar  cross 
section  to  the  current  state-of-the-art,  his  attack 
vehicles  will  be  extremely  difficult  to  acquire  at 
long  range.  The  affect  of  the  combined  high  speed 
and  high  agility  (l.e.,  high  g -maneuver  capability) 
with  low  radar  cross  section  yields  precious  little 
reaction  time  to  the  air  defense  system.  The  close- 
in  acquisition  will  seriously  degrade  the  existing 
air  defense  missile's  GSC  system’s  performance, 
since  most  are  based  an  a  proportional  navigation 
and  guidance  (PNG)  law.  An  environment  such  as 
presented  above,  however,  can  be  better  addressed 
in  terms  of  guidance  laws  explicitly  tailored  to 
this  type  of  threat.  Thus,  the  optimally  guidad 
and  controlled,  highly  maneuvering,  defensive  mis¬ 
sile  using  terminal  guidance  sensors  chosen  from 
across  a  wide  range  of  the  frequency  spectrum  must 
be  initiated  into  the  development  cycle.  As  a 
necessary  first  step,  the  research  described  in  the 
sequel  addresses  that  problem. 


Development  Program  Plan 


From  an  overall  systems  viewpoint,  this  pro- 
gran  shall  address  the  issue  of  creating  new  theory 
in  :he  GSC  area  to  meet  the  high  performance  threat 
of  the  future  as  one  of  the  leading  technology 
items.  Closely  associated  with  it  and  in  parallel 
with  the  GSC  effort,  weapon  system  work  shall  be 
undertaken  to  modify  airframe  and  propulsion  to  be 
capable  of  engaging  the  threat  of  the  1990' s.  Gen¬ 
eral  support  weapons  shall  be  viewed  initially  as  a 
subset  of  the  air  defense  systtm(s),  whereas  pre¬ 
viously,  these  two  classes  of  weapons  were  develop¬ 
ed  independently.  This  research  shall  attempt  to 
view  them  as  potentially  similar  systams  that 


NO.  TASK 

1  DEFINE  PROGRAM 

2  COLLECT  GSC  ELEMENTS 


3  DEFINE  CANDIDATE  SrSTEMStZTJ 

4  EVALUATE  CANDIDATES  d 

5  DESIGN  S  FA3RICA7E 


6  PROOF  OF  CONCEPT 


FY/9  FY80  FY81  FY82  FY83  FY84 


Figure  1.  Program  Milestone  Chart. 


Task  1.  Define  Overall  Program 


In  this  first  task  it  is  necessary  to  define 
the  program  objectives,  goals,  and  constraints. 

This  is  followed  by  a  survey  of  representative  mis¬ 
sile  plants  and  the  characteristics  of  sensors  and 
effectors  that  sight  be  used  (either  available  or 
under  development).  Terget  dynamics  and  initial 
conditions  must  be  defined,  as  must  missions  for 
the  missile  system(s)  using  the  advanced  GSC  sys¬ 
tems  that  emanates  from  this  program.  Also  includ¬ 
ed  tn  the  first  program  element  is  the  beginning  of 
coordination  with  other  missile-developing  agencies 
and  services.  Finally,  the  time  frame  predicted 
usage  of  the  GSC  systam  must  be  defined.  As  des¬ 
cribed  in  this  paper,  most  of  Task  1  has  been  com¬ 
pleted. 


Task  2.  Collect  Elements  for  Candidate  GSC  Systems 


In  this  task  detailed  definitions  of  the 
characteristics  of  candidate  guidance  laws,  candi¬ 
date  autopilots  (tontrol  laws),  state  estimation 
techniques,  state  truncation  techniques,  and  Dis¬ 
turbance  Aeecncodation  Control  (DAC)  shall  be 
accomplished.  s  Further,  characteristics  of  all 


expected  disturbances ,  both  externel  (to  tha  mis¬ 
sile)  and  Internal,  oust  be  collected  and  evaluated. 
The  need  for  OAC  and  for  state  estimation  will  be 
evaluated.  Characteristics  of  all  expected  signify 
leant  system  nonlinearities  will  be  collected  and 
their  dynamic  Importance  evaluated. 

Task  3.  Define  Candidate  GBC  Systems 

The  components  and  subsystems  that  have  been 
identified  In  the  earlier  tasks  will  now  be  combin¬ 
ed  into  candidate  QIC  systems.  Missile  structural 
dynamics  wilt  be  determined  (if  not  already  com¬ 
pleted  within  tha  plant  definition  of  the  earlier 
tasks)  or  refined,  as  necessary.  Tha  simulation 
program  objectives  will  be  defined;  the  simulation 
program  will  then  be  defined;  and  finally,  develop¬ 
ment  of  the  final  system  simulation  will  begin. 
Investigations  of  the  dynamics  of  Che  candidate  SAC 
systems  will  begin  using  both  mathematical  analysis 
and  computer  simulation.  A  figure  of  merit  (cost 
functional)  will  be  developed  during  this  task,  as 
will  model  error  effects  and  criteria. 

Task  <■  Evaluate  Candidate  SAC  Systems 

The  development  of  the  system  simulation  and 
the  investigation  of  the  candidate  GAC  systems, 
begun  in  Task  3,  will  be  completed.  The  candidate 
GAC  systems  then  will  be  evaluated  with  respect  to 
the  Figure  of  Merit  developed  in  Task  3,  using 
mathematical  analysis  and  computer  simulation.  The 
strengths  and  weaknesses  of  distributed  versus  cen- 
tralized  controllers  will  be  assessed,  probably  as 
enhanced  by  the  use  of  microcomputers .  Finally, 
the  best  GAC  system  will  be  selected  and  Integrated 
into  a  missile  airframe. 

Task  5,  Design  and  Fabrication 

A  detailed  design  of  the  selected  SAC  system 
will  be  performed.  Tha  testbed(s)  selected  for  use 
in  the  ?roof-o f-Concapt  phase  will  be  fabricated 
and/or  assembled. 

Task  5,  Proof-of-Concept 

The  Proof-of-Concept  will  be  demonstrated  with 
hardware  firings  during  this  task.  These  firings 
will  be  augnented  by  simulations  (computer  and 
hardware-in-the-loop  as  deemed  necessary)  and  anal¬ 
yses  as  necessary. 

IV.  Progress  to  Date 

The  accompl ishnents  to  date  have  been  achieved 
by  elements  of  MtCOM  and  through  the  use  of  re¬ 
search  contracts.  It  has  been  supplemented  through 
coordination  with  the  US  Air  Force  Armaments  Labo¬ 
ratory  at  Eg! In  Air  Force  Base,  the  US  Army  Bal¬ 
listic  Missile  Defense  Systems  Command  at  Hunts¬ 
ville,  Alabama,  and  the  Office  of  the  Under  Secre¬ 
tary  of  Defense.  This  coordination  has  been  and 
will  continue  to  be  carried  out  to  eliminate  dupli¬ 
cated  development  effort  on  similar  projects  within 
the  Department  of  Defense.  To  augoent  the  in-house 
research  and  engineering  capability,  several  re¬ 
search  contracts  have  been  initiated  in  specialized 
areas  currently  including:  the  Computer  Sciences 
Corporation,  Huntsville,  Alabama;  the  Dynamic 
Systems  Aesearch  and  Training  Corporation,  Hunts¬ 
ville,  Alabama;  the  University  of  Florida;  Western 
Kentucky  University;  and  the  Control  Dynamics  Com¬ 
pany,  Huntsville,  Alabama.  Additional  contracts 


with  other  organizations  are  anticipated  as  the 
scope  of  the  program  grows . 

The  technical  work  comprising  Tasks  1  and  2 
has  been  apportioned  to  various  members  of  the 
MICOM  and  contractor  team.  As  indicated,  most  of 
Task  1  has  been  completed,  and  effort  on  Task  2 
work  is  underway. 

A.  Target  Definition 

A  comprehensive  investigation  of  predicted 
future  targets  and  their  dynamics  has  been  complet¬ 
ed.  This  investigation  included  reviewing  and  dis¬ 
cussing  matarial  available  within  the  sources  of 
the  US  Army  and  US  Air  Force  and  Included  Inputs 
from  several  industrial  organizations.  Finally, 
the  collected  information  was  discussed  with  ele¬ 
ments  of  the  Office  of  the  Under  Secretary  of 
Dafensa,  Research  and  Engineering.  While  most  of 
the  collected  information  is  classified  for  secu¬ 
rity  reasons,  tha  results  of  this  most  important 
phase  of  the  study  leed  to  emphasizing  targets  in 
three  categories:  highly  maneuverable  aircraft, 
cruise  missiles,  and  tactical  ballistic  missiles. 
Not  addressed  in  this  study  are  RPV's  and  surface 
targets.  It  is  felt  that  use  of  a  sophisticated 
missile  system  to  engage  numerous  RPV's  would  not 
be  cost  effective  and  that  missiles  currently 
under  development  will  be  able  to  combat  surface 
targets.  Except  for  periodic  updating,  this  sub¬ 
task  is  now  complete. 

B.  Missile  Plant  Characteristics 

In  order  to  get  the  other  subtasks  underway. 
Standard  equations  of  motion  for  missile  bodies 
with  some  Structural  flexibility  have  been  used. 

The  missiles  are  assumed  to  be  acted  upon  by  the 
usual  aerodynamic  forces  and  torques.  Currently, 
variations  of  the  SPRINT  missile  family  are  among 
the  principal  contenders  for  the  airframe  and  pro¬ 
pulsion  system. 

C.  Aerodynamics 

A  detailed  plan  is  being  generated  for  future 
actions  in  gathering  and  generating  aerodynamic 
data  to  be  used  in  this  program.  Novel  aerodynamic 
shapes  arc  under  consideration  and  evaluation  for 
use  in  developing  control  authority  for  missiles 
using  advanced  5SC  systems.  These  include  missiles 
using  configurations  shown  in  Table  1  (along  with 
thoir  cheracterlstics) . 

D.  Error  Sources 

Expected  error  sources  have  been  categorized 
into  five  detailed  groups,  to  Include  predicted 
mean  and  standard  deviation  values.  Several  are 
presented  as  representative  of  the  groups.  These 
will,  of  course,  change  as  the  choices  are  narrowed 
and  the  program  develops. 

1.  Missile.  Included  in  this  group  are  para¬ 
meter  uncertainties  such  as  thrust  magnitude  and 
misalignment,  C.S.  location  and  offset,  mass,  and 
transverse  and  axial  moments  of  inertia. 

2,  Aerodynamic.  Normal  and  axial  forces; 
pitch,  yew,  and  ro) 1  moments;  pitch  end  yaw  damping 
derivatives;  roll  damoing  coefficient  due  to  fin 
deflection;  and  axial  drag  are  factors. 
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Tibi*  1.  Aerodynamic  Configuration*  and  Characteristics. 


Configuration 

Advantages 

Disadvantages 

Mach 

NO. 

Altitude 

(ft.) 

Maneuverability 

(g'*) 

1 .  Deployable  wings 
with  all-moveable 
tails 

a.  Simple  con¬ 
trols 

b.  Good  roll 
control 

a.  Control  force 
in  opposite  direc¬ 
tion  from  maneuver 

b.  High  angle  of 
attack 

2-4 

0-100  K 

10-15 

2.  Low  aspect 
ratio,  long-chord 
delta  wing  with 
tab  controls 

a.  Simple  con¬ 
trols 

b.  Good  roll 
control 

a.  Control  force 

In  opposite  direc¬ 
tion  from  maneuver 

b.  High  angle  of 
attack 

4-7 

0-140  K 

20-30 

3.  Flared  skirt- 
stabilized  missile 
with  all  moveable 
wings 

a.  Control  force 
In  direction  of 
maneuver 

b.  Lower  angle 
of  Attack 

Possible  higher 
drag  and  hinge 
moment 

4-7 

0-140  1C 

20-30 

4.  Reentry  body 
shape  with  moveable 
wedge  controls  * 

a.  Simple  shape 

b.  Good  roll 
control 

Roll  control  sur¬ 
faces  separate 
from  maneuvering 
controls 

3-10 

7 

S.  Lifting  body  * 

a.  Lower  angle 
of  attack 

b.  Control  fore* 
In  direction  of 

maneuver 

Complicated  aero¬ 
dynamic  shape 

3-6 

0-90  K 

10 

*  Comment.  Must  uso  bank-to-tum  guldanct. 


3.  Instruments 

*  Accel  erometers  -  *ca1«  factor  stabll  1t>>. 
bias  stability,  non-orthogonality,  *g«"  -  seal* 
factor,  third-order  seal*  factor,  cross-axis  sensi¬ 
tivity.  cross-coupling,  seal*  factor  asymmetry,  and 
rectification  error. 

*  Gyro  scones  -  seal*  factor  stability; 
bias  stability;  non-orthogonality;  anlsoelastlc 
drift;  drift  rates  In  pitch,  yaw,  and  roll  due  to 
input  axis  and  spin  axis  accelerations  as  well  as 
those  (rates)  independent  of  acceleration,  due  to 
torquer  nonl Inearltles,  and  due  to  electronic 
noise;  and  mass  unbalance. 

*  Porro-orlsa  azimuth  alignment. 

*  Laser  Inertial  measurement  unit  mis¬ 
alignment  with  respect  to  missile  body  axes. 

*  Accelerometer  triad  origin  displacement. 

e  Uplink/downlink  bias  nnd  calibration 

errors . 

e  Optical  correlator  errors. 

t  Radar  -  rang*  track  and  angle  track 
noise  and  accuracy,  ground  clutter  noise,  and  tar¬ 
get  glint  noise. 

4.  External.  Wind  magnitude  and  direction 
(Initial  azimuth,  elevation,  and  roll  alignment  of 


the  missile;  target  velocity  and  Illumination 
Jitter;  samlactlve  laser  pointing  accuracy  and  beam 
divergence;  gravity  bias;  atmospheric  effects 
(e.g.,  upon  radio  range). 

5.  Subsystems 

*  Common  nonl inearltles  -  saturation, 
coulomb  (and  other)  friction,  backlash,  and  bang- 
bang  with  dead  zone. 

a  Computer  -  quantization,  truncation,  and 
fixed  word  length. 

e  Seeker  -  boreslght  error  (In  pitch  and 
yaw)  due  to  servo  noise;  channel  crosscoupling; 
coupling  between  the  seeker  head  and  the  airframe; 
effects  of  the  radome  and  Irdome  on  angle  linear¬ 
ity;  angle  bias  (the  electrical  equivalent  of  mech¬ 
anical  BSE);  gain  stability;  angle  noise;  and  bore- 
sight  error  In  pitch  and  yaw  due  to  clutter,  re¬ 
ceiver,  and  jamming  noise. 

e  Autopilot  -  bias  errors,  time  delays, 
and  gain  stability. 

*  Guidance  -  errors  In  Initial  position, 
velocity,  and  acceleration. 

As  the  program  progresses,  these  error  sources 
will  be  analyzed  further  to  determine  which  might 
be  amenable  to  cancellation  by  appropriate  disturb¬ 
ance  accofimodatlon  design  theory  and  techniques. 
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Z.  Guidance  Laws 

This  is  a  major  thrust  araa  within  the  program. 
As  indicated  above,  an  extensive  literature  search 
has  been  completed  and  documented.  This  was  follow¬ 
ed  by  placing  guidance  laws  in  five  categories  and 
describing  each  mathematically.  The  implementation 
and  predicted  performance  of  each  category  has  been 
initially  reported  in  Reference  1  and  subsequently 
updated  in  Reference  2.  A  summary  of  the  latter 
is  shown  in  Table  2. 

Investigations  continue  into  areas  of  optimal 
guidance,  In  particular,  emphasizing  digital  as¬ 
pects  in  anticipation  of  the  expected  use  of  on¬ 
board  digital  controllers.  Constant  effort  is  made 
to  reduce  implementation  complexity,  where  complex¬ 
ity  is  defined  as  requirements  for  hardware  and 
software.  Other  innovative  techniques  are  under 
Investigation.  They  Include  assessing  the  potential 
application  of  Singular  Perturbation  Theory,  Distur¬ 
bance  Accommodating  Theory,  and  means  of  determining 
or  predicting  the  very  important  (for  optimal 


applications)  quantity,  time-to-go,  i.e.,  remaining 
time  of  flight  at  any  Instant  (see  Table  3). 

1  .  Optimal  Control .  A  conventional  implement*  - 
tion  of  Linear-Quaoratic  Optimal  Control  is  por¬ 
trayed  graphically  in  Figure  2.  The  optimal  con¬ 
trol  authority,  uj*q,  is  selected  to  be 

u»Q  •  -R-’s^i  .  (1) 

The  matric  Riceati  equation  is  solved  to  obtain  the 
control  gain,  *x(t),  while  minimizing  a  quadratic 

performance  index. 

2.  Disturbance  Accorrodatinq  Control.  The 
Disturoance  Accommodating  Control  (DAC)  theory  is 
described  for  the  continuous-time  domain  in  Refer¬ 
ence  S.  The  general  nature  of  the  DAC  controller 
is  to  generate  a  real-time  on-line  estimate  of  the 
actual  (instantaneous)  disturbance  waveform  and 
Create  a  special  control  action  that  exactly 


Table  2.  Convention*!  Suidance 


Approach 

Advantages 

Limitations 

1.  Attitude  Pursuit 

a.  Simplest  Implementation 

b.  Fixed  targets 

Sensitive  to  target  velocity, 
disturbances 

2.  Velocity  Pursuit 

a.  Simple  implementation 

b.  Hon-maneuverl ng  targets 

Sensitive  to  target  accelera¬ 
tion,  disturbances 

1 

;  3.  Proportional  .'laviga- 
j  tlon  Guidance  (PNG) 

a.  Simole  implementation 

b.  Maneuvering  targets 

Sensitive  to  high  end  game 
maneuvers 

Table  3.  New  Guidance  Approaches 


Approach 

Advantages 

Limitations 

1.  Linear  Quadratic  (LQ) 
Regulator 

Better  than  PNG  against 
maneuvering  target 

a.  Tgg  estimate  required 

b.  Must  comoute  time-varying 
control  gain  (Kx) 

c.  Disturbances  Ignored 

2.  Linear  Quadratic 

Gaussian  (LQG) 

Regulator 

Better  than  LQ  against 
noise-type  disturbances 

a.  Tgg  estimate  required 

b.  Must  compute  time-varying 
control  gain  (1^) 

c.  Must  eomoute  Kalman  gain 
for  estimator 

3.  Disturbance-Utiliz¬ 
ing  Control  (OUC) 

Better  than  LQ  or  LQG 
against  waveform- type 
disturbances 

a.  Tgg  estimate  required 

b.  Must  conoute  2  time-vary¬ 
ing  control  gains  {<x>  KxI) 

4.  Singular  Perturbations 


Computational  efficiency 


An  approximation  to  optimal 
control 


cancel  s-out  th*  disturbance  effect  on  the  missile. 
The  DA C  theory  will  bo  extended  Into  •  discrete- 
time  domain.  To  date,  tho  dost  of  syst«s  and 
disturbances  amenable  to  a  discrete-time  version  of 
OAC  have  boon  defined.  Work  Is  underway  on  a  des¬ 
cription  of  the  state-reconstructor  that  will  be 
associated  with  a  digital  OAC  controller.  Also 
under  consideration  is  the  possibility  of  applying 
OAC  to  various  missile  subsystem  or  coeiponent  out¬ 
puts.  such  as  sensors,  whose  normal  outputs  have 
been  modified  by  the  Influence  of  the  disturbances 
upon  the  sensors. 

An  innovative  modification  of  OAC  is  01s- 
turbance-UtlllzIng  Optimal  Control  (OUC)7.  In  this 
case  waveform-type  disturbances  are  exploited  opti¬ 
mally.  Examples  of  such  disturbances  are  drag, 
target  maneuvers,  wind  gusts,  any  effects  of  the 
gravitational  field.  A  graphical  portrayal  of  OtIC 
implementation  is  shown  In  Figure  3.  In  this  ease, 
the  optimal  control  authority  Is  specified  as 

u0UC  *  *R  9  (Kxx  *  Kxz*>  •  ^ 

where  Kx(t)  Is  found  by  solving  a  metric  Riceatl 
equation,  and  1C  (t)  Is  found  by  solving  a  linear 
differential  equation.  To  date,  computer  simula¬ 
tions  have  shown  OUC  to  be  quite  effective  when 
compared  to  the  performance  achieved  by  using  con¬ 
ventional  LQ  controllers. 


Disturbance 


3.  Singular  Perturbations.  Appl  teat  ion  of 
Singular  Perturbation  Tneory  to  missile  control  may 
be  attractive  if  It  is  damned  necessary  for  the  3 
control  law  to  account  for  high  order  model  terms. 
The  standard  approach  is  to  approximate  the  model 
with  relatively  lower  order  equations.  However, 
the  neglected  higher  order  terms  may  be  dynamically 
significant.  Because  their  inclusion  might  create 
computational  problems,  a  possible  alternative 
might  be  the  application  of  Singular  Perturbation 
Theory.  The  application  is  particularly  amenable 
to  controller  design  where  the  open-loop  plant  has 
a  wide  eigenvalue  dispersion,  slow  and  fast  modes, 
or  parasitic  parameters. 

F.  Sensor  Characteristics 

Various  existing  and  predicted  sensors  have 
been  characterized.  Trade  studies  to  aid  In  their 
selection  have  been  identified.  It  now  appeers 
that  a  sensor  using  some  form  of  pattern  recogni¬ 
tion  may  be  required.  The  problem  is  to  find  or 
begin  developing  a  sensor  that  can  provide  a  guid¬ 
ance  signal  with  a  superior  signal -to-nolse  ratio 
from  data  that  has  been  del iberately  modified,  for 
example,  by  high  powered  jamming  equipment. 

6.  Computer-Aided  Desien  Tools 

A  number  of  computer  programs  are  available 
(such  as  root  locus  and  other-  frequency  domain 


Figure  2.  Conventional  Linear  Quadratic  Optimal  Control. 


Disturbance 

w(t) 


Figure  3.  Olsturbence-Utlllzlng  Optimal  Control. 
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techniques)  to  «1d  in  design  of  control  systems. 

Most  of  them  ere  In  the  continuous-time  domain, 
although  there  does  exist  a  z-transform  manipula¬ 
tion  program.  An  executive  routine  currently  is 
being  coded  to  manipulate  these  programs  effi¬ 
ciently. 

H.  Microcomputer  State-of-the-Art 

A  continuing  assessment  of  microcomputer  state- 
of-the-art  has  been  Instigated  within  this  program. 
Microcomputers  will  became  an  Integral  part  of 
advanced  SAC  systems,  with  their  small  size,  light 
weight,  and  relatively  low  cost.  Presently,  micro¬ 
computers  software  capabilities  are  being  investi¬ 
gated  with  a  Tektronix  3002  Microprocessor  Labo¬ 
ratory.  This  investigation  will  be  followed  by  an 
analysis, of  the  computational  capability  to  meet 
the  requirements  Imposed  by  an  advanced  SAC  system. 

I.  Oiqital  Design  Tool  Development 

The  probability  that  an  advanced  SAC  system 
will  be  implemented  digitally  seems  to  be  nearly 
assured,  tt  is  for  this  reason  that  the  foregoing 
work  has  been  oriented  strongly  in  the  direction 
of  digital  implementation  and  its  consequences.  A 
review  of  analysis  and  design  tools  available  to 
the  SAC  system  designer  is  being  conducted.  Where 
the  need  for  a  new  tool  appears  to  be  warranted,  it 
will  be  developed  within  this  program  when  possible. 
Emphasis  is  placed  on  both  simplicity  of  applica¬ 
tion  and  being  able  to  draw  on  practising  engineers 
existing  engineering  training  and  experience. 

One  example  of  such  a  tool  is  “SAM"  (acronym 
for  Systematic  Analytical  Method).  SAM  provides  an 
alternative  to  the  use  of  signal  flow  graphs  and 
the  application  of  Mason's  Sain  Rule  to  determine 
selected  states  of  the  missile.9  The  technique  is 
particularly  useful  for  analysis  of  complicated 
sampled-data  control  systems.  An  advantage  of 
using  SAM  is  that  the  cumbersome  application  of 
Mason's  Gain  Formula  can  be  avoided.  Further,  the 
entire  method  of  constructing  signal  flow  graphs 
may  be  circumvented.  Since  only  the  equations  des¬ 
cribing  the  system  are  needed  for  SAM,  even  the 
customary  block  diagram  is  not  needed.  The  tech¬ 
nique  Is  analytical  in  nature  and  makes  use  of  a 
systematic  manipulation  of  the  system  algebraic 
equations.  These  manipulations  follow  prescribed 
rules  set  forth  in  the  technique. 

A  second  simplified  technique  is  the  determina¬ 
tion  of  digital  control  system  response  by  cross- 
multiplication.10  This  technique  permits  the 
analyst  to  obtain  the  response  (at  the  samp! in g  in¬ 
stants)  of  any  system  state  from  its  closed-loop 
transfer  function  expressed  in  the  complex  z-domain. 
If  It  is  desired  to  know  the  response  between 
sampling  Instants,  either  the  submul ti pi e  method  or 
the  modified  z-transform  method  may  be  adapted  to 
the  cross-multiplication  technique. 

A  third,  more  sophisticated  technique  1$  the 
Parameter  Space  Method.  It  is  being  developed  for 
determining  the  stability  and  dynamic  characteris¬ 
tics  of  a  digital  control  system  in  terms  of  sever¬ 
al  selected  system  parameters."  The  method  re¬ 
quires  that  the  system  characteristic  equation  be 
available  in  the  complex  z-domain.  Although  not 
necessary,  its  application  is  facilitated  by  aug¬ 
menting  the  analytical  results  with  graphical  por¬ 
trayals  in  a  selected  multi  parameter  space.  The 


method  is  based  on  analysis  and  synthesis  methods 
for  linear  and  nonlinear  control  system  design 
which  are  amply  described  in  Siljak's  excellent 
monograph  on  the  subject. in  essence,  the  para- 
meter  space  method  permits  the  designer  to  evaluate 
graphically  the  effects  of  the  locations  of  the 
roots  of  the  characteristic  equation.  Hence,  he 
may  design  the  control  system  in  terms  of  his  se¬ 
lected  performance  criteria;  e.g.,  absolute  stabil¬ 
ity,  damping  ratio,  settling  time.  He  is  able  to 
set  the  effects  on  the  characteristic  equation 
roots  (and  hence  on  system  dynamics)  of  changing 
several  adjustable  parameters.  The  method  has  been 
extended  to  portray  the  effect  of  varying  the 
sampling  period,  thereby  permitting  one  to  observe 
the  effect  of  the  choice  of  various  values  assigned 
to  the  sampling  period  on  absolute  and  relative 
stability.  Also,  simple  recursive  formulas  have 
been  derived  so  that  the  resulting  formulation  is 
deliberately  cast  in  a  form  particularly  amenable 
to  solution  by  a  digital  computer  or  a  desk  cal¬ 
culator,  emphasizing  the  interplay  between  analysis 
and  computing  machines. 

J.  Documentation 

A  major  portion  of  the  progress  reported  above 
has  been  documented  in  a  comprehensive  US  Army 
Missile  Coraiand  report. 13 

V.  Future  Plans 

Plans  for  the  next  fiscal  year  revolve  primar¬ 
ily  about  impl ementating  Task  3  (define  candidate 
GAC  systms)  of  the  development  program  as  well  as 
completing  any  partially-completed  portions  of  Task 
2  now  underway,  dose  communication  will  be  main¬ 
tained  with  the  intelligence  community  to  become 
aware  of  any  changes  to  the  presently  predicted 
targets  as  contrasted  to  the  attribution  mentioned 
above.  Detailed  analytical  models  of  the  missile 
plant(s),  effectors,  and  autopilot(s)  will  be  devel¬ 
oped.  The  possibility  of  modifying  plant  charac¬ 
teristics  by  making  Innovative  use  of  aerodynamics 
is  to  be  Investigated.  Projected  advances  in  the 
field  of  propulsion  will  also  be  investigated.  The 
comparison  of  various  digital ly-impl  emented  guid¬ 
ance  laws  will  be  continued,  including  those  that 
incorporate  DAC.  Applications  (other  than  to  guid¬ 
ance  laws)  of  DAC  theory  to  improve  system  perform¬ 
ance  will  be  investigated.  Development  of  a  mod¬ 
ular  guidance  simulation  to  implement  this  investi¬ 
gation  has  already  begun.  Trade  studies  concerning 
identified  existing  and  future  sensors  will  be  con¬ 
ducted,  with  the  possibility  of  developing  a  new 
sensor  with  caoabtlities  not  yet  In  existence.  Ad¬ 
vances  in  microcomputer  state-of-the-art  will  be 
watched  closely.  Finally,  tools  to  aid  the  GAC 
system  designer  to  handle  digital  implementation 
will  continue  tc  be  both  assessed  and  developed. 

VI.  Conclusions 

As  indicated  within  this  paper,  there  is  a  cler 
need  for  the  development  of  an  advanced  GAC  system 
for  US  Army  future  tactical  missiles.  This  need  is 
dictated  by  the  predicted  targets  of  the  future, 
the  anticipated  battlefield  of  the  future,  and  the 
characteristics  of  tactical  missiles  that  are  ei¬ 
ther  in  the  Inventory  or  under  development  at  this 
time.  The  nature  of  the  future  GAC  system  will  be 
digital  so  that  the  overall  missile  system  may  be 
availed  of  present  and  predicted  advantages  that 
are  implicit  in  digital  controllers.  The  form  of 


tht  guidance  law  probably  wltl  ba  optimal,  sinca 
tha  parfonaanea  crltarla  that  must  ba  minlmlzad 
must  taka  »ora  Into  account  than  miss  distance. 
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ABSTRACT 

Tha  engagement  of  maneuvering  land  vehicles  with  gun  systems  place  extreme  performance  re¬ 
quirements  on  the  fire  control  system  designs.  The  effectiveness  of  a  gun  fire  control  sys¬ 
tem  depends  on  the  capability  to  provide  an  accurate  fire  control  solution,  i.e.  predict  the 
future  position  of  the  target  a  projectile  tlme-of- flight  later.  Mon  linear  predicticn  is 
shown  to  not  only  improve  performance  but  to  also  increase  available  doe  for  firing  igamse 
aaneuverlng  targets.  Sub-optimal,  multi -variable,  adaptive  estimation  approaches  are  shown 
to  improve  the  effectiveness  of  predictive  fire  control  svstems. 

Sensitivity  analyses  are  presented  that  relate  system  Induced  errors  and  target  motion  in¬ 
duced  errors  to  tracking  noise  and  predictor  order.  Relationships  between  system  stability 
and  performance  for  two  basic  types  of  fire  control  systems  are  presented. 

INTRODUCTION 

This  paper  discusses  the  fire  control  system  problem,  the  nature  of  land  vehicle  mobility  and 
agility  and  the  ability  of  predictive  fire  control  systems  to  effectively  engage  maneuvering 
vehicles-  Existing  performance  specification  do  not  satisfactorily  describe  the  level  of 
maneuverability  expected  in  a  tactical  situation.  Rather,  present  specifications  define  per¬ 
formance  requirements  for  fixed  vehicle  speed  and  heading  movement  which  has  resulted  in  the 
development  of  fire  control  system  designs  chat  are  significantly  degraded  in  i  maneuvering 
target  environment.  The  problem  is  addressed,  in  general,  for  the  four  cases  of  firing 
vehicle- target  vehicle  movements.  The  processes  required  in  the  fire  control  solution  st* 
identified  and  the  sensitivity  of  system  performance  to  the  propagation  of  cracking  errors 
is  discussed.  The  stability  and  performance  characteristics  of  two  generic  fire  control 
system  configurations  are  analyzed  in  some  detail. 

GUN  EIRE  CONTROL  SYSTEM  PROBLEM 

The  purpose  of  gun  fire  control  systems  is  to  have  a  projectile,  that  has  been  fired  a  time 
of  flight  previously.  Impact  the  target  that  was  sighted  a  time  of  flight  earlier.  The 
critical  taction  parameters  that  degrade  the  performance  of  predictive  fire  control  systems 
have  been  identified  as  cyclic  oscillations  exhibiting  frequencies  that  are  within  the 
motion  capabilities  of  tactical  land  vehicles  (1).  Tracking  error,  defined  as  the 
difference  between  target  and  reticle  position,  does  rtoc  in  itself  cause  the  performance 
degradation.  The  Inability  of  the  fire  control  system  to  determine  the  motion  deriv¬ 
ative!  of  the  line-of-eight  (LOS)  to  the  target,  and  predict  the  future  position  of  the 
target  are  the  two  main  factors  that  cause  fire  control  systems  degradation. 

The  error  in  the  ability  of  a  fire  control  system  to  cause  the  projectile  to  Intercept  the 
target  a  time  of  flight  later  Is  referred  to  as  total  gun  pointing  (7GP)  error.  TGP 
error  is  defined  as  the  offset  between  the  actual  gun  pointing  direction  at  round  exit  and 
the  Location  of  the  target  centroid  at  round  impact.  The  TCP  error  is  the  sum  of  the  propa¬ 
gated  system  induced  (SI)  irrors  and  target  induced  (TI)  errors  (i.e.  TSP  error  •  SI  errors 
+  TI  enrors).  The  SI  errors,  considered  in  this  study,  art  the  tracking  error  (difference 
between  cha  tracker  LOS  and  true  LOS  to  the  target)  at  the  tine  of  firing  (lay  error)  and 
the  estimation  errors  (difference  between  estimated  LOS  states  and  true  LOS  states).  The 
SI  errors  propagated  chrough  a  projectile  time  of  flight  result  in  a  kinematic  lead  error. 

The  TI  error  is  caused  by  the  target  motion  during  the  time  of  flight  of  the  projectile. 

It  is  dependent  on  the  order  of  the  lead  solution  in  the  fire  control  system.  For  a  first 
order  lead  system  the  TI  error  is  the  difference  between  the  actual  LOS  movement  during  a 
projectile  time  of  flight  and  the  propagated  LOS  movement  assuming  perfect  LDS  rate  at  the 
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CIO*  of  fir*.  Tha  first  ord«r  predictor  lyttn  TI  «rror  ignore*  ch*  pr*»*oc*  of  actual  cargac 
acceleration  *t  time  of  firing  tad  during  projactlla  flight  tin*.  For  a  aacond  ord«r  load 
system  th*  TI  *rror  is  th*  diffarane*  batvaan  tha  actual  LOS  motion  during  a  projectile  tiaa 
of  flight  and  tha  propagated  LOS  movement  asauaing  parfact  LOS  rate  and  acceleration  at  the 
tiaa  of  firing.  Tha  aacond  ordar  predictor  system  TI  error  accounts  for  target  acceleration 
at  tiaa  of  firing  but  Ignores  the  actual  target  acceleration  during  projectile  flight  elm*. 

This  dlacrlbution  of  errors  is  shown  in  Figure  1.  Th*  ballaitlc  flight  characteristics  of 
tha  projectile  ara  ignored. 

Tha  fir*  control  solution  occurs  during  a  short  tiaa  interval  which  is  related  to  th*  tiaa 
of  flight  of  tha  projectile.  Tha  motion  condition*  of  both  ch*  flrar  and  tha  cargac  ara 
naadad  to  understand  and  aolv*  tha  fir*  control  syitam  problem.  Four  motion  conditions 
exist:  stationary  flrar-atacioaary  cargac,  stationary  flrer-ooving  target,  moving  firar- 
scaclonary  target,  and  moving  firar-moviag  target.  Tha  stationary  firar-statiooary  cargac 
la  tha  Isaac  dynamic  situation  and  is  tha  Isaac  complex  caaa,  and  th*  moving  flrer-moving 
target  1#  ch*  moat  complex  caaa.  For  each  of  tha  eaaaa,  tha  LOS  batvaan  ch*  flrar  and 
target  la  tha  key  to  which  of  tha  four  flra  control  processes  are  being  celled  upon  in  * 
demanding  manner. 

FIRE  CONTROL  STSTSt  PROCESSES 

A  fir*  control  syatem  may  ba  broken  down  into  four  distinct  processes.  Each  of  chase  pro- 
<:•*•••  era  prasanc  in  all  rypas  of  fir*  control  system*.  They  in:  tracking,  aatiaation. 
prediction,  and  gun  pointing.  In  specific  designs  thee*  four  processes  ars  accomplished  In 
different  manner*. 

Th#  cracking  process  Is  important  in  all  four  eaaaa.  For  ch*  moving  flrar  casaa,  tracking 
become*  more  crldcal  because  the  base  motion  of  ch*  flrer  must  ba  comoanaacad  and  it  n av 
b*  effected  in  a  secondary  manner  by  target  motion.  Tracking  is  usually  accomplished  man¬ 
ually  and  is  concerned  with  tha  alignment  of  th*  sight  reciclt  with  th*  target.  Tha  gunner 
la  involved  directly  at  this  stags  and  accuracy  of  cracking  will  be  a  characterization  of 
ch*  ability  of  any  given  gunner  to  perform  tha  task.  Test  data  obtained  from  experimental 
investigations  can  ba  uaad  to  determine  tracking  error  mean*,  standard  davletions.  end 
correlation  time  constants  useful  for  building  models  of  the  cracking  errors. 

Th*  estimation  process  is  ch*  lncaroadlaca  stage  between  the  tracking  process  and  the  pre¬ 
diction  process  and  lei  configuration  is  dependant  upon  tha  ordar  of  th*  prediction  procaaa. 
Estimation  la  th*  process  of  filtering  th*  tracking  data  to  provide  th*  nacaaaary  target 
motion  information  required  in  th*  prediction  procaaa.  Th*  accuracy  of  ch*  tracking  data 
will  influence  th*  performance  of  th*  estimation  proc***.  Th*  syatem  error  Induced  by  Ch* 
estimation  process  dacraaaaa  with  improvement  in  cracking  accuracy. 

Prediction  of  target  future  position  to  obtain  intercept  batvaan  projactlla  and  cargac  Is 
dependant  upon  an  aaclmaca  of  cha  prasanc  motion  of  the  target  and  claa  of  flight  of  the 
projectile.  The  output  of  the  estimator  la  not  a  complete  description  of  ch*  prsaanc  motion 
of  cha  chraac,  therefore,  cha  predictor  does  not  have  the  nacaaaary  Information  to  calculate 
cha  threat's  future  position  exactly.  If  restrictions  are  placed  on  the  allowable  threat 
motion*,  then  ch#  predictor's  ability  to  dstermin*  iti  future  position  is  improved.  Over¬ 
simplification  of  allowable  threat  motions  has  placed  unrealistically  simplified  requirements 
on  ch*  operation  of  ch*  estimation  and  pradlction  proccaaaa.  Realistic  chraac  mocions  ara 
datarmlnad  by  ch*  mobility  capabilicias  of  tactical  vehicles.  In  the  pete,  ch*  majority  of 
threats  Chat  hav*  bean  studied  have  bean  aonaccalaratlng.  Tha  requirements  of  an  ascimator 
and  a  pradlccor  for  chit  type  of  motion  ara  to  combine  ch*  apparent  threat  velocity  estimate 
and  projectile  elm*  of  flight  for  ch*  lead  solution.  The  required  lead  is  constant  and  can 
be  realized  after  some  fettling  claa.  Tha  existence  of  accelerating  cargac*  requires  cha 
estimator  and  predictor  to  develop  constantly  changing  lead  angles,  hanca,  tha  need  for  non¬ 
linear  pradlction. 

An  important  point  to  obsarv*  is  that  for  tha  stationary  flrar-movlng  target  caaa,  tha 
pradlction  procaaa  la  required  to  provide  gun  contend  orders  that  orient  tha  gun  to  account 
for  cargac  motion  during  tha  projactlla' a  time  of  flight,  wharaaa  In  tha  moving  firsr- 
acaclonary  target  caaa  chi*  pradlction  process  is  not  required  because  ch*  LOS  existing 
batvaan  cha  firar  and  cargac  at  Instant  of  firing  do**  not  mova  during  ch*  projactila's 
da*  of  flight.  For  cha  moving  flrar-movlng  target,  LOS  alto  movas  aftac  projactlla  firing. 

Tha  gun  pointing  procaaa  is  required  to  align  and  scabillza  tha  gun  along  ch*  predicted  LOS 
to  tha  target.  Th*  stabilization  and  raaponsa  of  th*  gun  pointing  loop  is  a  major  concern 
for  flra  control  aystam  performance  against  maneuvering  targets.  Stsblllzadon  of  cha  gun 
pointing  procaaa  could  hav*  an  tdvars*  afftet  on  ovsrsll  system  performance.  Th*  moving 
flrar  eaaaa  will  stress  ch*  gun  pointing  procaaa  most  severely  but  It  is  possible  that  tha 
gun  pointing  procaaa  will  b*  squally  straatad  for  th*  stationary  firtr-moving  carget  caaa  with 
non  linear  prediction. 
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FIRE  CONTROL  SYSTEM  CONFIGURATIONS 

The  three  currently  uaad  firs  control  conflgursclons  ara  known  as  manual  or  iron  tight*  dis¬ 
turbed  reticle  and  atablllzad  sight-dirsctor  aystana.  A  fourth  mehtod  callad  closed  loop 
rafar*  to  projectile  spotting  to  ad jute  tha  fira  control  aolutlou  and  la  not  conaldarad  In 
Chla  dlacuaalon.  Tha  oanual  fira  control  ayttan  uaaa  tha  bmta  force  approach  and  concantratea 
on  atablllzlng  tha  gun  poaition  exclusively.  In  this  syscan  tha  lead  la  introduced  manually, 
therefore,  chara  la  no  autoaatlon  of  tha  fira  conf.  1  aatlnatlon  and  prediction  process.  Tha 
disturbed  reticle  systsa  stabilizes  tha  gun  position  and  disturbs  tha  position  of  tha  track¬ 
ing  reticle  from  the  gun  line  position.  In  this  scheme  tha  tracking,  estimation  and  predic¬ 
tion  processes  ara  inseparable  and  tha  fira  control  solution  Is  automated.  Tha  rejection  of 
flrar  vehicle  base  notion  Is  difficult  to  accomplish  in  this  type  of  system.  Tha  last  system 
to  be  considered  is  tha  stabilized  sight -director  system.  Tha  tracking  process  is  accomplish¬ 
ed  by  a  tracker  which  is  rotated  from  a  stabilized  base  chat  has  signals  applied  to  isolate 
firar  vehicle  base  movement.  The  resulting  LOS  orientation  is  referenced  to  incarlal  space, 
as  contrasted  to  tha  gun  line  for  disturbed  reticle  systema.  Tha  estimation  process  is  the 
incarnadines  link  between  tha  cracking  process  and  the  prediction  process.  The  prediction 
process  uses  tha  estimation  process  outputs  combined  with  projectile  time  of  flight  to  deter¬ 
mine  tha  gun  pointing  cossiands.  Tha  gun  pointing  process  uses  the  estimated  LOS  co  the  target 
summed  with  tha  calculated  lead  offset  of  the  gun  from  the  tracker  LOS  to  position  the  gun 
line. 

How  wall  a  fire  control  syacem  configuration  performs  is  a  function  of  target  movement, 
flrar  movement  and  fira  control  system  design-  Tha  analytical  mathodology  required  to  study 
this  problem  should  be  constrained  to  real  time  solution  aechanlsms.  Another  way  to  say  this 
la:  post  dare  analysis  techniques  using  data  obtained  from  field  tests  will  cot  provide  the 
insight  chat  is  required  to  obtain  an  understanding  of  the  relative  performance  of  different 
fira  control  systems.  Probability  of  hit  information  is  useful  for  an  assessment  of  systems 
chat  have  bean  fialdad  but  is  not  applicable  for  crade  off  studied  of  the  type  required  in 
this  SCudv.  Analytical  methodologies  such  as  servo  mechanism  synthesis  and  modern  filtering 
technology  ara  required  to  study  this  problem. 

MANEUVERING  TARGET  DESCRIPTION 

A  quantitative  description  of  the  threat  is  required  to  evaluate  cha  performance  of  fire 
control  systema  operating  against  maneuvering  targets.  To  develop  this  description,  it  is 
necessary  to  consider  tha  Tsohility  and  agility  characteristics  of  threat  vehicles  in  a  real¬ 
istic  combat  environment.  A  thorough  description  of  anticipated  maneuvering  seema  to  defy 
identification  because  threat  maneuvers  constitute  a  large  set  of  possibilities  even  when 
constrained  by  tactical  doctrine,  driver  policy,  terrain  and  vehicle  capabilities.  Two 
approaches,  analytical  and  empirical,  are  available  for  consideration  in  the  attempt  to 
identify  the  maneuver  characteristics  of  land  vehicles.  An  analytic  approach  would  view  each 
maneuver  aa  being  composed  of  elements  from  an  idealized  group  of  movements.  An  empirical 
approach  would  view  the  maneuvers  as  having  actually  occurred  during  limited  tests  for  diff¬ 
erent  types  of  maneuvering  vehicles.  Neither  of  these  approaches  provide  a  comolete  maneuver 
description,  but  a  combination  of  these  two  approaches  offers  some  advantages  and  Is  the 
rationale  adopted.  The  analytic  approach  will  partially  overcome  the  incompleteness  of  the 
empirical  daca  base  while  the  empirical  data  will  offset  the  mathematical  idealizations  of 
Cha  analytic  mathodology. 

Empirical  Approach 

When  using  empirical  daca  to  demonstrate  the  performance  of  a  gun  fire  control  system,  base¬ 
line  performance  can  be  determined  with  no  concerns  arising  from  idealization  of  the  man¬ 
euvers.  Sines  tha  number  of  maneuvers  will  bs  rather  small,  chav  neither  provide  sufficient 
Information  about  cha  robusesasa  of  a  fira  control  design  mehcodology  nor  the  pathology  whan 
cha  fira  control  ayscam  begins  to  dtgrsde.  When  demonstrating  tha  performance  of  a  fira 
control  ayacaa  againac  experimental  data,  caution  must  be  exercised  co  asaurc  that  tha 
aaqilrlcal  daca  is  properly  inputted  co  tha  fira  control  svstes  modal.  Matching  of  tha  data 
racaa  and  noise  levels  often  requires  soma  preprocessing  of  exparlaantal  data  co  prepare  it 
for  use  in  simulation  studies. 

Analytic  Approach 

Aa  a  supplemant  co  tha  empirical  approach  tha  analvtlc  approach  is  uaad  to  invastlgac*  sen¬ 
sitivity  effects  for  a  larger  group  of  movements.  Simulating  nev  or  pathological  maneuvers 
require  that  cha  analytic  capability  superimpose  maneuvers  arising  from  random  disturbances 
ar.d  Intentional,  voluntary  vehicle  driver  coraands.  The  random  disturbances  may  be  represented 
in  terms  of  time  histories  or  power  spectral  densities.  The  time  history  approach  is  basac 
on  tha  development  of  a  mathematical  model  of  vehicle  movement  influenced  bv  terrain  effects 
and  arbitrary  driving  habits  of  individual  drivers .  It  is  assumed  that  for  no  random  effects 
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caused  by  terrain  irregularities  or  driver  input,  the  vehicle  would  follow  a  straight  line- 
conatent  speed  peth.  Maneuvers  are  viewed  as  perturbations  on  this  straight  line-constant 
spaed  path.  Apparent  acceleration,  a  (t),  which  la  correlated  in  time,  accounts  for  the 
vehicle's  deviation  froa  a  atralght  line  peth.  Maneuver  capability  is  expressed  by  three 
quantities:  the  variance,  or  aagnltuda  of  a  (t),  the  cyclic  naneuver  frequency  and  the  else 
constant  of  the  maneuver . 

Intentional,  voluntary  vehicle  driver  commanded  notion  of  land  vehicles  over  terrain  is  a 
complicated  subject  in  Itself  and  will  not  be  investigated  in  this  study.  It  is  recognized 
however,  that  an  interaction  between  vehicle  horsepower,  weight,  suspension,  and  locoaotion 
concepts  do  coabina  with  terrain  over  which  It  is  aovlng  to  provide  different  levels  of 
nobility  with  respect  to  a  fixed  reference  fraae.  Therefore,  different  vehicle  designs  will 
have  different  nobility  levels  defined  in  tens  of  notion  and  derivatives  of  notion.  Agil¬ 
ity  is  closely  related  to  nobility  end  yet  it  is  a  slightly  different  description  of  inten¬ 
tional  vehicle  notion.  Where  nobility  describes  the  covenant  of  a  vehicle  froa  one  location 
to  another  location  in  a  given  period  of  tine,  agility  describes  the  vehicle's  ability  to 
alter  Its  naan  path. 

SENS ITIV ITT  OF  FISE  CONTSOL  PROCESSES 

Degradation  in  gun  pointing  accuracy  results  froa  two  najor  error  sources,  systen  and 
c*rget  Induced  errors.  The  target  induced  errors  are  caused  by  the  notion  of  the  target 
during  the  tlaa-of-fllght  of  the  projectile.  Since  the  target  haa  the  capability  to  man¬ 
euver  within  constraints  of  the  terrain,  vehicle  characteristics  and  driver  policy  during  a 
projectile's  tlae-of-f light,  there  Is  no  such  thing  as  a  correct  (perfect)  lead  solution. 

The  lead  solution  is  based  on  the  projected  target  position  using  the  present  target  states 
and  projectile  tlae-of-f light.  Therefore,  the  target  induced  error,  in  general,  cannot  be 
reduced  to  zero  for  a  maneuvering  target.  However,  it  can  easily  be  shown  that  the  pre¬ 
diction  procasa  is  capable  of  reducing  the  gun  pointing  error  due  to  target  notion. 

The  systaa  induced  errors  ere  made  up  of  bias  and  randoa  errors  emanating  from  specific 
components  and  subsysteas.  The  propagation  of  these  errors  degrade  the  performance  of  the 
fire  control  systaa.  The  systea  induced  errors  of  major  concern  are  chose  occurring  in  the 
cracking  process.  Sensitivity  analyses  have  been  performed  to  evaluate  the  degradation  of 
gun  pointing  coaaenda  to  tracking  process  errors. 

The  analysis  considers  the  fire  control  systaa  processes  to  be  interfaced  in  tandem  wlch 
no  feedback  of  outputs  to  a  previous  process.  The  analysis  is  further  lialted  to  a  segment 
of  a  maneuvering  target  path  which  was  generated  by  a  maneuvering  target  path  simulation 
prograa.  This  analytically  generated  peth  provides  an  exact  tine  history  of  the  target  states 
(position,  velocity  and  acceleration). 

The  cracking  process  is  modelled  by  ousting  randoa  errors  of  known  variance  wlch  the  output 
of  a  perfect  LOS  sensor.  The  output  of  the  tracking  process  is,  by  definition,  in  IDS 
coordinates,  however,  cartesian  coordinates  ere  used,  by  choice  and  not  a  limitation  of  the 
methodology,  in  the  estimation  processes.  For  s implication,  the  transformation  from  LOS  to 
cartesian  coordinates  is  accomplished  prior  to  adding  cracking  noise. 

A  sub-optimal,  adaptive  Kalman  filter  (K7)  is  used  for  the  estimation  process  in  the  generic 
fire  control  system  under  consideration.  The  noisy  tracking  process  signal  is  processed  by 
the  KF  to  provide  e  "best"  estimate  of  the  target  states  (poslclon,  velocity  and  acceleration) . 
The  aetlmatlon  errors  are  minimized  by  providing  the  filter  with  the  correct  variance  of  the 
observation  noise.  In  practice,  this  perfect  match  of  noise  variance  is  not  achleveable  but 
can  be  approached  with  detailed  error  analysis  of  the  cracking  process  or  with  software 
methodology  to  estimate  the  noise.  The  latter  is  probably  desirable  and  necessary  because 
Che  variance  of  the  tracking  process  error  is  not  cine  invariant  in  a  combat  environment. 

The  KT  Is  the  generic  fire  control  system.  The  KF  equations  end  theory  are  wall  known  and 
era  presented  elsewhere  (2,J).  However,  the  adaptive  feature  of  the  designed  KF,  which 
requires  on-line  computation  of  the  filter's  gain,  is  outlined  (A).  The  adaptive,  time 
varying  gain  is  obtained  by  changing  the  variance  of  the  uncertainty  of  the  embedded  target 
dynamics,  as  a  function  of  the  estimated  path  geometry.  The  forcing  fimecion  for  the  target 
dynamics  is  modelled  as  a  random  (Gaussian  noise)  race  of  change  of  acceleration.  The 
varlenca  of  u  is  defined  in  the  body  coordinates  of  the  target  as  constant,  diagonal  elements 
of  the  q  matrix.  The  Q  matrix  is  rotated  as  the  target  maneuvers  to  provide  a  time  varying 
Q  matrix  in  the  filter's  coordinate  systea. 

The  sensitivity  of  the  estimates  to  the  tracking  process  noise  is  evaluated  for  a  typical 
maneuvering  target  path.  The  ground  crack  of  the  maneuver  is  shown  in  Figure  I.  The  max¬ 
imum  speed  and  lateral  acceleration  are  10  m/sac  and  2  a/ sec*,  respectively.  Figure  3  shove 
the  degradation  in  velocity  estimates  as  the  standard  deviation  of  the  cracking  process  noise 
on  the  assumed  position  observation  is  increased  from  O.OS  mater  to  1.0  mater.  The  degrade- 


cion  is  the  estimate,  of  lateral  acceleration  for  cha  same  nolaa  lavala  1*  shown  Is  Figure  4. 
A  comparison  of  chaaa  evo  figures  shows  chat  cha  valoclty  estiaatts  art  not  as  senslcivt  to 
cha  propagation  of  cracking  nolsa  as  cha  acceleration  eseiaaeas. 

Ths  prediction  process  provides  cha  command  for  pointing  cha  gun  co  cha  predicted  target 
position.  The  estimated  future  position  of  cha  target  depends  on  the  order  of  the  prediction 
process.  Ideally,  ona  would  Ilka  co  foracssc  cha  target  position  so  chat  a  projectile  fired 
a  clna-of-f light  earlier  would  arrive  at  a  point  In  space  slsulcaneoualy  with  the  target. 
Unfortunately,  only  the  presaac  states,  which  are  never  known  exactly,  are  available  for  use 
la  computing  future  cargec  position. 

With  knowledge  of  cha  true  future  position  of  the  target,  available  from  the  target  motion 
modal,  the  degradation  in  the  gun  pointing  commands  can  be  evaluated  for  different  tracking 
errors.  Target  lnducad  errors  and  the  propagation  of  the  tracking  process  noise  are  analysed 
to  evaluate  chair  effect  on  gun  pointing  commands . 

The  carpet  induced  errors  are  functions  of  target  maneuver  characteristics,  projectile 
else  of  flight  end  prediction  order.  For  a  given  prediction  order  and  with  perfect  know¬ 
ledge  of  cha  present  target  states  and  time  of  flight,  the  resulting  target  Induced  errors 
are  lower  bound  prediction  errors.  Effects  of  time  of  flight  and  order  of  prediction  are 
shown  in  Figure  5  for  a  maneuvering  cargec  whose  maxima  speed  and  lateral  acceleration  is 

10  m/sec  and  3.S  m/sec'.  Prediction  errors  are  improved  for  decreases  in  time-of-f lighc  and 
higher  order  of  prediction. 

Tlrse  order  prediction  is  linear  and  requires  only  accurate  estimates  of  velocity  to  approach 
the  lower  bounds  of  prediction  error.  Second  order  prediction  requires  not  only  accurate 
velocity  but  also  acceleration  estimates  to  minimise  the  prediction  errors.  Figure  6  shows 
the  standard  deviation  of  prediction  error  for  the  target  maneuver  shown  in  Fig.  2  as  a  func¬ 
tion  of  tine-of- flight  and  variances  of  tracking  process  noise  for  first  order  prediction, 
‘h***  results  indicate  chat  the  degradation  in  prediction  error  is  minimized  as  the  quality 
or"  backing  improves.  However,  the  existence  of  the  lower  bound  curve  for  second  order  pre¬ 
diction  provides  additional  improvement,  not  realised  by  first  order  prediction.  Assuming 
position  observations  (input  co  Che  KF>  with  e  lo  noise  of  1.0  meter.  Figures  6  and  '  show 
that  chera  is  no  large  difference  between  first  and  second  order  prediction.  However,  second 
order  prediction  with  a  reduction  in  the  tracking  process  error  to  0.05  meter  (15  microradians 
at  a  range  of  2000  maters)  provides  a  significant  improvement  In  the  lead  solution.  Unlike 
f-rst  order  prediction,  second  order  prediction  is  not  only  more  sensitive  co  the  tracking 
process  noise  but  also  to  the  observation  stata.  Figure  7  shows  that  Improvements  are  re- 
/  *  Ch*  obMrv*«°n»  rates  rather  then  position.  If  tracking  accuracies  of  0.04 

m/sec  (20  microradian/ sec  at  2000  meters)  ace  achieved,  the  prediction  erTor  is  within  about 
t«n  percent  of  the  lower  bound  for  second  order  prediction. 

The  lead  errors  discussed  above  era  the  differences  between  the  predicted  and  actual  target 
positions  for  an  estimated  clna-of-f light.  Targets  are  not  point  sources  and  a  more  mean- 
-ngful  criteria  for  evaluating  tha  syscam  processes  is  the  percent  time  on  target  for  a 
spec it -ed  angagamanc  time.  Assuming  a  targat  size  of  2.3  maters  X  2.3  meters,  independent  of 
target  orientation,  the  percent  time  on  target  for  the  same  tracking  accuracy  in  Figure  '  is 
depicted  in  Figure  a  for  times  of  flight  between  1.0  seconds  and  2.5  seconds. 

STABILITT  ANALYSIS  OF  GENERIC  FIRE  CONTROL  SYSTEMS 

C«n«r*i  Discussion 

configurations  in  exi.cence:  manual,  disturbed  reticle  and 
in  tarae  of  how  the  fire  control  pro- 

e!!“?aa!!.??fhf 1**!  J11  Operational  systems  utilize  che  human  operator  to  null 

^?Aitf;r*aC»b*tV~°  observed  c.rg.t  and  th.  reticle  position.  The  d.gra.  of  p.rtici- 
di*fer«nf  ^  eseh  of  the  three  types  of  fire  control  systems  Is  considerably 

=aMU««ion  ?n  stability  of  the  closed  loop  man-mmehin.  system  is  an  Important 

consideration  in  determining  performance  end  is  one  of  the  primary  dlstinguishlna  feature, 
that  ch.racc.rlt.,  th.  tff.cclv.ns,.  of  the  three  tvpa.  of  fir.  fontr"  s^tlmi?  I**the** 
th*  t^‘Cl‘ln*•  ®  s timet ion  end  prediction  processes,  sr.  D.rformad  bv  tha  man 
"dad  bv^  rJI^I  !?' Y  °rlrC  th*  lln*  ln  ‘“ordanc.  with  th.  information  pro- 
Urtt-dir^r  b^  ch«  “n  K>  th.  disturbed  reticle  and  stabilized 

ayatema,  hovev.r  it  is  accomplished  differently.  The  estimation  and  predl-- 
tion  processes  slso  mac  halted  differently  In  th...  two  tvpa.  of  f  irilontrolsyst*., 
^dl«urb!d^Scn?nt  ‘“J""?*  •dvent.ga,  of  «  stabilized  sight-director  system  compared  to 
prediction  nrL  r?  decoupling  of  th.  tricking  proc.ss  from  th.  estimation  and 

r«i 2‘.  T  V  tUrr,nC  *nd  s«tv.  at  th.  r.f.renc.  from  which  the 

rhl^t-!  dn  th*  il'curb*d  "title  system.  Involvement  of  th.  humen  gunner  In 

the  .urret  loop  for  the  dlstrubed  reticle  system  end  hi,  absence  from  the  turret  loop  for  th. 
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stabilized  sight-director  system  Is  s  distinguishing  feature  of  the  systems.  The  cracking 
process  is  therefore  aore  isolated  from  the  estimation,  prediction,  and  gun  pointing  processes 
in  the  stabilized  sight-director  system. 

Disturbed  Reticle  Fire  Control  System 

One  cite  control  configuration  in  current  use  is  the  disturbed  reticle  concept.  The 
following  discussion  is  lntsnded  to  describe  in  detail  the  functioning  or  the  disturbed 
recicle  fire  control  system  and  identify  the  four  processes,  showing  how  each  is  related  to 
the  other.  Figure  9  describes  the  signal  flow  in  such  a  system  and  the  four  major  processes 
have  been  identified  in  terns  of  where  in  the  system  each  is  accomplished. 

The  input  to  the  system  is  the  LOS  from  the  target  to  the  reticle  of  the  tracking  system. 

The  human  operator,  present  in  most  current  systems,  moves  the  handle  bar  controller  to 
align  the  reticle  of  the  tracking  system  to  be  coincident  with  the  target.  The  ability  of 
any  human  controller  to  accomplish  this  task  defines  the  quality  of  the  tracking  process. 
Handle  bar  controller  output,  which  is  dirsctly  related  to  the  LOS  rate,  is  ueed  to  drive  two 
independent  subsyateoe.  The  first  is  the  turret  servo  which  is  contended  to  rotate  at  a  rate 
directly  proportional  to  the  handle  bar  controller  deflection.  The  second  subsystem  driven 
by  the  handle  bar  controller  is  a  lead  screw  3arvo  and  reticle  system.  The  displacement  of 
the  lead  screw  servo  is  directly  proportional  to  the  filtered  handle  bar  controller  deflection 
multiplied  by  the  projectile  time  of  flight.  The  lead  screw  displacement  is  used  to  position 
tha  recicle  of  the  tracking  system. 

These  are  Cwo  distinct  feedback  signal  paths  in  the  disturbed  reticle  conf igureticn  and  the 
human  le  a  series  subsystem  in  both  paths.  Another  important  observation  is  to  note  that 
the  signal  loop  made  by  the  turret  servo-man-'nandle  bar  controller  is  a  degenerative  feedback 
loop  because  of  the  negative  summing  Junction.  The  signal  loop  made  by  the  filter-time  of 
flight  lead  servo-reticle  servo-man  handle  bar  controller  is  a  regenerative  feedback  loop 
because  of  cwo  negative  surra  in g  junctions.  During  norael  operation  of  the  disturbea  reticle 
system,  Che  performance  of  these  cwo  feedback  paths  give  rise  to  a  dynamical  system  chat 
exhibits  soma  undesirable  performance  characteristics.  Without  further  crossfeed  compensa¬ 
tion.  the  closed  loop  performance  of  the  disturbed  reticle  system  is  at  best  marginally  stable 
and  at  worst  unstable.  To  overcome  this  condition,  compensation  signal  oaths  are  added.  The 
basic  compensation  is  a  tachomacer  generator  signal  from  the  lead  screw  servo  which  is  com¬ 
bined  with  Che  turret  servo  error  signal.  This  composite  signal  is  fed  to  the  turret  servo 
and  the  recicle  servo  to  compensate  for  cha  dynamical  mismatch  that  occurs  In  the  recicle  ar.d 
turret  servos.  However,  there  is  no  such  thing  as  a  perfect  compensation  and  the  undesirable 
performance  characteristics  alluded  to  earlier  can  never  be  complete!-'  nullified,  not  to  men¬ 
tion  the  potentially  precarious  situation  that  might  occur  if  any  failure  or  shift  occurs  in 
tha  compensation  paths. 

Ths  important  thing  to  observe  about  the  root  locations  in  the  figures  10  and  11  is  that 
there  are  numerator  roots  in  the  righc  half  of  the  S  plant.  This  arises  from  the  basic 
disturbed  reticle  configuration  and  muse  be  considered  a  fixed  element  phenomenon  in  this 
type  of  system.  The  poles  or  denominator  roots  describe  che  system  operating  polnc  for  a 
system  gain  of  zero.  The  zeros  or  numerator  roots  describe  the  svseem  operating  polnc  for  a 
system  gain  of  infinity.  The  dotted  trajectories  connecting  these  two  extremes  ate  a  pictor¬ 
ial  description  of  che  operating  point  loci  for  all  intermediate  gains.  These  systems  exhibit 
conditional  stability  because  of  the  presence  of  positive  feedback  in  the  equivalent  transfer 
function  between  B  and  A.  These  are  different  closures  than  exist  for  a  negative  feedback 
chat  occurs  when  both  the  reticle  and  turret  crossfeeds  are  present  as  shown  in  Figure  IZ. 

The  existence  of  these  simultaneous  crossfeeds  from  che  lead  screw  servo  ar.d  turret  servo 
error  to  the  turret  servo  and  reticle  servo  tend  to  offset  che  non-minimum  phase  root  con¬ 
dition  shown  in  Figures  10  and  11. 

In  summery,  it  is  the  location  of  chi  oparsting  points  that  decerslnt  the  system  stability 
characteristics.  The  frequency  content  of  the  cracking  error  is  directly  releced  to  the 
operating  points,  but  squally  Important  is  the  magnitude  of  the  tracking  error  which  is  in¬ 
fluenced  by  cha  location  of  che  numerator  roots  of  che  closed  loop  transfer  function.  These 
effect*  are  interrelated,  but  the  fundamental  underlying  requirement  is  to  achieve  an  adequate 
stability  margin  of  che  closed  loop  system.  This  stability  consideration  is  important  for 
fire  control  system  performance  and  the  designers  must  cake  these  fsetors  into  account.  The 
end  result  Is  system  performance  which  may  be  acceptable  or  not  acceptable. 

It  can  be  asked  why  so  such  concara  abouc  this  situation  because  disturbed  reticle  svstems 
have  performed  satisfactorily  in  ths  past.  Psrhsps  this  is  so,  but  with  the  introduction 
of  maneuvering  targets,  cha  performance  of  this  type  of  svstea  aav  be  adversely  effected. 

When  the  target  LOS,  shown  in  Figure  9  moves  at  a  constant  race,  che  .-lumen  operator  is 
required  to  move  the  handlebar  controller  a  nominal  fixed  amount.  The  turret  servo  develope 
a  fixed  nominal  rate  and  tha  Lead  screw  servo  assumes  a  fixed  nominal  position.  It  then 


become*  the  task  of  eh*  human  co  perturbat*  eha  hand la  bar  conerollar  about  this  nominal 
position  la  ordar  eo  mlnlalzs  tha  tracking  arror.  When  tha  targat  LOS  rat*  la  not  constant, 
which  la  tha  situation  for  msnauvaring  targats,  tha  hand la  bar  controllar  must  b*  moved 
consistent  with  tha  changing  targat  LOS  rata.  Tha  nominal  handle  bar  controllar  position 
la  not  tha  only  dlffaranca  In  eha  sratarn  oparatlon  for  manauvarlng  targats.  Tha  turret 
servo  accelerates  and  decelerates  and  tha  lead  screw  servo  Is  constantly  being  driven  to  a 
new  position.  Tha  position  of  eha  raelda  Is  a  result  of  thas*  two  signals  paths.  The 
dynamic  performance  mismatches  era  quarantaad  eo  ba  greater  than  for  tha  non-manauvarlng 
targats  and  tha  tracking  performance  will  ba  degraded.  This  degradation  occurs  from  the 
Inability  of  the  closed  loop  system  eo  accurately  null  tha  constantly  changing  target  LOS 
rata.  Tha  extant  of  this  degradation  nay  not  b*  1  Mediately  obvious  to  eha  casual  analyse, 
but  eha  oscillatory  nature  of  this  degradation  will  ba  observed  one*  a  sufficiently  close 
survey  of  eha  tracking  arror  la  made.  It  la  Imperative  that  tha  resulting  stability  aargln 
of  tha  closed  man-aaehlna  system  b*  large  to  Insure  acceptable  performance  against  maneuv¬ 
er  lag  targats. 

Recant  work  has  shown  that  tactical  esrgaea  can  execute  maneuvers  of  such  a  nature  that 
when  projectile  elmas  of  flight  of  l.S-2.0  sac  era  considered,  target  Induced  notion  after 
projactll*  firing  will  causa  excessive  miss  distances  when  linear  predictor  fir*  control 
systems  era  assumed  and  more  over  thas*  miss  distances  can  be  significantly  reduced  when  non¬ 
linear  or  higher  ordar  predictor  fir*  control  systems  ara  employed.  Thas*  observations  Indi¬ 
cate  lower  boundary  alas  distances  ara  possible  for  non  linear  lead  systems.  When  this  sit¬ 
uation  Is  presented  to  the  fire  control  designer,  his  inclination  will  be  to  consider  the 
possibility  of  Including  nonlinaar  prediction  in  tha  fir*  control  system.  In  the  disturbed 
raticl*  configuration  shown  In  Figure  9,  this  may  b*  a  design  Impossibility  because  of  the 
level  of  cracking  performance  obtainable  from  tha  oparatlon  of  the  disturbed  reticle  systems. 
To  be  more  specific,  the  tracking  error  to  develop  usable  target  accelerations  for  nonlinear 
prediction  Is  required  to  ba  smaller  than  tha  tracking  error  for  first  order  prediction. 

Tha  trade-off  between  tha  propagated  system  induced  tracking  errors  for  the  nonlinear  esti¬ 
mation  process  must  ba  offset  by  the  target  induced  prediction  arror  luprovementa  realized 
by  the  higher  order  prediction.  The  key  Ingredient  for  this  situation  to  exist  In  a  fire 
control  system  Is  to  have  high  quality  tracking  errors. 

If  the  human  tracker  Is  replaced  by  an  automatic  tracker  tha  performance  limitations  imposed 
by  tha  loop  structures  In  a  disturbed  reticle  system  may  negate  the  potential  improvement 
attainable  from  the  improved  tracking.  It  is  the  coupled  nature  of  the  tracking,  estimation, 
prediction,  and  stabilization  process  occurring  in  the  disturbed  reticle  configuration  that 
restrict  Its  growth  to  better  fire  control  syetam  performance,  especially  agelnsc  maneuvering 
targets. 

Stabilized  Sight-Director  Fir*  Control  System 

A  stabilized  sight-director  fire  control  system,  shown  in  Figure  13  is  actually  tvo  distinct 
system*  that  are  brought  together  to  accomplish  the  trecklng,  estimation  and  prediction  pro- 
ceaees  of  a  fir*  control  system.  Stabilization  of  tha  tracking  system  Is  independent  from 
stabilization  of  the  turret.  The  stabilized  sight  Is  decoupled  from  turret  and  hull  motion 
by  the  reverse  corquing  of  the  outer  glmbal  of  the  cracker  to  account, for  distrubances  of 
the  tracker  base  which  is  mounted  on  Che  turret.  This  decoupling  enhances  the  ability  of  the 
cracker  eo  maintain  coincidence  between  the  sight  reticle  and  the  target  LOS.  The  stabilized 
reticle  position  can  utilize  both  position  and  rat*  feedback  to  augment  the  stability  of  the 
sight.  The  orientation  of  tha  sight  reticle  is  therefore  an  independent  process  from  the 
turret  motion. 

Position  and  race  of  the  LOS  art  fed  to  a  filter  or  ettlnatlon  process  to  determine  the 
nacettary  Information  about  the  LOS  to  the  target  that  will  be  needed  to  offset  the  turret 
servo  from  the  stabilized  tracker.  Multi-variable,  sub-optimal  technology  can  be  applied  to 
further  Improve  the  quality  of  trecklng  that  can  be  realized  from  the  stabilized  sight- 
tracker.  Therefor*  either  linear  or  non  linear  estimates  of  LOS  movement  can  ba  considered 
aa  poaalhilltlea.  If  LOS  accelerations  sr*  to  be  estimated,  the  appropriate  modeling  of 
target  dynamics  and  tracker  uncertainties  will  be  required  to  insure  chat  the  degree  of 
suboptlmellty  is  not  excessive.  One  very  significant  plus  that  couples  ths  estimation  and 
trecklng  process  in  s  favorable  manner  is  the  utilization  of  sight  line  rate  siding  feedback 
to  the  tracker  obealned  from  estimation  of  the  target  rates  end  acceleration.  This  concept 
relax**  the  teak  of  the  human  tracker  or  auto- cracker  and  will  improve  the  minimization  of 
tracking  error. 

Output  of  tha  target  state  estimator  Is  used  in  tvo  separate  oaths.  Tha  first  path  uses 

5^  and  to  drive  the  turret  servo  as  a  director  to  follow  the  trecker  LOS.  The  second 
signal  path  combine*  target  state  estimates  with  projectile  elm*  of  flight  and  offsets 
the  gun  from  the  cracker  LOS  by  the  appropriate  velue  to  permit  Intercept  of  projectile  and 
targat  a  time  of  flight  later. 
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Parformanca  of  cho  stabllixad  eight-director  system  should  not  bo  compromised  by  maneuvering 
targets  to  ch«  «xtoat  chat  cha  disturbed  ratlclo  oyataa  is  compromised.  The  baalc  reaaon 
far  eh  la  la  that  Cha  cracking  system  la  aaaaatlaUy  decoupled  from  the  lead  prediction  aye- 
tea.  However,  there  are  aome  inherent  atabllizatlon  probleaa  chat  can  occur  In  thia  config¬ 
uration  and  they  are  accentuated  by  the  ceapeation  Co  obtain  high  praforaanee  of  the  gun 
pointing  procaaa.  The  argument  goea  aa  follove:  with  Increased  cracker  perforaaace,  the 
gun  atabllixatlon  servo  can  be  aada  to  perfora  more  rapidly,  thereby  Increasing  the  overall 
capability  of  the  ayataa.  However,  with  increased  performance  being  required  of  the  turret 
servo  to  follow  the  turret  coaaands,  the  stability  of  the  turret  servo  may  be  compromised 
because  of  the  high  gains  in  the  director- follower  loop.  Experience  with  slallar  types  of 
systaaa  has  shown  thac  because  of  non-rlgld  gun  cube  and  hull  structures,  Che  follower  loop 
systaa  must  be  phase  stabilized  and  not  gain  stabilized,  as  Is  the  case  for  less  responsive 
systaSM  such  as  disturbed  reticle  ayateas.  This  requires  sophisticated  cementation 
circuits  to  overcome  system  instabilities. 

The  stabilized  sight  is  identified  between  cha  target  input  and  the  sight  output  is  Figure  13. 
It  in  turn  drives  the  gun  turret  servos  which  are  used  to  position  the  base  of  the  stabilized 
sight.  The  signal  flow  diagraa  and  root  loci  for  the  stabilized  sight-director  systea  is 
shown  In  Figure  14.  This  la  the  saaa  basic  root  locus  obtained  in  the  disturbed  reticle 
system  whan  the  cross  fseds  ware  included.  The  dotted  lines  show  the  movement  of  the 
stability  as  the  gain  is  Increased.  The  addition  of  series  compensation  circuits  In  the 
tracker  transfer  function:  such  as  T1  5+1  *  whlch  c“  b*  »dd»d  in  «  straight-forward 

T,  S+I 

manner  will  alter  the  shape  of  the  loci  to  obtain  an  optimized  operating  point,  which  would 
be  difficult  in  the  disturbed  reticle  systea.  The  fundamental  purpose  of  che  tracking  pro¬ 
cess  is  to  align  with  Simultaneously  any  disturbances  on  che  stabilized  sight  are 

compensated  by  orientation ‘of  the  sight  base  thereby  simplifying  the  tracking  teak. 

CONCLUSIONS 

The  inherent  ability  of  a  stabilized  sight-director  fire  control  system  to  decouple  the  track¬ 
ing  estimation,  prediction  and  gun  pointing  proedsses  may  ba  exploited  to  Improve  effective¬ 
ness  whan  angaglng  maneuvering  targets.  Accurate  tracking  is  necessary  for  non  linear 
prediction  and  tail ti-var table,  sub-optimal  design  tschnology  is  required  to  schieve  the 
needed  accuracy  of  the  target  state  estimates  for  mechanizing  nonlinear  pradlctlon.  Further 
studies  are  required  to  identify  che  specific  detail!  of  the  resulting  systea  design.  A 
complementary  methodology  employing  stability  and  performance  analyses  will  assist  in  this 
quest. 

ACKNOWLEDGEMENT 

The  authors  thank  Mr.  Warren  Muehlberger  for  drawing  the  figures  and  Ms.  Sharon  Taylor  for 
typing  the  manuscript. 


REFERENCES 


1.  Burke,  Harold,  "Fire  Control  Svstea  Performance  Degradation  Vhan  A  Tank  Gun  Engages  A 
Maneuvering  Threat"  XVII  US  Army  Operations  Research  Symposium,  Ft.  Lea,  VA  i-9  Nov  19'8. 

L**t!,rua»  James  F.,  "An  Approach  To  Firs  Control  Svstam  Computations  end  Simulations". 
US  Army  Materiel  Systems  Analysis  Activity  Technical  Report  No.  126,  April  1975. 

3.  Gelb,  Arthur.  "Applied  Optimal  Estimation".  MIT  Press,  1974. 

4.  Burke,  Harold  H.,  Perkins,  Tony  R. ,  Uathrua,  James  F . ,  "Stats  Estimation  of  Maneuvering 

Zpi&il**  jLk*  Jllcgrlag" .  US  Army  Materiel  Systems  Analysis  Activity  Technical  Rsport 

No.  186,  October  1976. 


♦chi-  ;  •  ;  . 

IW .  ■* 


16 


VtlCCU 


Figure  9  Disturbed  Reticle  Fire  Control  System 
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Fi&*a  11.  Disturbed  Reticle  Root  Locus, 
Turret  Servo  Cross  feed. 


Figure  12.  Disturbed  Reticle  Locus  with  Turret 
Servo  and  Reticle  Goss  feeds 
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Figve  13  Stabilized  Sight  -  Director  Ti re  Control  System 


Figure  14  Stabilized  Sight  -  Director  Root  Locus. 
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DETECTION  AND  INITIATION  OF  FIRING  COMMAND  FOR  AN  ACCELERATION  PREDICTOR 
FIRE  CONTROL  SYSTEM  ENGAGING  MANEUVERING  TARGETS 


H.  H.  Burke  -  AMSAA 


Maneuvering  targets  degrade  the  performance  of  all  state  of  the  art  predic¬ 
tive  fire  control  systems.  Anti-aircraft  guns,  and  tank  guns  are  the  two  main 
weapon  systems  whose  effectiveness  is  lessened  by  the  presence  of  maneuvering 
targets.  The  application  of  the  methodology  described  in  this  paper  will 
significantly  increase  the  performance  for  these  systems,  especially  when  the 
jinking  or  agile  movement  of  the  target  increases. 

A  gun  fire  control  systems  function  is  to  offset  the  gun  line  from  the 
target  line-of-sight  to  cause  a  projectile  to  intercept  the  target  a  time-of- 
flight  after  firing  the  gun.  Two  classifications  of  target  motion  occur;  non¬ 
maneuvering  and  maneuvering.  The  projectile-target  closest  approach  is  a 
measure  of  miss  distance  and  the  performance  of  the  gun  fire  control  system. 
Non-maneuvering  targets  require  a  constant  offset  between  the  target  line-of- 
sight  and  the  gun  line.  The  magnitude  of  the  required  offset  is  the  product  of 
the  target  velocity  perpendicular  to  the  line-of-sight  and  the  projectile  tine- 
of-flight.  Maneuvering  targets  require  a  time  varying  offset  having  a  magnitude 
related  to  the  non-maneuvering  target  offset  plus  an  additional  amount  related 
to  the  target  acceleration,  acceleration  rate,  etc.,  perpendicular  to  the  line- 
of-sight  combined  with  the  appropriate  functions  of  projectile  time-of-flight. 

The  two  most  familiar  fire  control  systems  use  (1)  target  velocity  and  (2)  target 
velocity  and  acceleration  combined  with  time-of-flight  to  determine  the  gun  line 
offset.  Recent  application  of  sub-optimal  estimation  methodology,  specifically 
Kalman  filtering  has  resulted  in  the  development  of  velocity  plus  acceleration 
controlled  offsets  that  have  significantly  increased  the  performance  of  gun  fire 
control  systems  when  maneuvering  targets  are  engaged.  A  penalty  is  associated 
with  the  estimation  of  accelerations.  Whereas,  it  is  relatively  simple  to 
estimate  velocities,  such  is  not  the  case  for  accelerations.  As  projectile 
flight  times  increase  the  benefits  derived  from  acceleration  estimates  are 
penalized  by  the  errors  introduced  in  the  offset  calculation  by  the  "noise". 
Current  tactical  mobility  indicates  that  maneuvering  targets  require  target 
acceleration  estimates  in  addition  to  target  velocity  estimates  to  calculate 
effective  gun  line  offsets.  As  projectile  time-of-flight  increases  the  benefit 
of  this  offset  philosophy  is  degraded. 

When  the  firing  rate  of  gun  is  relatively  low  with  respect  to  the  cyclic 
motion  of  the  maneuvering  target,  it  is  possible  to  correlate  the  firing  times 
with  the  zero  crossing  of  the  velocity  perpendicular  to  the  line-of-sight.  The 
closest  approach  of  the  projectile  to  the  target  a  time-of-flight  later  cresses 
from  a  positive  to  a  negative  miss  distance  in  this  region  and  both  velocity  and 
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acceleration  are  used  to  determine  the  gun  line  offset.  The  LOS  velocity  zero 
crossing  times  can  be  sensed  with  rate  gyros  on  fire  control  system.  The 
apparent  displacement,  velocity  and  acceleration  of  a  maneuvering  target  is 
viewed  from  a  fire  control  system  situated  at  some  distance.  The  LOS  rate  zero 
crossings  are  detected  by  a  rate  gyro,  causing  the  gun  to  fire  a  projectile. 

The  closest  approach  miss  distance  between  the  target  and  projectile  a  time-of- 
flight  later  indicate  the  zero  crossing  character  of  the  closest  approach  miss 
distance  in  this  region.  Burst  fire  will  be  more  effective  in  this  region. 


DETECTION  AND  INITIATION  OF  FIRING  COMMANDS  FOR  AN  ACCELERATION  PREDICTOR 
FIRE  CONTROL  SYSTEM  ENGAGING  MANEUVERING  TARGETS 

H.  H.  Burke  -  AMSAA 

'  DEVELOPMENT  OF  A  FIRE  CONTROL  SYSTEM  CONCEPT  TO  IMPROVE  PERFORMANCE  AGAINST 

MANEUVERING  TARGETS! 

An  investigation  of  the  performance  of  different  order  predictors 
for  gun  fire  control  systems  engaging  maneuvering  targets  was  conducted  in 
1977  and  reported  on  in  Reference  1.  The  intent  of  that  study  was  to 
demonstrate  the  theoretical  reduction  in  miss  distance  that  was  obtainable, 
as Slimming  perfect  tracking  of  a  maneuvering  target.  The  notion  that  fire 
control  systems  should  consider  second  order  or  acceleration  terms  in  the 
prediction  of  the  gun  line  was  introduced  as  a  possible  alternative  to 
the  present  first  order  or  velocity  term  used  in  current  fire  control 
systems.  In  1978  a  paper  presented  at  AORS  ICVTI  (Ref  2)  extended  this 
work  and  discussed  progress  in  the  technical  area  of  analytical  describing 
maneuvering  targets  for  the  purpose  of  studying  fire  control  system  performance. 
The  three  processes  taking  place  in  all  fire  control  systems  were  identified 
as  tracking,  estimation  and  prediction  i There  is  a  fourth  process,  pointing 
which  was  added  later).  The  three  basic  known  types  of  fire  control  systems; 
manual,  disturbed  reticle,  and  stabilized  sight-director  were  described 
and  compared  in  a  basic  manner.  Prediction  characteristics  of  first  and 
second  order,  systems  were  discussed  and  the  relative  magnitude  of  fire 
control  syttea  errors  from  the  tracking,  estimation,  and  prediction 
processes  were  described.  A  brief  mention  of  a  firing  doctrine  utilizing 
second  order  prediction  to  maximum  advantage  was  presented  'This  idea 

ll;  AMSAA  TR  234,  Aug  78. 

'2)  Given  at  XVI I  AORS,  Nov  78. 


23 


was  not  exploited  until  recently  and  later  in  this  paper  it  will  be  shown 
that  the  fundamental  utilization  of  this  firing  doctrine  is  the  concept 
that  can  be  easily  mechanized  and  applied  to  tank  gun  fire  control  systems, 
with  results  that  offer  major  performance  improvement  over  existing  first 
order  prediction  concepts)  A  paper  describing  the  application  of  sub-optimal 
state  estimation  design  methodology  to  the  development  of  improved  fire 
control  systems  was  the  subject  of  Reference  3.  Some  of  the  earlier  work 
was  contained  in  this  paper  along  with  the  introduction  of  the  stability 
characteristics  of  the  two  main  fire  control  concepts;  disturbed  reticle 
and  stabilized  sight -director.  It  was  argued  that  the  inherent  performance 
of  a  stabilized  sight  was  superior  to  a  disturbed  reticle  because  the  track¬ 
ing,  estimation  and  prediction  processes  were  more  decoupled  and  therefore 
each  process  could  be  designed  to  function  better.  Some  results  presented 
in  Reference  3  indicated  that  the  so  called  "lower  bound"  performance  for 
second  order  predictions  was  superiorxto  first  order  predictors  for  an 

v 

V 
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analytically  generated  path  when  tracking  'ftoi3e  was  introduced  in  the  form 
of  uncorrelated  noise.  This  advantage  decreased  and  for  larger  values  of 
tracking  noise  an  inversion  in  performance  between  the  two  prediction  orders 
occurred.  This  trend  was  used  to  argue  that  the  superior  tracking  per¬ 
formance  of  a  stabilized  sight-director  system  would  enhance  the  benefit  of 
second  order  prediction.  Another  method  of  deciding  on  the  benefits  of  the 
two  orders  of  predictors  was  offered;  time  on  target.  It  was  also  shown  that 
for  the  analytically  modeled  target  motion  being  studied,  that  the  basic 
shape  of  the  prediction  error  for  first  and  second  order  systems  was  dif¬ 
ferent,  and  it  was  believed  that  this  was  the  fundamental  reason  that  the 
time  on  target  indication  showed  that  for  low  tracking  errors,  second  order 

Given  at  Univ  of  ?itt  Modeling  i  Simulation  Symposium,  Apr  79. 
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prediction  was  superior  to  first  order  prediction.  A  separate  AMSAA  Interim 
Note  (Reference  4)  was  written  in  May  79.  The  purpose  of  this  paper  was 
to  more  fully  develop  the  argument  that  predictor  orders  we re  all  iaprotant 
in  detemining  the  distribution  of  errors  in  a  fire  control  system  when 
maneuvering  targets  are  engaged.  The  target  induced  errors  were  claimed 
to  be  much  larger  then  the  tracking  and  estimation  errors  in  themselves , 
but  it  was  admitted  that  the  tracking  and  estimation  errors  had  to  be 
small  to  utilize  the  second  order  prediction,  thereby  reducing  the  target 
motion  induced  errors.  It  was  believed  and  still  is  that  the  disturbed 
reticle  fire  control  configuration  can  not  meet  the  demands  placed  on  it 
that  second  order  prediction  requires.  This  will  have  to  be  studied  in 
detail  before  a  positive  conclusion  can  be  reached. 

Throughout  the  period  of  this  work  we  have  communicated  with  many 
different  groups,  both  within  and  outside  the  government.  One  of  the  items 
we  have  shared  with  them  is  a  data  tape  having  six  time  histories  of  man¬ 
euvering  targets.  Initially  our  work  centered  on  these  paths  but  in  order 
to  remove  the  instrumentation  noise  (the  data  gathering  system  had  no 
relationship  to  a  fire  control  tracking  process,  but  was  rather  a  ranging 
system)  the  paths  were  smoothed  to  surpress  the  high  frequency  noise, 
leaving  frequencies  no  greater  than  1/8  HZ  in  most  cases.  These  data  were 
then  used  to  determine  the  "lower  bound"  performance  measures  mentioned 
earlier.  Based  on  these  smoothed  paths,  the  analytically  generated 
maneuvering  target  generator  was  developed.  For  the  fire  control  studies 
conducted,  to  determine  the  relative  performance  of  first  and  second  order 
predictors,  the  analytical  generated  paths  were  used. 

l4,AMSAA  Interim  Note  C-82,  May  79 
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Some  ocher  groups,  using  the  data  types  supplied  them  chose  to  operate 
directly  on  che  data  following  the  removal  of  some  obvious  data  timing 
irregularities.  This  approach  was  more  in  keeping  with  real  world  data, 
in  that  there  were  certainly  frequencies  greater  than  1/8  HZ  in  the  data. 

The  filtering  algorithms  used  were  therefore  required  to  cope  with  correlated 
high  frequency  noise  superimposed  on  the  basic  path  motion  that  the  analytically 
generated  paths  contained.  The  results  obtained  from  these  studies  were 
not  as  favorable  toward  second  order  prediction  as  the  AMSAA  studies.  In  an 
attempt  to  sort  out  the  reasons  for  this  difference  in  reuslts  and  to 
determine  the  mechanism  that  causes  the  performance  trends  to  shift  when 
empirical  data,  such  as  the  empirical  data  tapes,  are  used  compared  to  the 
analytically  generated  paths,  the  same  data  were  analyzed  by  us,  using  the 
optimal  filtering  designs  we  had  developed  using  the  maneuvering  path 
generator.  Before  discussing  the  findings  of  this  effort  it  should  be 
mentioned  that  in  Reference  1,  for  the  "lower  bound"  studies,  it  was  found 
that  if  the  cyclic  period  of  the  target  motion  was  twice  as  rapid  as  the 
time-of-flight  of  the  projectile,  an  inversion  between  first  and  second 
order  predictors  occurred.  Three  of  the  paths  reported  on  in  Reference  1 
had  this  characteristic  and  the  results  are  described  in  Reference  1. 

It  was  argued  that  realistic  maneuvering  targets  would  not  exhibit 
frequencies  in  this  region _  ( '«  HZ).  The  empirical  data  contained  in  the 
subject  tapes  supplied  to  the  other  researchers  and  now  being  studied  by 
us  probably  exhibited  such  frequencies,  not  from  actual  vehicle  movement, 
but  from  instrumentation  imperfections.  This  was  the  basic  reason  for 
development  of  the  analytical  generator,  but  one  of  its  short  comings  may 
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be  that  the  correlated  high  frequency  tracking  errors  are  not  realisticly 
represented.  A  counter  argument  to  this  idea  is  that  the  data  contained 
on  the  tapes  does  not  reflect  the  correlated  tracking  error  signal  for  a 
given  fire  control  system;  be  it  a  disturbed  reticle  or  stabilized 
sight-niirector  system. 

However  several  of  the  paths  studied  by  some  of  the  other  researchers 
have  recently  been  processed  by  AMSAA  using  the  same  estimation  node!  math 
developed  using  the  analytical  path  generator.  Some  extremely  interesting 
and  exciting  results  have  emerged.  Initially,  time— on-carget  was  used  as 
a  figure  of  merit.  Indeed,  as  some  of  the  other  investigators  have  pointed 
out,  the  second  order  predictor  does  worse  than  the  first  order  in  some 
regions  and  better  in  other  regions  of  the  paths.  Where  is  the  cue  that 
controls  this  performance  reversal?  Inspection  shows  that  as  the  cyclic 
movements  referred  to  in  Reference  1,  are  less  than  twice  the  flight  time 
the  second  order  system  out  performs  the  first  order  system  and  in  periods 
where  the  cyclic  motion  becomes  more  prononunced  at  higher  frequencies , 
{i.e.  higher  apparent  accelerations )  the  first  order  system  out  performs 
the  second  order  system.  This  explains  the  reversal  of  performance 
obtained  using  the  empirical  data.  The  smoothed  data  having  only  the 
basic  path  motion,  with  the  instrumentation  noise  removed,  did  not  exhibit 
a  reversal  (the  frequencies  in  the  data  were  such  that  the  apparent  periods 
were  always  greater  than  twice  the  time-of-flight. ) 

It  was  still  believed  that  in  spite  of  this  turn  of  -vents, seen  from 
analyzing  the  empirical  tapes,  that  some  advantageous  trade-off  existed  for 
second  order  predictors  over  first  order  predictors  using  the  empirical 
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path  data.  This  was  based  on  the  fact  that  the  miss  distance  time  histories 
for  the  two  predictors  were  entirely  different  for  a  given  path.  This  trend 
is  demonstrated  graphically,  in  Figure  S  of  Reference  3  for  the  analytical 
path  generator  studies  and  is  repeated  for  all  of  the  empirical  data  analyzed 
from  the  data  tape  mentioned.  Examples  of  this  are  shown  in  Figure  1  and 
appendices  2  £  3.  As  discussed  in  Reference  3,  the  trend  shown  in  Figure  1 
is  not  as  clear  cut  when  larger  values  of  noise  are  introducted,  but  it 
does  still  exist  in  practice  theory  according  to  the  empirical  paths  we  hava  analyzed 
Thus  far  this  trend  does  persist  with  suprisingly  accurate  timing.  This 
situation  speaks  well  for  the  model  we  are  using  for  the  sub-optimal  estima¬ 
tor  and  for  the  fact  that  it  was  developed  using  an  analytical  path  generator 

model  rather  than  empirical  data.  The  main  effects  described  in  Figure  5 

* 

of  Reference  3  and  Figure  1  of  this  paper  are  (i)  the  miss  distance  for  a 
second  order  predictor  passes  through  zero  a  time-of-flight  increment  after 
the  apparent  velocity  sensed  by  the  fire  control  system  passes  thru  zero 

(2)  the  miss  distance  for  a  first  order  fire  control  system  passes  thru  zero 

a  time-of-flight  after  the  apparent  velocity  passes  thru  its  maximum  value  and 

(3)  when  the  time-of-flight  exceeds  the  time  remaining  before  the  apparent 
acceleration  reverses  direction,  the  miss  distance  errors  for  a  second  order 
predictor  becomes  very  large. 

These  findings  have  been  found  to  exist  is  the  empirical  data  and  can 
be  utilized  to  develop  a  fire  control  system  concept  that  will  greatly 
improve  performance  against  maneuvering  targets,  over  existing  fire  control 
systems  using  no  prediction  or  I'st  order  prediction.  The  concept 
has  the  inherent  built  in  ability  to  revert  to  first  order  automatically 


and  become  second  order  when  required.  It  is  based  on  tbe  firing  doctrine 
concept  mentioned  in  Reference  2  which  wee  ignored  by  us  for  a  long  time 
in  favor  of  the  time  on  the  target  advantage  that  second  predictor  provided 
when  analytical  paths  were  studied.  Noise  detracts  from  the  comparison. 

The  cue  or  sensing  device  required  is  a  rate  gyro  to  sense  LOS  rate.  As 
the  LOS  rate  cresses  zero,  the  fire  control  system  will  fire  and  the  projec¬ 
tile  will  impact  the  target  a  time  of  flight  later.  The  timing  accuracy  of 
these  "firing  windows"  is  very  accurate,  even  when  the  empirical  data  tapes 
are  used.  In  some  cases  there  is  less  time  interval  in  the  "firing  window" 
than  others,  but  for  the  majority  of  opportunities,  based  on  the  logic  of 
the  LOS  rate  passing  thru  zero,  the  time  intervals  vary  between  100  ms  and 
800  ms.  When  the  firing  rate  is  in  the  vicinity  of  1  round/8  sec  for  tank 
guns,  this  rationale  is  reasonable,  in  that  the  best  opportunity  to  hit  the 
maneuvering  target  is  determined  by  utilizing  these  "firing  windows”, 
when  the  target  is  not  maneuvering  too  much,  the  accleration  levels  occurring 
at  these  regions  where  the  LOS  rate  passes  thru  zero  are  not  significant 
and  do  not  corrupt  the  basic  goodness  of  first  order  prediction.  Analysis 
of  the  empirical  data  seems  to  indicate  this  effect. 

It  is  believed  that  these  findings  justify  the  incorporation  of  second 
order  prediction  in  fire  control  systems,  and  that  instead  of  assuming  a 
"fire  hose"  or  arbitrary  firing  time  to  compare  different  predictors,  that 
a  firing  doctrine  similar  to  the  ome  juat  described  be  adopted.  This 
concept  degrades  gracefully  to  the  manual  offset  or  no  predictor  case 
at  the  maneuver  level  lesaens  and  both  velocity  and  acceleration  are  small. 


The  inherently  better  tracking  capability  of  a  stabilized  sight-director 
systea  will  provide  a  more  precise  sensing  of  the  "firing  window"  than  can 
be  obtained  with  a  disturbed  reticle  systea.  Also  this  firing  doctrine  can 
be  used  to  advantage  for  a  rapid  burst  fire  gun  by  firing  the  vurst  in  this 
"firing  window". 

For  cyclic  frequencies  (be  they  real  maneuvering  vehicle  movements  or 
tracking  line-of-sight  motion)  there  is  a  penalty  for  higher  frequencies  and 
extended  time-of-flight  as  shown  in  appendix  i.  This  has  implications  on  gunner 
sight  line  movement  frequency  characteristics  and  long  engagement  ranges  if  there 
is  validity  to  the  serpentine  or  cardoid  model  in  representing  real  world 
vehicle  maneuvering.  The  results  of  the  analysis  reported  in  appendices  2  £  3 
for  ATMT  empricial  paths  433  and  315  indicate  there  is  such  validity  and 
a  second  order  parabolic  predictor  can  predict  serpentine  vehicle  movement  in 
this  region  of  apparent  velocity  nulling  (region  1)  with  errors  described  in 
appendix  1.  When  the  first  order  predictor  outperforms  the  second  order 
predictor,  close  inspection  of  the  region  reveals  that  the  cyclic  frequencies 
are  larger  (-  HZ).  It  is  believed  that  suppression  of  these  tracking 
frequencies  will  be  a  minor  problem  to  cope  with  by  the  fire  control  designer. 
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APPENDIX  1 


2'nd  ORDER  PREDICTOR  AND  SERPENTINE  MANEUVERING  VESICLE 

MOVEMENTS 


The  maneuvering  path  generator  is  built  to  generate  movement  of  a  maneuvering 
vehicle  from  point  A  to  3  as  though  the  vehicle  moves  on  arcs  of  circles 


LO*> 


V 


If  this  movement  is  viewed  from  a  distance,  such  that  it  is  approaching 
the  FCS,  and  if  the  velocity  of  the  vehicle  is  constant  thru  arcs  3t  and  3,, 
the  apparent  positions,  velocities  and  accelerations  perpendicular  to  the  LOS 
can  be  assumed  to  be  simple  harmonic  motion.  Points  1  and  2  are  the  places 
where  the  apparent  position  reverses  direction,  apparent  velocity  crosses 
zero  and  apparent  acceleration  peaks.  A  detailed  discussion  of  this  set 
of  conditions  is  contained  in  a  paper  delivered  at  AORS  XVII,  Nov  79.  This 
"serpentine"  math  model  has  been  in  a  sense  validated  by  the  fact  that  for 
the  AUCT  empirical  data  tapes,  at  least  in  thq  regions  surrounding  points 
l  and  2,  these  apparent  position,  velocity  and  acceleration  trends  exist. 

The  next  question  to  address  is  how  well  does  a  second  order  or  velocity 
plus  acceleration  offset  predictor  work  for  this  region  about  15  2. 

Assuming  the  analytic  model  of  apparent  vehicle  motion  from  the  maximum 
position  point  at  t«0,  the  equation  relating  apparent  movement  to  time 
is  y  *  -a  *  acoswt 
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rCS 


Representing  cos  ujt  in  a  series,  we  have 


where  a  ■  amplitude  or  radius  of  circle 
u  ■  cyclic  frequency  of  movement 
t  =*  time 


COS  uC 


1  - 


(  Mt) 
2  ! 


r  lit)  1  a(  lilt)  *  1 

Then  y  »  -a  +  [a——— -  +  — — -  ] 


2! 


4! 


or 


7 


/  i\  t 1  <a*t : 

•(au*)  — —  **•  .  i  > 

2!  2!  12 


3ut  au*  *  aagnitude  of  apparent  acceleration  at  point  l  (A  ) 


Then 


Aapp  tj  Aapp  tj  w*t * 
2  +  '2  12 


or 


y  ■» 


Aapp 


t*  Jl  -  u*t»  [ 

2  l  12  J 


The  apparent  movement  of  y  for  the  serpentine  math  model  of  the  path 
generator  from  point  1  ( t*0 )  approximates  the  parabolic  second  order 
predictor  concept,  ie: 
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Then  y  •  -App  - 

with  the  tern  representing  the  error  between  the  simple  harmonic  and 

12 

parabolic  curves. 

To  obtain  some  insight  into  the  magnitude  of  this  error  the  following  table 
shows  the  error  term  for  various  cyclic  frequencies  and  tines  of  flight 


t  (tine  of  flight)  1  sec  l.S  sec  2.0  sec 

a  »  1/4  HZ  (1.5  r/s)  0.1875  0.421  0.75 

u  *  1/8  HZ  (.785  r/s)  0.05  0.078  0.21, 

a)  »  1/10  HZ  (.828  r/s)  0.033  0.049  0.13 

Plotted,  this  demonstrates  the  degrading  relationship  that  u( cyclic 
frequency)  has  on  the  ability  of  a  second  order  predictor  to  predict 
simple  harmonic  motion  in  this  localized  region  (point  1). 
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APPENDIX  2 
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Inclosure  la 


ATMT  Path  433 

Tp  *  1.3  sec.;  Filter:  Q  ■  6.,1.,  R  ■  ,01 


Firing  Occurs  At 

2'nd 

Order 

Miss  Distance 

1st  Order 

(M) 

LM) 

2.9s 

1. 3tof 

4.2 

3.09 

4.66 

4.4s 

4.3 

1.62 

3. 98 

4.4 

0.19 

3.25 

4.  3 

-0.48 

2.30 

4.6 

-0.17 

2.90 

4.7 

-0.23 

2.78 

Vx*3.2  M/S 

Vy»4.3  M/S 

3*  tan  VY  ,  39.5° 

vx  ‘ 

(2) 

Firing  Occurs  At 

12.7 

12.3 

+0.196 

-1.94 

11.3s 

12.9 

-  .484 

-2.44 

1 . 5tof 

13.0 

-  .69 

-2.38 

13.0s 

13.1 

-  .31 

-2.59 

13.2 

-1.23 

-2.91 

13.3 

-1.69 

-3.26 

13.4 

-1.94 

-3 . 35 

V7-5 

Vy«2.2 
a  -  i«.3° 
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ATMT  PATH  433 


(  3  ) 

Firing  Occurs  At 

16.2s 

1.3tof 

17.7s 


VY  =  3.3 
3  =32.5° 


Miss  Distance 

2nd  Order  1st  Order 


IM) 

;m) 

17.4 

-  .338 

2.81 

17.3 

-  .491 

2.64 

17.5 

-.  438 

2.57 

17.7 

-  .233 

2.85 

17.8 

+  .148 

3.17 

17.9 

4.  .766 

3.32 

18.0 

+  .346 

3.57 

(  4  ) 

Firing  Occurs  At 


24.3s 
1. 5tof 
23.3s 


26.0 

-4.891 

-7.21 

26.1 

-4.692 

-7.20 

26.2 

-4.216 

-6.95 

26.3 

-2.12 

—3 . 77 

26.4 

-  .072 

—4.  53 

26.3 

-t-2.06 

-3.17 

26.6 

+3.66 

-1.98 

V,^  *  6.6 

V,,  *  S. 4 

3  =  44.1° 


(  3  ) 

Firing  Occurs  At 

30.0s 
1 . 3tof 
31.3s 


4.  7 


*  *  43.3° 


31.1 

2.11 

4.56 

31.2 

2.35 

4.96 

31.3 

2.12 

4.38 

31.4 

1.08 

4.03 

31.5 

.25 

3.51 

31.6 

-  .56 

2.94 

31.7 

-  .98 

2.70 

31.8 

-1.52 

2.34 

31.9 

-2.36 

1 . 7  i 
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4 


(  a  ) 

firing  Occurs  At 
39.3s 
l.Stof 

40.  as 

( acceleration  reversal  during  tof) 


vy  »  .  121 

3  «  1.3° 


Miss  Distance 


2nd  Order 

1st  Order 

i 

IMJ 

IM) 

40.2 

-1.95 

-3.25 

40.3 

-1.12 

-2.84 

40.4 

.45 

-1.96 

40.5 

1.50 

-1.28 

40.6 

1.63 

-  .996 

40.7 

1.81 

-  .732 

40.8 

2.42 

-  .239 

40.9 

3.80 

-  .681 

(  7) 

firing  Occurs  At 

40.7s 

1.4tof 

42.2s 


41. a  *  .434  1.42 

41.9  +  .199  .98 

42.0  *  .197  8.34 

42.1  -  .041  .87 

42.2  -  .850  .74 

42.3  -1.07  .32 

42.4  -1.22  .00 

42.5  -1.18  -  .16 

42.5  -1.01  -  .19 

42.7  -1.04  -  .13 


(  8) 

firing  Occurs  At 
57.8s 


l.Stof 

59.3S 


58.9 

59.0 

-2.85 

-3.13 

59.1 

-2.50 

-3.04 

59.2 

-2.40 

-3.12 

59.3 

-1.81 

-2.86 

59.4 

-  .67 

-2.27 

59.5 

*1.09 

-1 . 22 

59.6 

*2.50 

-2.55 

39.7 

-3.26 

-  .343 
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ATMT  PAIS  433 


(  12) 

Firing  Occurs  At 
ai.2s 
l.stof 
82.7s 


vv  -  1.4 


3  *  10.3° 


(  14) 

Firing  Occurs  At 


VX  *  a'3 

Jr  "  «’.9f 


2nd  Order 


Miss  Distance 

1st  Order 


imJ 

l  M) 

82.2 

-1.89 

-3.41 

82.3 

-1.71 

-3.36 

82.4 

-1.38 

-3.30 

82.3 

-  .44 

-2.78 

82.8 

f  .23 

-2.34 

82.7 

.14 

-2.30 

82.8 

-  .147 

-2.70 

82.9 

-  .35 

-3.26 

83.0 

-1.87 

-3.43 

83.1 

-2.06 

—3. 51 

Firing  Occurs  At 

83.4 

.321 

2.35 

84.4S 

83.3 

—  .360 

1.88 

1.5tof 

83.6 

-2.26 

.771 

33.9s 

83.7 

3.43 

-  .019 

83.8 

3.03 

-  .174 

(Velocity  reversal  during  tof) 

83.9 

-2.33- 

.012 

86.0 

1.195 

.695 

86.1 

.044 

1.18 

VY  «  7.7 

86.2 

1.028 

1.69 

X 

86.3 

1.18 

1.74 

91.6s 

92.8 

-  .74 

-2.69 

l.Stof 

92.9 

-  .88 

-2.75 

93.1s 

93.0 

-1.32 

-3.16 

93.1 

-1.82 

—3 . 0  3 

93.2 

-1.66 

-3.37 

93.3 

-1.12 

-3.14 

93.4 

-  .847 

-3.04 
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APPENDIX  3 


ATMT  Path  315,  Twi*f*r  Vahid* 


Inclosure  lb 


ATMT  313 
Tp  •  1.3  sec 


fl) 


Firing  Occurs  At 

Miss  Distance 

7.8s 

2nd  Order 

1st  Order 

1 . 3tof 

im) 

. 

IM) 

9.3* 

8.9 

9.0 

-2.241 

-  2.98 

9.1 

-1.996 

-  2.88 

9.2 

-1.821 

-  2.36 

9 . 3 

-1.022 

—  2.53 

9.4 

—  .374 

—  2.31 

9.3 

-  .247 

-  2.26 

9 .  S 

-  .382 

-  2.42 

9.7 

-  .  356 

-  2.58 

V,.  3  4*8 

A 

V  *  2  •  9 

V* 

(2) 

Firing  Occurs  at 

15.63 

16.8 

l.Stof 

16.9 

-1.99 

.703 

17.1s 

17.0 

-2.57 

.258 

17.1 

-2.52 

.220 

17.2 

-1.72 

.613 

17.3 

-1.48 

.  633 

17.4 

-1.29 

.714 

17.3 

-  .601 

.969 

17.6 

4-  .449 

1.43 

vx  *  8‘2 

VY  -  1.4 

..  9.7° 

17.7 

*1.692 

2.021 
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ATff  PATH  313 


(  3> 

Firing  Occurs  At 


24.3s 

l.Stof 

25.8 


V 


X 


7.2 


V„  =  3.2 
3  =  24.0° 


2nd  Order 

(T) 

25.3 

25.4 

25.5  -2.259 

25.6  -1.520 

25.7  -1.210 

25.8  -1.042 

25.9  -1.099 

26.0  -  .812 
26.1  -  .866 
26.2  -1.47 


Miss  Distance 

1st  Order 

r*n 


+3.46 

—3.09 

+3.01 

+2.96 

+3.00 

+2.86 

+2.892 

+3.17 


V  4  ) 

Firing  Occurs  At 


29.4s 
1. 5tof 
30.9 


9.6 

.2 


30.4 

30.5 

-  .969 

1. 165 

30.5 

•+  .884 

1.092 

30.7 

.557 

.892 

30.8 

+  .208 

.  721 

30.9 

-  .611 

.234 

31.0 

-1.65 

31.1 

-2.218 

-  .  707 

31.2 

-2.22 

—  7<39 

31.3 

31.4 

31.5 

(5  > 


38.9s 


40.4s 


X 


l  occurs  At 

40.0 

-1.226 

-1.92 

40.1 

-1.077 

-1.86 

of 

40.2 

-1.452 

-2.06 

40.3 

-  .945 

*1 .79 

40.4 

-  .313 

-1.39 

40.3 

+  .143 

-1.12 

40.6 

+  .875 

-  .73 

7*2 

40.7 

+1.560 

.  25 

1.3 

40.3 

+2.75 

2 

—  ■ 

*-• 

o 

• 

to 

o 
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ATHT  PATH  313 


(6  ) 

Firing  Occurs  At 
44.2s 
1. 3tof 
45.7s 


V.c  »  1.3 
V. .  3  .25 


3  =*9.5 


'  7  ) 

Firing  Occurs  At 
59.1s 
1 . 5tof 
60.6 


V  *  4.7 

V,.  a  1.9 

a  =  22. 0C 


(  8  ) 

Firing  Occurs  At 
74.0s 
l.Stof 
75. is 


2nd  Order 

73*7 

45.3  1.28 

45.4  1.17 

45.5  .647 

45.6  -  .235 

45.7  -  .966 

45.8  -1.4 

45.9  -1.53 

46.0  -1.77 


Miss  Distance 

1st  Order 


-2.3 
-2.7 
-2.3 
-2.1 
-1.3 
-  .411 
+2.12 
+  .472 


2.365 
2.369 
2.254 
1.34 
1.46 
1.12 
.308 
.  655 


75.0 

.514 

1.86 

75.1 

.324 

1.74 

75.2 

1.16 

2.21 

75.3 

1.51 

2.47 

75.4 

1.24 

2.44 

75.5 

.974 

2.37 

75.6 

-  .089 

1.73 

75.7 

-  .786 

1.31 

75.8 

-1.25 

1.05 

*  8.0 


ATMT  PATH  315 


(  9  ) 

Firing  Occurs  At 
76.6s 
l.Stof 
78.1s 


V,,  =  4.4 
3  =35.4° 


Miss  Distance 

2nd  Order  1st  Order 


IMJ 

(M) 

77.6 

-6.35 

6.30 

77.7 

-3.23 

6 . 00 

77.8 

-4.15 

5 . 77 

77.9 

-2.29 

4.92 

78.0 

-  .740 

4.11  ■ 

78.1 

-  .025 

3.71 

78.2 

-  .244 

3.36 

78.3 

-  .656 

3.96  1 

78.4 

-  .  665 

3.30 

78.5 

>  .020 

3.51 

( 10> 

Firing  Occurs  At 
81.1s 
1.3tof 
82.6s 


Vx  =6.0 

VY-  6.0 

9  *  45.0° 


82.0 

3.42 

6  •  3  < 

82.1 

o .  06 

6.32 

82.2 

2.41 

6.71 

82.3 

1.48 

6.21 

82.4 

-  .410 

5.33 

82.5 

-1.30 

4.68 

32.6 

-1.58 

3.89 

82.7 

-1.13 

4.04 

82.8 

-  .314 

4.08 

82.9 

83.0 

(U) 

Firing  Occurs  At 
84.3s 


l.5tof 
85. 8s 


8.6 


VY 

3 


3.0 

30.2 


o 


85.4 

—3 . 84 

-7.302 

85.5 

-2.011 

-6,251 

85.6 

-  .673 

-  5.40 

85.7 

>1.711 

-4.01 

85.8 

>3.46 

-  2.37 

85.9 

>4.71 

-  1.90 

86.0 

>5.049 

-  1.23 

86.1 

4.93 

-  .937 

86.2 

* 
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ATMT  Path  315 

(  12) 

Firing  Occurs  At 

37.5s 

1.3tof 

8970 


2nd  Order 

“  "  (M) 


38.6 

3.628 

38.7 

2.963 

38.8 

1.981 

88.9 

.850 

89.0 

-  .678 

89.1 

-1.49 

89.2 

-1.88 

39.3 

—2.13 

89.4 

89 . 5 

—o .  4  ( 

Miss  Distance 

1st  Order 
— TmT - 

5.938 
5.595 
5 . 3ol 
4.722 
3.768 
3.197 
2.817 
2.542 
1. 630 


=  o.o 


3=  31.2 
(.13) 

Firing  Occurs  At 
91.9s 
1 . 3tof 
93.4s 


9o.0  -2.75 

93 • 1  —1 .68 

93.2  -1.29 

93.3  —  ,743 

90. 4  1-1.72 

93.5  3.03 

93.6  4.79 


-5.36 
-5 » 2o i 
— o. Ooo 
-4.603 
-3.088 
-2.227 
-  .316 


V.{  =  8.4 

Vv  =4.3 
27.1° 


(  14) 

Firing  Occurs  At 
96.1s 
l.otof 
S^7 . 6s 


V...  =  3.0 

.  -0 

i  2.  j 


97.1 

97.2 

97.3 

97.4 
97.3 

97.6 

97.7 
97.3 
97.9 
98.0 
98.1 


5.2 

5.6 
5.4 

4.6 
3.8 

2.3 
.934 
.323 

-  .678 
-1.67 


4.98 

5.46 

5.49 

5.28. 

5.02 

4.23 
3.43 

3.23 
2.32 
2.26 
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HYBRID  COMPUTER  SIMULATION  OF 
COMBAT  TANK  DRIVEN  RETICLE  FIRE  CONTROL 

By  John  N.  Groff 

US  Army  Materiel  Systems  Analysis  Activity 


1.  INTRODUCTION 

The  US  Army  Materiel  Systems  Analysis  Activity  (AMSAA)  and 
the  Ballistic  Research  Laboratory  (8RL)  have  jointly  developed  a  man-in- 
the-loop  hybrid  computer  simulation  of  a  combat  tank's  turret/weapon 
stabilization  drives  and  those  portions  of  the  fire  control  that  are 
necessary  for  lead  angle  generation.  The  model  that  has  been  imple¬ 
mented  on  the  BRL  EAI  690  Hybrid  Computer  is  linear  although  it  does 
possess  the  following  non-linear  features: 

•  Hand  Control  Dead  Zones 

•  Response  Rates  versus  Handle  Control  Deflection 

Currently,  the  simulation  is  still  undergoing  development  of 
the  gunner's  oscilloscope  display,  scaling  of  the  input/output  variables, 
and  programming  of  the  time  series  software  necessary  to  reduce  the  data. 
However,  qualitative  rather  than  quantitative  information  concerning 
gunner  tracking  performance  is  being  provided  by  the  simulation. 

2.  BACKGROUND 

In  1978  AMSAA  and  BRL  jointly  undertook  the  responsibility  of 
investigating  whether  proposed  modification  to  a  driven  reticle  continuous 
lead  insertion  system  had  a  significant  effect  in  reducing  weapon 
pointing  errors.  The  basic  problem  was  a  response  mismatch  between  the 
reticle  projection  unit  (RPU)  of  the  gunner's  sight  and  the  turret/weapon 
stabilization  drive  which  resulted  in  the  weapon  lagging  the  RPU.  The 
resulting  offset  appeared  in  the  sight  as  apparent  tracking  error  which 
the  gunner  attempted  to  null  out.  The  problem  tended  to  be  further  ag¬ 
gravated  if  evasive  or  sinusoidal  target  motion  was  introduced.  This 
tended  to  produce  an  oscillatory  instability  or  "rubberband”  tracking 
error  effect,  fhe  investigation  took  the  form  of  a  paper  analysis 
using  availab’e  system  description  documentation.  The  results  of  this 
paper  analysis  are  presented  in  Section  3  and  served  as  the  basis  for 
development  of  the  AMSAA/BRL  Hybrid  Simulation.  Further  motivation  for 
the  hybrid  work  was  provided  by  AMSAA' s  delivery  accuracy  efforts. 

AMSAA  has  the  responsibility  of  determining  quasi -combat 
hitting  probabilities  for  armored  vehicles,  and  the  use  of  engineering 
simulation  techniques  have  proved  to  be  useful  in  performing  this 
work.  Currently,  engineering  simulations  for  the  M60A1,  M60A3,  and 
XM1  tanks  are  being  developed  in  which  all  major  components  of  these 
tanks'  fire  controls,  weapon  stabilization  drive  loops,  suspension 
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system  as  well  as  terrain  profile  and  human  characteristics  will  be 
modeled  and  integrated  into  a  single  simulation. 

Previous  AMSAA  man-in-the-loop  studies  have  indicated  that 
the  gunner  model  used  in  these  simulations  may  not  be  entirely  satis¬ 
factory  for  delivery  accuracy  work.  Part  of  the  problem  stems  from 
the  fact  that  the  McRurer/Krendel  gunner  model  which  will  be  used  does 
not  possess  adaptable  features,  i.e.,  the  ability  of  the  gunner  model 
to  make  internal  adjustments  in  gain,  bandwidth,  and  neuromuscul ar 
reaction  time  in  response  to  target  motion  and  characteristics  of  the 
tracking  system  in  which  the  gunner  is  operating.  Another  inadequacy 
of  the  present  model  is  that  it  is  a  purely  deterministic  model  lacking 
any  residual  randomness. 


3.  RES'uLTS  CF  THE  PAPER  ANALYSIS 


Figure  1  depicts  the  block  diagram  of  the  driven  reticle 
system  that  is  being  used  in  this  study. 
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FIGURE  1.  DRIVEN  RETICLE/CLOSED  LOOP  CCNFIGURATI ON 

A  McRurer/Krendal  man-model  was  interfaced  with  this  system  ana  us 
in  the  analysis. 

The  general  human  transfer  function  used  is  of  the  form: 


-TS 

Kg  (S  +  s>2)  e 

G(s)  =  - - 

(S  +  W]_)  (s  +  013) 

Where,  t  is  a  neurological  time  delay;  K,  ^3,  oj 2  are  values 
that  tho  human  adjusts  for  the  task  at  hand.  Typically,  he  s' 

these  to  achieve  a  loop  crossover  frequency  of  3  rad/sec  wrn  : 
phase  margin  of  35  to  45  degrees.  u»3  is  associated  with  tne 
human  bandwidth. 

For  this  investigation,  the  following  nominal  values 
used;  sjj  »  0.0  rad/sec,  W2  *  0.667  rad/sec,  013  =  6.667  rad/ sec, 
and  r  *  0. 1  sec.  The  human  gain,  K,  was  established  on  the  basis 
subsequent  analysis. 

Specifically,  open  loop/root  locus  techniques  were  a,..'! 
the  overall  system  for  the  purpose  of  determining  relative  scaur 
this  system  with  the  gunner  model  included.  On  the  basis  of  thes. 
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analyses,  the  open  loop  gain  was  adjusted  to  provide  a  gain  margin  of 
3DB  at  the  crossover  frequency.  Figure  2  depicts  a  further  simplifi¬ 
cation  of  the  overall  system  shown  in  Figure  1.  Essentially,  the  TACH 
feed-ahead  loop  shown  in  Figure  1  has  been  eliminated  and  replaced  with 
two  parallel  branches;  an  RPU  loop,  and  a  Turret  drive/stabilization 
loop.  In  addition,  the  common  gain  in  both  these  loops  has  been  com¬ 
bined  with  the  human  gain,  optical  sight  gain,  and  the  handle-bar  gain 
to  yield  a  total  system  gain,  Ks.  The  open  loop  transfer  function  for 
the  system  depicted  in  Figure  2  is  defined  by  Equation  3.2. 


Figure  2.  SIMPLIFIED  AZIMUTH  BLOCK  DIAGRAM 

Where, 

Gl  (s)  »  Ks  (s  +  0.667)  EXP  (-ts)/  s(s  +  6.667) 

G2  (s)  3  (tf/(l  ♦  0.39s)  Q  (.3,  21.9))  ((1/  D  (s)  Q  (.45,  18.0)-1) 
,'1/(1  +  0.23s)  First  Order  Demod  1 

o  (s)  \ 

Cl/Q  (.7,  100)  Second  Order  Demod/ 

S01  (s)  •  Gx  (s)  (G2  (s)  +  G3  (s))  3.2 

The  subscript,  ol ,  denotes  open  loop.  The  open  loop  gain,  Ks,  has  been 
computed  for  the  cases  where  D(S)  was  either  a  first  order  or  second 
order  demod  filter,  i.e., 


P-  (S) 

First  Order  Filter 
Second  Order  Filter 


Open  Loop  Gain 
6.3 


12.5 

With  these  gains,  the  respective  closed  loop  frequency 
responses  were  generated  for  the  systems  having  either  a  first  order  or 
second  order  demod  filter.  Figure  3  depicts  the  bode  frequency  responv' 
for  these  systems.  It  is  apparent  that  substantial  improvement  woul i 
be  realized  from  a  system  using  a  second  order  demod  filter  in  place  :: 
the  first  order  filter. 
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figurti  3  Driven  Retide/CIowd  Loop  fr«CM»ocy  Rssponj* 


In  order  to  quantify  this  improvement,  the  following  star 
techniques  employing  power  spectral  density  of  target  motion  were  use:. 
Equation  3.3  represents  a  one  sided  power  spectral  density  model  of 
target  lateral  linear  acceleration,  y. 


r..  («») 
x 


2o^a  (u£  +■  (a  ^  +  wc^) 
w4  +  2 (a  ^  +  {a  ^  ♦  u^)2 


3.  j 


51 


tAG.  GfcGkfcfcS 


Where, 

a  and  ut  are  positive  constants 
1  -  <3  •  1,2,  3  meters/sec^ 
w  *  Angular  frequency 

The  basic  technique  Is  to  play  the  target  PSD  through  an  error 
transmissabil ity  relationship,  l.e. , 


9g  -  0L  .  S0i  (s)  „ 

-  -  G  ,(s) 

®t  1  +  So]  (s) 

where  Gcj  (s)  *  Closed  Loop  Transfer  Function 
Sol  (s)  *  Open  Loop  Transfer  Function 


e  *  <h  *  (Qg  -  ®l)  “  ®t  -  Gd  (s)«t 
M  (1  -  Gci  (s))0t 
.  Sol(s) 

a  1  ..... .  Qf 

l  *  Sol (s) 


1  +  S01  (s) 

Equation  3.8  may  be  rewritten  in  the  following  form: 


S*~  1  *  Sol (s) 


For  sufficiently  small  0^, 

9t  -  X  . 

TT 

Equation  3.9  now  becomes 


7  i  +  s0l(s) 
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3.11 


In  terms  of  PSD  relationships.  Equation  3.10  becomes, 

r  (“»)  *  r..(w)/  (1  +  G0]  (ito)  G0i ( i u») )  w4  R^ 
e  X 

Integration  of  r  (u)  over  the  frequency  domain  yields  the  variance 
e 


of  the  error,  i.e. , 


a  ^  m  ____  j  r  (u)d(ii. 

ec  2*  e  3.12 

Table  3.1  presents  the  error  one-sigma  values  for  the  two  different 
systems  as  a  function  of  lateral  acceleration.  A  1500  meter  range  was 
used. 


TABLE  3.1 

CLOSED  LOOP  ERROR,  a 

(R  =  1500  METERS 

ee 

Lateral 

First  Order 

Second  Order 

Accel erati op 

Demod  Filter, 

Demod  Filter, 

meters/sec. 

a 

a 

a 

X 

c 

e 

1.0 

0. 280 

0.135 

2.0 

0.559 

0.270 

3.0 

0. 839 

0.406 

4.  DESCRIPTION  OF  HYBRID  SIMULATION 
4. 1  Oscilloscope  Display. 

The  AMSAA/BRL  simulation  presents  a  fixed  base  display  to  the 
gunner  whereas  the  actual  system  has  the  ability  to  slew  the  turret. 

In  terms  of  task  description,  the  actual  systems  might  be  viewed  as  a 
pursuit  tracking  task  while  the  simulation  is  compensatory.  Figure  4 
shows  the  oscilloscope  display  presented  to  the  gunner  which  consists 
of  a  point  target  and  a  cross  hair  reticle.  Both  the  point  target  and 
reticle  possess  the  ability  to  be  offset  from  the  center  of  the  scope 
display.  This  is  accomplished  by  commands  received  from  the  EAI680 
analog  consoles.  These  commands  may  be  generated  in  a  variety  of  dif¬ 
ferent  ways  depending  on  how  the  various  angular  offsets  depicted  in 
Figure  1  are  formed  into  different  display  commands. 
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FIGURE  4.  DRIVEN  RETICLE  DISPLAY 
Consider  the  residual  error  9r  which  is  presented  to  the  gunner, 

9R  *  9t  +  ®L  '  9g  1 

In  terms  of  the  control  handle  output,  ©g^,  Equation  4.1  takes  the  form; 

s)  D(s)  GRpu(s)  hn  -  G,  (s)  : 


5R  =  0t  +  GLt  GRPU  ^  0g  "  Gg  ^ 
zf  c 


4.2 


Where, 


Or  *  Angular  Offset  between  Point  Target  and  Reticle, 
©t  *  Angular  Offset  of  Target,  Motion 
»  Commanded  Gun  Rate 


GLIs)  -  Tf  Kf 

- ,  where  t«  is  time  of  flight, 

T^+l 


gtac($) 

D{S) 
GftPU  (*) 


Ktacs 

9 

*Gear 

Demodulator  Transfer  Function, 
Kr/Q  (;r,  «r),  and 
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Cl  -  ©t, 

c2  -  0L  *  Gi  (s)  Grpu  (s)  QgC, 
f 


C3  -  Gg  (s)  Gjac  (s)  Gtt  (s)  D(S)  Qgc 

C4  3  Cg(s)  OgC. 

Let  Cs  represent  the  spot  command  and  CCf,  the  cross  hair  command.  Var¬ 
ious  ways  exist  of  summing  the  various  commands  (C1.C2.C3.C4)  into  the 
spot  command.  Cs  and  cross  hair  command  Cch-  Three  ways  of  summing  the 
commands  are  presented  in  this  paper. 

Method  1. 

Cs  3  Ci  +  C2  -  C3  -  C4 
Cc^  -  0.0,  and 
Method  2. 

Cs  3  Ci 

Cch  3  C2  ■  C3  -  C4. 

Method  3. 

Cs  *  Ci  -  C4 
Cch  3  C2  -  C3 

The  first  method  results  in  a  purely  compensatory  type  of  tracking  task 
presented  to  the  gunner.  This  proved  unsatisfactory,  since  none  of  the 
"rubberband"  phenomenon  was  exhibited.  The  gunner  simply  nulled  out 
the  residual  error. 

The  second  mechanization  presents  a  pursuit  tracking  task  to 
the  gunner  in  which  he  must  chase  the  point  target  over  the  display 
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screen:  whereas,  in  the  actual  tank  the  pursuit  track  task  results  in 
the  turret  rotating  and  the  reticle  staying  close  to  the  center  of  the 
sight  display. 

Method  two  has  not  yet  been  mechanized  and  checked  out  and 
may  prove  unfeasible  for  certain  types  of  target  motion. 

Method  three  is  currently  mechanized  and  will  be  demon¬ 
strated.  This  mechanization,  while  compensatory,  does  simulate  the 
" rubberband"  phenomenon  and  can  accommodate  all  types  of  target  motion. 

4.2  Overview  of  Hybrid  Computer  Simulation. 

Figure  5  provides  the  reader  with  an  overview  of  the  AMSAA/ 
BRL  Hybrid  Computer  Simulation. 


FIGURE  S  OVERVIEW  OF  HYBRID  SimilATIOH 


All  of  the  major  computing  components  have  been  identified  and  for  the 
most  part  their  functions  are  self-evident.  However,  some  classification 
concerning  the  model  that  has  been  programmed  and  the  planned  use  of 
the  output  data  is  needed.  It  is: 
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a.  The  turret  drive  and  fire  control  model  shown  in  Figure  1 
of  this  paper  has  been  programmed  on  the  EAI  580  analog  computer.  An 
elevation  drive  model  has  also  been  implemented.  For  the  time  being  a 
very  simple  target  path  generator  has  also  been  programmed  on  the  680. 
The  following  equation  defines  the  evasive  part  of  the  target  path 
generator . 

0t(s)  *  G(s)  Wn(s)  4.4 


Where 


G(s)  -  K/(ts  +  1) 

Wn  *  White  Noise 
0-t  =  Target  Angular  Velocity 

The  constants  K  and  r  have  been  chosen  so  as  to  provide  zero  to  peak 
amplitude  of  12  mils/sec  with  a  bandwidth  of  0  to  .2  Hz.  The  simula¬ 
tion  also  can  be  exercised  against  a  constant  velocity  target. 

b.  The  analytical  objectives  of  this  effort  are: 

(1)  Confirm  the  results  obtained  from  the  paper  study 
of  Section  3.0. 

(2)  Create  a  data  base  suitable  to  refine  the  gunner 
model  currently  being  used  in  digital  engineering 
simulation. 

Initially,  it  is  planned  to  develop  a  residual  error  model  which 
describes  the  random  portion  of  the  human.  Box  Jenkins  and  Parameter 
Identification  techniques  appear  suitable  for  this  purpose.  Ultimately, 
a  simplified  version  of  the  Kleinman  (Kalman  Filter/Predictor)  model 
would  be  the  desired  goal  of  this  work. 
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ROBUST  AUTOREGRESSIVE  MODELS 
FOR  PREDICTING  AIRCRAFT  MOTION  FROM  NOISY  DATA 

Stephen  F.  Huling  and  Max  Mintz  Walter  Dziwak  and  Stanley  Goodman 

Dept,  of  Systems  Engineering  US  Army  Armament  Research  and 

University  of  Pennsylvania  Development  Command 

Philadelphia,  Pa.  19104  Dover,  New  Jersey  07801 

INTRODUCTION ; 

Traditionally  fire  control  prediction  algorithms  have  been  based  upon  third- 
order  system  models  in  each  coordinate  (1) .  In  a  previous  paper  by  the  authors 
(2),  the  concept  of  using  higher  order  autoregressive (AR)  models  was  introduced. 
It  was  found  that  in  a  noiseless  environment  these  models  provided  robust 
predictors  which  could  significantly  improve  the  capabilities  of  an  anti-aircraft 
(AA)  artillery  weapon  system  against  a  large  class  of  aircraft  maneuvers  at 
extended  times  of  flight.  This  follow-up  study  examines  the  case  where  the 
observations  made  by  the  weapon  system  are  corrupted  by  noise.  The  achieved 
results  show  that  robust  higher  order  AR  models  still  yield  considerable 
improvement,  especially  in  the  filtering  of  the  sensor  signals.  In  addition, 
an  artificial  flight  path  was  constructed  by  finding  those  autoregressive 
coefficients  which  maximize  the  variance  of  the  estimated  aircraft  velocity. 

It  was  found  that  AR  predictors  modeled  after  real  aircraft  data  work  well 
against  this  "worse  case"  flight  path. 

TEST  DATA  AND  SIMULATION  PROGRAM: 

The  system  models  discussed  here  and  in  Reference  2  were  determined  using  time 
series  analysis  from  test  data.  This  data  was  collected  by  the  US  Navy  from 
flight  paths  flown  by  an  A-7E  high  performance  aircraft  simulating  bombing  a 
defended  ground  target.  The  recorded  data  was  manipulated  and  smoothed  to 
yield  consistent  position,  velocity,  acceleration  and  acceleration-dot (the 
first  derivative  of  acceleration)  data  at  0.1  second  intervals  in  the  XYZ 
coordinate  system.  The  following  descriptions  should  make  it  clear  that 
leaving  the  fourth  and  higher  derivatives  of  position  all  undefined  does  not 
cause  any  of  the  complexities  of  the  maneuvers  to  be  lost.  Although  the 
complete  data  set  consist  s  of  eleven  different  flight  paths,  results  for 
only  three  representative  ones  will  be  presented  here.  These  will  be  the 
same  passes  that  were  used  in  Reference  2  and  the  reader  is  encouraged  to 
refer  to  this  paper  for  a  schematic  figuration  of  these  passes. 

Flight  path  #1  is  representative  of  a  general  class  of  maneuvers  known  as  a 
"dive  toss".  That  is, the  payload  is  released  while  the  aircraft  is  rolling 
and  pitching  which  gives  the  effect  of  the  bomb  being  tcssed.  Hence,  the 
weapon  release  point  can  be  some  distance  from  the  target;  in  this  case, 
the  aircraft  is  1400m  downrange.  The  tracking  data  for  this  pass  is  initiated 


at  a  range  of  5180m — 4020m  downrange  and  3270m  in  elevation.  In  the  middle 
of  the  pass  there  is  a  sharp  3.4g  diving  turn  followed  by  ten  seconds  of 
rolling  back  and  forth  and  then  a  5.5g  maneuver  away  from  the  target.  Less 
than  25  seconds  elapse  between  the  initial  point  and  the  instant  of  minimum 
altitude  (560m),  which  occurs  just  after  the  aircraft  passes  the  target. 

Flight  path  #2  bears  some  similarity  to  pass  #1.  The  significant  differences 
are  that  the  rolling  back  and  forth  in  the  middle  of  the  pass  is  not  as 

pronounced  and  the  weapon  release  point  is  closer  to  the  target  since  this 

pass  is  representative  of  the  common  "dive"  maneuver. 

Flight  path  #3  is  called  a  "pop-up"  maneuver.  Here,  the  pass  is  initiated 
with  the  aircraft  5800m  downrange,  at  an  altitude  of  260m  and  pulling 
almost  6g's  as  it  starts  a  climb.  As  the  aircraft  climbs  over  the  next 
10  seconds  to  an  altitude  of  1440m,  it  rolls  over  so  that  at  the  peak  of 
its  climb  it  is  on  its  back.  It  continues  to  roll  as  it  begins  to  dive  at 

the  target.  It  performs  a  3g  turn  one  way  followed  quickly  by  a  4g  turn 

the  other  before  it  steadies  to  deliver  its  payload.  Then  it  executes  a 
6g  bank  away  from  the  target.  All  this  maneuvering  occurs  within  30  sec. 

To  compare  the  effectiveness  of  the  various  prediction  algorithms,  the  simulation 
program  discussed  in  Reference  2  is  used.  This  program  generates  a  two- 
dimensional  histogram  (prediction  times  vs  miss  distance) ,  which  will  be 
the  means  of  comparison.  The  ballistics  (i.e.  range  as  a  function  of  time) 
of  the  AA  rounds  are  assumed  to  be 

R  =  V  t/(l+0.129t)  (1) 

m 

where  ^(assumed  to  be  1175m/sec)  is  the  muzzle  velocity.  The  results 
presented  in  Reference  2  used  constant  velocity  ballistics) . 

PREDICTION  WITH  NOISELESS  OBSERVATIONS: 


Heretofore,  models  in  fire  control  predictors  have  been  rather  simple.  The 
current  field  peice  utilizes  a  linear  or  position-plus-rate-times-time 
algorithm.  Previous  studies  have  suggested  the  use  of  a  quadratic  algorithm  which 
is  the  above  plus  an  extra  term  comprised  of  acceleration  multiplied  by  one- 
half  of  the  prediction  time  squared.  A  variation  of  this  approach  to 
prediction  is  based  on  a  first-order  Markov  model  of  acceleration  in 
continuous  time ,  i.e., 

a(t)  *  -(*t  a(t)  +  u  (t)  (2) 

Solving  this  general  equation  gives  the  following  discrete  predictors 
m  2 

x  ,.  *  X  ,  +  Tx  +  1 /cj  (exp(-uiT)  +WT-D* 
n+k  n-1  n-1  n-1 
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(3) 


A 

y 

y  ,  + 

•  2 

Ty  +  1  /u>  (exp  (-  u/  T) 

+  WT-I)y 

n+k 

n-1 

n-1 

n-1 

A 

z  =* 

z  + 

Tz  ,  +  l/aj2(exp(-UJT) 

+  U»T-1)  z 

n-t-k 

n-l 

n-1 

n-1 

where  T  = 

(k+l)d 

=»  (k+l)*0.1  second 

This  prediction  algorithm  will  be  used  as  the  benchmark  for  comparison  since 
it  demonstrated  a  prediction  capability  better  than  the  other  predictors  mentioned 
above.  The  results  obtained  using  this  predictor  in  the  simulation  program 
assuming  no  noise  in  the  observations  are  presented  in  Table  1.  Notice  that  for 
prediction  times  of  less  than  one  second  all  the  rounds  for  all  three  passes  score 
"hits"  (defined  as  a  miss  distance  less  than  5m) .  For  prediction  times  of  one 
to  two  seconds,  the  fractions  of  hits  ranges  from  one-third  (passes  #1  and  #2) 
to  one-fifth  (pass  #3) .  The  number  of  rounds  within  fifteen  meters  is  important 
too,  for  close  rounds  can  have  an  effect  on  the  pilot's  resolve  to  carry  out 
his  mission.  However,  Table  I  does  not  show  any  hits,  or  even  a  consistent 
number  of  close  rounds,  for  prediction  times  much  greater  than  two  seconds.  The 
similarity  of  the  effectiveness  of  this  algorithm  across  the  different  flight 
paths  should  be  noted. 


The  authors  in  Reference  2  proposed  that  acceleration-dot  can  be  modeled  as  a 
fifth-order  autoregressive  process  (see  difference  3  and  4),  i.e.. 


a  =  A,  a  ,  +  _  +  )3,a  ,  +  A  a  .  +  £ca  -  +  u 

n  /  1  n-1  r 2  n-2  V  3  n-3  ^4  n-4  P5  n-5  n 


(4) 


where  the  residuals  (u  )  are  uncorrelated  and  zero-mean.  This  AR  model  can  be 
combined  with  the  stanSard  expansions,  (for  X-direction) , 


X  = 

x  i 

+  d  V 

n 

n-1 

n-i 

v  = 

v  , 

+  4a 

n 

n-1 

n-1 

• 

a  = 

a 

+  d  a 

n 

n-1 

n- 

n-1' 


(where  d  =*0.1  second),  to  yield  a  matrix  one-step  predictor 

*A  B 


s 

~n 


-  n-1 


where  s 


n 


S-  k 

(x  ,  v  ,  a  ,  a_,a„  . ,  a 
=■>  n  n  n  n-i  = 

i  a  £/i 

Old 
0  0  1 
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(5) 


(6) 


61 


d3/6  0  0  0  0 

4^2  0  0  0  0 

4  0  0  0  0 


and  £  is  the  zero  matrix.  Although  Akaike's  AIC  statistic  (5)  indicates  that  the 
optimal  model  order  is  larger  than  five,  it  was  found  that  models  with  lower  orders 
provide  a  significantly  poorer  fit  to  the  data  while  models  with  higher  orders  do 
not  give  substantial  improvements.  Therefore,  a  fifth-order  AR  model  is  used. 
Predictions  for  longer  times  of  flight  can  be  made  by  raising  the  partitioned 
matrix  in  Equation  6  to  the  appropriate  power.  For  position  prediction  there  are 
similar  matrices  for  the  Y  and  Z  directions. 


Implementation  of  this  prediction  algorithm  leads  to  the  results  in  Table  II.  The 
robustness  of  this  type  of  predictor  that  was  reported  in  Reference  2  can  be  seen 
in  these  results  since  a  single  set  of  coefficients 

=*  4.029471 

—6.536585 

=  5.313477 

—2.148406 

=  0.340802 

was  used  for  all  three  directions  and  for  all  three  flight  paths  .  This  single 
set  was  derived  via  a  least-squares  criterion  from  the  X-direction  acceleration- 
dot  of  flight  path  #2.  Notice  that  now  almost  all  the  rounds  score  hits  for 
prediction  times  up  to  two  seconds,  twice  as  long  as  for  the  benchmark  predictor. 
Also,  there  is  a  reasonable  amount  of  success  for  prediction  times  of  two  or 
three  seconds  --  about  one-fifth  score  hits  for  passes  #2  and  #3  while  one-third 
doe  for  pass  #1,  plus  more  than  half  the  rounds  are  within  15m  for  all  three 
passes.  For  longer  prediction  times  the  results  differ  from  pass  to  pass. 

Pass  #1  still  has  a  good  number  of  close  rounds  while  pass  #3  has  almost  nothing. 
But  again,  for  prediction  times  up  to  three  seconds  at  least,  the  results  for 
this  algorithm  over  very  different  flight  profiles  demonstrate  a  consistent 
performance  which  is  a  great  improvement  over  the  benchmark  algorithm. 


h 

h 

h 

h 
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PREDICTION  WITH  NOISE-CORRUPTED  OBSERVATIONS 


When  the  observations  are  corrupted  by  noise,  the  problems  become  two-fold;  first 
the  observations  must  be  filtered  to  give  estimates  of  the  various  quantities  at  the 
appropriate  time  instants,  and  second,  these  estimates  must  be  manipulated  to  yield 
reasonable  position  predictions.  Since  the  theory  says  that  a  Kalman  filter  is 
optimal  for  a  linear  system  (6)  and  since  the  assumption  that  acceleration-dot  is 
an  AR  process  allows  a  linear  definition  of  the  aircraft  motion,  a  Kalman  filter 
will  be  used  to  provide  the  state  estimates  in  each  direction. 

Given  a  model  of  the  aircraft's  kinematic  motion  in  one  coordinate 


-  Fs 
— n- 1 


(7) 


where  the  model  residual  vectors  (un)  are  zero-mean  and  uncorrelated  with  covariance 
matrix  Q.  The  position  and  velocity  observations  are  given  by 


where 


(8) 


Then  the  Kalman  filter  can  be  written  as  follows: 


S  =  Fs  +  K  f  Z 

n  —  n-l  -n  L  -n 


where  the  Kalman  filter  gain  matrix  is 

K  -  P  H  '  (HP  h'  +  R  )_1 
-n  -n  —  — n~  -n 


(9) 


(10) 


and  the  error  covariance  matrix  is 


P 


-n 


E  (I  -K  H) P  f'  +  0 
—  ~n-l - n-l~ 


(11) 


The  quantity  in  the  square  brackets  in  Equation  9  forms  the  innovations  process  which 
should  be  a  white  noise  process  if  the  filter  is  properly  tuned.  The  filtering 
algorithm  to  be  part  of  the  fire  control  simulation  program  will  include  three 
such  Kalman  filters  —  one  for  each  direction. 


The  sensors  employed  by  AA  weapons  systems  measure  range  R,  the  elevation  angle  £, 
the  azimuth  angle  A,  range-rate  R  and  the  angle  rates  E  and  A.  Each  of  these 
measurements  will  have  some  independent  noise  in  it.  The  results  to  be  presented 
were  produced  assuming  gaussian  noises  with  the  following  one-sigma  values : 
range-2m,  elevation  and  azimuth  angles  -  0.5  milrad,  range  rate  -  2m/sec  and 
angle  rates  -  1%  of  the  actual  rate.  The  observations  £  which  include  position 
and  velocity  in  the  X,  Y  and  Z  direction  of  a  particular  Kalman  filter  are  simply 


transformations  of  the  noisy  (R,  E,  A,  R,  E,  A) 


Since  the  work  in  this  study  was  performed  in  the  XYZ  coordinate  system  while 
the  noise  occurs  in  REA,  the  2x2  R  matrices  that  enter  into  the  three  XYZ 
Kalman  filters  must  be  derived  from  the  ocvariance  matrix  S  for  (R,  E,  A,  R, 
E,  A) .  Since  the  noise  added  to  the  different  quantities  is  assumed  to  be 
uncorrelated,  this  matrix  will  be  zeros  everywhere  except  on  the  diagonal  which 
will  be  made  up  of  the  appropriate  variances.  Notice  that  since  the  angle  rate 
variances  are  a  function  of  the  data,  ^  (and  hence  the  R  '  s)  will  change 
with  time.  Now,  the  full  6x6  covariance  matrix  in  XYZ  is  given  by 


444 


where  A  is  the  Jacobian  matrix  of  the  transformation  at  time  tn>  In  order  to 
simplif^  the  calculations,  a  decoupling  of  the  XYZ  direction  filters  is  forced 
on  R  .  For  example,  the  R  matrix  for  X  has  the  variance  of  the  X-position 
noise  and  the  variance  of- She  X-velocity  noise  on  the  diagonal  and  the  covariance 
between  the  two  on  the  off-diagonals.  The  covariances  between  directions  are 
lost.  Preliminary  work  indicates  that  filtering  in  REA  which  would  eliminate  the 
need  for  this  decoupling  will  provide  better  estimates.  However,  it  still 
seems  that  prediction  is  best  done  in  XYZ. 

This  study  will  present  results  obtained  using  two  different  filters.  The 
first  is  based  on  the  type  of  model  used  in  the  benchmark  predictor.  In  fact, 
a  first-order  AR  model  of  acceleration  with  the  coefficient  ot =  0.995  will  be 
used.  Tests  show  that  the  performance  of  this  filter  is  rather  insensitive  to  the 
exact  value  of  this  coefficient.  The  vectors  and  matrices  in  the  Kalman  filter 
formulation  defined  above  are  the  following  for  the  X-direction  case: 

5n  =  (xn'  V  an>’' 

1  d  d2/2" 

F  =  0  1  A 

0  0  «* 


H  -  f1  0  °1 

-  0  1  0 


u  *  (0,  0,  u  ) 
~h  n 
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2 

where  E  Cun  J  =  Only  position,  velocity  and  acceleration  are  being  estimated 

so  only  predictors  that  require  these  variables  can  be  used  with  this  filter. 

The  results  in  Table  III  were  produced  by  the  benchmark  predictor,  which  is  based 
on  a  first-order  Markov  model  of  acceleration,  operating  on  the  estimates 
produced  by  three  third-order  Kalman  filters.  The  results  in  these  histograms 
form  the  benchmark  in  this  scenario  of  noisy  measurements. 

The  second  filter  is  based  on  the  fifth-order  AR  model  of  acce) eration-dot.  This 
means  that  the  state  transition  matrix  £  (now  an  8x8  matrix)  is  the  same  as  the  pre¬ 
dictor  matrix  of  Equation  6.  The  results  presented  here  were  generated  by 
three  Kalman  filters  —  a  different  set  of  coefficients  for  each  direction 
(produced  from  pass  #2) .  However,  the  same  filters  were  used  on  all  three 
flightpaths.  Equation  6  also  defines  the  state  vector  £.  The  8x8  residual  co- 
variance  matrix  will  be  all  zeros  except  for  the  element  q(4,4)  which  is  equal  to  the 
variance  of  the  residuals  of  the  X,  Y  or  Z  AR  model  in  pass  #2.  Again  the 
observations  are  of  position  and  velocity  so  the  H  matrix  is  now  2x8.  Since  this 
filter  produces  estimates  of  position  through  acceleration-dot,  the  benchmark 
and  higher  order  predictors  can  both  be  used. 

Figure  1  contains  the  spectral  density  estimates  of  the  velocity  innovations 
processes  for  the  X-direction  of  pass  #2  for  this  and  the  third-order  Kalman  filters. 

As  was  previously  stated,  if  a  filter  is  properly  tuned,  these  processes  should 
be  white.  It  is  clear  from  the  figure  that  the  higher  order  filter  has  a  flatter 
spectrum  which  indicates  that  it  is  more  closely  tuned  to  the  data  and  hence 
that  it  will  provide  better  estimates. 

Table  IV  contains  the  results  obtained  using  this  eigth-order  Kalman  filter  with 
the  benchmark  predictor.  A  large  increase  in  accuracy  for  short  prediction  times 
(less  than  one  second)  is  evident  in  all  three  flight  paths.  Also,  except  for 
pass  #2,  there  is  a  significant  over-all  improvement  in  effectiveness  for 
prediction  times  between  one  and  three  seconds.  It  seems  clear  that  given 
reasonable  noise  levels  little  can  be  expected  for  prediction  times  greater 
than  three  seconds. 

A  single  third-order  AR  model  of  acceleration-dot  with  coefficients 
fl  *  1.451057 
=-0.725526 
=  0.125C00 

provides  the  basis  for  a  position  predictor  that  was  used  against  all  directions 
and  flight  paths.  This  predictor  together  with  the  three-eigth  order  Kalman 
filters  used  above  gives  the  results  in  Table  V.  Note  there  is  a  further  increase 
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in  the  accuracy  for  short  prediction  times-  But  here,  for  the  longer 
prediction  times,  only  pass  #3  shows  any  significant  improvement  over  the 
results  in  Table  IV.  The  results  for  the  other  flight  paths  are  only 
minimally  better. 

PREDICTION  AGAINST  A  "WORSE  CASE  FLIGHT  PATH” 


One  can  associate  with  a  pfc^  order  autoregression  process  a  characteristic 
equation  whose  roots  lie  inside  the  unit  circle  of  the  complex  plane  if  the 
process  is  stable.  Thus,  for  (4) ,  the  corresponding  characteristic  equation 
is 


^-0 


(14) 


When  the  roots  associated  with  (14)  for  each  of  the  flight  passes  are  plotted 
in  the  complex  plane,  it  is  seen  that  they  congregate  in  clusters.  One  can 
identify  five  clusters  per  coordinate.  Furthermore,  the  groupings  of  the 
clusters  is  remarkably  similar  in  each  coordinate.  This  is  suggestive  of 
a  robust  property  of  the  models  which  resulted  in  the  single  set  of  coefficients 
for  each  coordinate  exhibited  above. 

One  can  partition  each  cluster  and  design  a  fifth  order  A-R  process  which  has 
"worse  case"  properties  from  the  point  of  view  of  the  A-A  artillery  weapon. 

The  resulting  class  of  "worse  case  flight  paths"  formed  by  determining  the 
fi  coefficients  from  the  root  locations  can  then  be  used  to  test  the 
effectiveness  of  the  fifth  order  predictors  based  on  real  flight  data. 

One  can  also  use  the  WCFP  predictor  to  see  how  well  it  performs  against  the 
real  aircraft  data. 

The  WCFP  was  designed  in  the  following  way:  One  root  per  cluster  was  chosen 
to  produce  a  set  of  coefficients  such  that  the  steady  state  variance 
of  the  velocity  estimate  is  maximized.  The  condition  of  maximum  variance 
is  achieved  by  roots  located  on  the  boundary  of  each  cluster.  A  magnified 
view  of  two  clusters  and  the  WCFP  root  locations  is  shown  in  Figure  2.  A 
computer  search  over  each  boundary  determined  the  desired  roots.  The  ground 
track  of  a  WCFP  is  exhibited  in  Figure  3. 

Table  VI  shows  the  performance  of  three  predictors  against  the  WCFP.  Observe 
that  the  a-dot  predictor  with  the  WCFP  coefficients  achieves  the  highest 
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average  hit  probability,  with  the  fifth  order  in  a-dot  with  x  coefficients 
from  flight  pass  10  a  close  second.  The  benchmark  predictor  is  a  distant 
third. 

It  was  found  that  the  WCFP  predictor  also  performs  well  against  the  real  attack 
flight  paths. 

CONCLUSIONS: 


Higher  order  AR  Models  of  acceleration-dot  provide  robust  models  of  high 
performance  aircraft  maneuvers  which  are  much  improved  over  the  models 
based  on  a  first-order  model  of  acceleration.  The  filtering  capabilities 
of  a  fire  control  system  could  be  greatly  enhanced  by  the  incorporation  of 
the  higher  order  models.  Further  enhancement  can  be  achieved  by  the  use 
of  the  higher  order  model  in  the  predictor  as  well. 

These  improvements  are  noise  level  dependent,  of  course.  Higher  noise  levels 
mean  a  bigger  payoff  by  changing  to  the  higher  order  filter  with  less  of 
a  return  for  the  higher  order  predictors.  Lower  noise  levels  mean  the 
reverse  with  the  payoff  for  the  higher  order  predictors  approaching  the 
noiseless  case.  However,  in  the  case  of  lower  noise  levels,  the  higher 
order  filter  is  still  expected  to  yield  significantly  better  estimates  than 
the  lower  order  filter.  In  addition,  it  is  needed  to  provide  the  higher 
order  predictor  with  the  estimated  quantities  that  it  requires. 
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FIGURE  t.  SPECTRAL  DENSITY  ESTHETES  OP  X-DIRECTION 
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TABLE  I.  BENCHMARK  PRF.DICTOR 


Prad. 

tin* 

blCVMD 

(NOISELESS) 

FLicjrr  path  i 
Nunbar  of  round* 
v/  nisa  batwaan 

0-5  -10-15-20  T 

■  loan 

tisa 

0  -“I 

45“ U  0  a  -  45 

TK 

1-2 

15  11  16  0  >  42 

7.66 

2-3 

0  5  6  3  -  37 

29.77 

3-4 

0  0  3  2  -  39 

40.32 

4-5 

0  0  0  0  -  30 

94.36 

To cal- 194 

JI7IT 

FLICirr  PATH  2 

Prad. 

Nunbar  of  round* 

tin* 

w/  ml**  batwaan 

Ha  an 

batwaan 

0-5  -10-15-20  T 

Mi*  a 

U  -  i 

1*9  U  (J  3  -  49 

T75T 

1-2 

16  17  5  3  -  41 

6.61 

2-3 

0  4  2  4  -  37 

49.69 

3-4 

0  0  0  0  -  33 

77.48 

4-5 

0  0  0  0  -  34 

141.31 

Tocal-193* 

49.23 

FLIGHT  PATH  3 

Prad. 

Nunbar  of  round* 

tin* 

w/  ala*  batwaan 

Maan 

batwaan 

0-3  -10-13-20  T 

Mia* 

0-1 

<w*  u  cl  (J  -  44 

1.23 

1-2 

9  29  7  0  -  43 

7.32 

2-3 

i  l  2  2  -  38 

47.00 

3-4 

o  0  0  0  -  38 

117.06 

4-3 

0  0  0  0  -  34 

200.84 

Tot* 1-1 99 

67 .53 

TABLE 

Pred. 

tlna 

batwaan 

II.  AR  A-DOT  PREDICTOR 
(NOISELESS) 

FLIGHT  PATH  1 

Nunbar  of  round* 
w/  ml**  batwaan  Van 

0-5  -10-15-20  T  Miss 

0-1 

IT 

mt 

M3  -  46 

0  07 

1  -  2 

42 

0 

0 

0-42 

1.67 

2  -  3 

12  16 

3 

2-37 

11.97 

3-4 

0 

1 

3 

5  -  37 

33.35 

4  -  5 

2 

5 

3 

4-33 

52  21 

Total-195 

16 .  ho 

FLIGHT 

PATH  2 

Pred. 

Number 

of 

rounds 

tin* 

w/ 

nla* 

batwaan 

lean 

batwaan 

0-5 

-10-15-20  T 

Miss 

0  -  I 

59“ 

0 

mt 

M3  -49 

3 .  .  0 

1-2 

37 

4 

0 

0-41 

2  . 7-* 

2  -  3 

7  ' 

15 

3 

1-37 

13.71 

3  -  4 

0 

1 

4 

3-35 

44 .19 

4-5 

0 

0 

0 

2-31 

109  50 

Total-193 

~43“T4 

FLICHT 

PATH  3 

Prad. 

Nunbar 

of 

round* 

tim* 

w/ 

mi** 

batwaan 

•  lean 

batwaan 

0-5 

-10-15-20  T 

Miss 

0-1 

44 

MT 

0  -  “55 

3.07 

l  -  2 

44 

1 

0 

0-45 

1.86 

2-3 

6 

10 

6 

12  -  39 

12  .  71 

3-4 

0 

0 

1 

1  -  36 

61  54 

4-3 

0 

0 

0 

0-33 

102.17 

Total-r97 

~vrsz 

69 


TABLE  III.  THIRD-ORDER  FILTER  UITU 
BEIK3DIAWC  FRED  I  CTO?. 


TABLE  IV.  EIGHTH -OfcDER  FILTER  WITH 
BEIJCIHARK  PREDICTOR 
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FIGURE  2  Roots  of  Fifth-order 

AR  models  of  X-dircction  Acceleration-dot 
for  all  Flight  Paths. 
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FIGURE  3 


Ground  Trace  of  WCFP 


TABLE  iv  WORST  CASE  FLIGHT  PATH  (NOISELESS) 


FIFTH-ORDER  AR  IN  A-DOT  PREDICTOR  (WCFP  COEFFICIENTS) 


Prediction 

Number 

of  Rounds 

with  Closest 

Miss 

Time 

Approach 

Miss 

(in 

meters) 

between 

: 

RMS 

between 

0-5 

5-10 

10-15 

15-20  20-25 

>25  Total 

Error 

0.0  -  1.0 

31 

0 

0 

0 

0 

0  = 

31 

0.01 

1.0  -  2.0 

43 

0 

0 

0 

0 

0  = 

43 

0.48 

2.0  -  3.0 

30 

8 

1 

0 

0 

0  = 

39 

4.88 

3.0  -  4.0 

2 

12 

8 

4 

5 

4  » 

35 

16.40 

4.0  -  5.0 

1 

0 

0 

1 

3 

24  = 

29 

64.27 

(seconds) 

— 

- - 

— 

— 

- - 

- - 

— 

Total  Rounds  =  177 

Average  Hit  Probability  ■  0.61578;  Total  RMS  Error  =  27.116 


FIRST-ORDER  IN  A  PREDICTOR  (BENCHMARK) 

Prediction  Number  of  Rounds  with  Closest  Miss 

Time  Approach  Miss  (in  meters)  between:  RMS 
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(seconds)  -  - - - - -  - -  - 

Total  Rounds  *  178 
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38 
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3 
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A  DESIGN  METHODOLOGY  FOR  ESTIMATORS  AND  PREDICTORS 
IN  FIRE  CONTROL  SYSTEMS 
BY 

JAMES  F.  LEATHRUM 
CONSULTANT  TO  AMSAA 

ABSTRACT 

The  state-of-the-art  in  the  design  of  Kalman  filters  for  fire  control 
systems  leaves  the  designer  with  several  parameters  to  be  used  to  over¬ 
come  the  effects  of  modelling  errors.  The  fixing  of  these  parameters 
usually  requires  extensive  simulation  and  trial -and-error  searching  for 

satisfactory  operating  conditions.  In  the  process,  the  effects  of  the 
various  modelling  errors  are  easily  confounded  and  the  intuitive  under¬ 
standing  of  target  behavior  often  lost. 

The  purpose  of  the  analysis  reported  here  is  to  establish  a  design 
methodology  which  begins  with  the  allowable  variances  in  miss  distance 
and  leads  directly  to  filter  parameters  for  an  optimal  filter.  Structural 
mismatching  between  the  filter  and  the  actual  target  are  left  to  analysis 
by  simulation.  The  advantage  of  this  methodology  is  that  the  performance 
of  the  fire  control  system  in  terms  of  miss  distance  enters  the  design 
process  at  the  outset  rather  than  as  a  "take  what  we  get"  outcome.  Math¬ 
ematically,  the  process  Involves  explicit  solution  of  the  steady  state 
matrix  Riccatl  equation. 


Conventional  Design  Methodolo 


The  conventional  approach  to  the  design  of  estimators  and  predictors 
for  fire  control  systems  is  best  illustrated  by  the  following  development 
of  models  and  parameters.  One  would  start  by  formulating  target  and 
observer  models  of  the  form. 

(a)  Target  Model 

Vi  *  Vic  *  Vk 

(2)  Observer  Model 
Yk  3  HkXk  *  Vk 

These  models  inuiediately  involve  a  linearization  approximation .  The  target 
model  captures  the  well  defined  motion  in  the  state  transition  matrix,  , 
and  leaves  the  less  defined  part  of  the  motion  to  a  noise  term,  B^U^.  The 
observer  is  usually  a  statement  that  not  all  the  state  components  are  visible, 
and  that  the  observations  are  corrupted  by  an  error,  V(<.  (The  index,  k,  is  a 
discrete  time  index). 

If  one  can  further  approximate  and  by  white  gaussian,  zero  mean 
processes,  an  estimator  of  can  be  formulated  as: 


Xk+1  3  }k  Xk 


xk  *  xk  +  Kk  (Yk  -  Hk  V 


Which  is  the  Kalman  Filter  wherein 


Kk 3  pkHYRk +  Hkp  H V 


P  3  Wk  *  8kQk8k 
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pk 5  V  Wk 

«k  I  E<W> 

Ok  =  E‘VV 

In  the  most  sophisticated  fire  control  systems,  the  target  noise  is 
represer ted  in  a  target  oriented  coordinate  system.  Thus,  the  given  Qk 
will  rotate  as  the  target  moves  which  in  turn  leads  to  a  nonsteady  K^. 

The  Kalman  gains  tend  to  change  throughout  the  estimation  process.  In 
addition,  R,  may  be  range  dependent  which  leads  to  further  variability 
in  Kk. 

In  designing  such  a  filter,  the  implementor  is  left  with  choices  of 
the  maqnitude  of  and  R^  (i.e. ,  |J Qt<U  and  A  conventional  design 

process  would  require  assessing  j(  R|J|  from  the  accuracy  of  the  instrumental 
used  by  the  observer.  Since  llQ^II  represents  unmodelled  behavior,  it  is 
usually  adjusted  to  achieve  some  other  objective,  such  as  white  innovation, 
minimum  ensemble  miss  distances.  Whatever  the  objective,  the  last  pnase  is 
unguided  by  the  theory  and  thus  usually  requires  extensive  simulation.  The 
design  process  for  a  fil ter-oredictor  in  tandom  is  illustrated  by 


-Re-design 


HQk»  i|Rki( 
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A1 lowable 
Miss  Distance 
(i.e.  la  miss) 

Predictor  Algorithm 

The  principles  of  the  design  will  be  illustrated  by  restricting  the 
discussion  to  a  single  dimension  and  further  restricing  models  to 
}  :  Upper  Triangular 

Q^:  Scalar  constant,  q 
Rk:  Scalar  constant,  r 

These  restrictions  do  not  limit  variability  of  the  gains  in  the  final 
design,  but  only  allow  one  to  focus  attention  on  the  magnitudes  of  the 
parameters  in  each  direction.  The  design  process  requires  solution  of 
the  steady  state  filter  operations  which  become 

P  =■  iPi  +  BqB'  -$P H'(R  +  HPHT^Po 

The  solution  for  P  in  terms  of  r  and  q  requires  iteration.  However 
a  closed  form  solution  for  P/q  and  r/q  in  terms  of  band  width  of  the  target 
motion  is  possible.  It  requires  the  observation  from  the  analogous 


78 


continuous  models  that 


p33/,  -  <-v/S3)  pt 

Where  wv  is  the  bandwidth  for  velocity. 

is  the  bandwidth  for  acceleration 

d 

it  =  t^i  -  t^;  the  time  increment  between  observations 
The  required  bandwidths  can  be  assessed  from  the  power  spectrum  observed 
in  field  tests  of  generic  targets.  The  ratios  (P/q  and  r/q)  completely 
specify  the  one  dimensional,  steady  state  design,  but  they  do  not  product 
the  parameters  needed  for  |{  Q^Hand  ll  in  a  multidimensional  design.  The 
required  magnitudes  are  obtained  from  the  variance  of  miss  distance 

-~5.  *fPn/q  *\fP22/q‘tf't'fP33/q  mtfZ/2 


Where  in  an  optimal  design,  the  P  is  interpreted  as  the  variance  of  the 
estimator  error.  The  above  equation  represents  the  variance  propagation 
through  a  second  order  predictor  (t^  is  the  time  of  flight.)  Other 
predictor  algorithms  could  be  used  at  this  point.  Since  the  ratio  P/  , 
is  computed,  by  asserting  a  la  miss,  one  can  directly  determine  q.  Further, 
knowledge  of  r/q  leads  directly  to  r. 

The  logical  outcome  of  this  process  is  the  question  of  whether  an 
observer  with  an  accuracy  on  the  order  of  J*r  is  achievable.  The  power 
spectrum  of  the  target  motion,  and  the  limits  of  miss  distance  in 
conjunction  may  force  the  enhancement  of  instrumentation  technology. 

A  Typical  Design 

The  purpose  of  this  section  is  to  illustrate  the  design  process 
by  starting  with 

1  )  A  Second  order  predictor. 
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2)  8andwidths  of  the  target  motion 

u)y  »  0.128hz 
ujk  =  o.ieohz 

3)  1  amiss  distance  =  1  meter. 

4)  Time  ofHight  =  2.0  sec. 

In  the  order  which  they  would  be  computed,  the  variance  ratios  are 


i,  J 

Pi,j/q 

3.3 

.0737 

2,3 

.02716 

1,2 

.00367 

1,3 

.00534 

2,2 

.0146 

1,1 

.00126 

and  then  r/q  =  1.59X10'3 

In  this  computation,  the  model  coefficents  are 


where  the  i  and  8  parts  of  the  model  are  determined  by  the  data  rate  of 
At  »  0.1  sec.  From  the  la  miss  distance,  the  q  is  fixed  at 
q  *  1.88 

thence  a  2.99  X  1q"3 

Given  the  parameters  of  this  example,  the  technological  conclusion  is 
that  an  observer  with  a  la  accuracy  of  r  *  0.0547  meters  is  needed  to 
achieve  a  1  a  miss  distance  of  1  meter. 

The  observer  accuracy  vs.  miss  distance  is  summarized  in  the  following 

table. 
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la  Miss  Distance,  meters  q  fr,  meters 


1.0 

1.88 

0.0547 

1.5 

4.23 

0.0820 

2.0 

7.52 

0.1094 

2.5 

11.75 

0.1367 

3.0 

16.92 

0.1641 

Conclusions 

The  design  process  is  completely  characterized  by  assessing  the 
proportionality  constant  in 

r  =  C  .  amiss 

where  C  is  determined  by  the  bandwidths  of  the  target  motion  and  by  the 
predictor  algorithm.  The  computation  of  C  may  be  tedious,  but  it  is 
straight  forward  and  free  of  iteration. 

Although  the  impact  of  this  methodology  is  clear  from  the  models 


used  here,  its  full  generalization  is  yet  to  be  worked  out.  More  general 
forms  of  i,  8,  and  H  need  to  be  considered  in  the  interest  of  establishing 
the  theoretical  limits  of  the  design  process. 


AN  ADAPTIVE  Ll'.AD  PR  I'D  I CT  ION  ALCOK  T.T1IM  FOR  JrtANF.UVKK  1  NO  TAUGKT  FXGAGFM1  '.T 

Pak  T.  Yip  &  Norman  P .  Co  Leman 
USA  ARRADCOM 
Dover,  NJ  07801 

ABSTRACT .  An  algorithm  concept  which  processes  target  bearing 
and  range  input  data  and  provides  "optimal"  estimates  of  Lare,et  position, 
velocity  and  acceleration  a  time-of -flight  in  the  future  is  discussed. 

Since  the  algorithm  concept  involves  certain  important  statistical  assump¬ 
tions  about  target  acceleration  dynamic  models,  these  assumptions  will  be 
discussed  in  detail  along  with  several  important  methods  used  in  the  model 
identification  process.  Secondly,  the  filter  algorithm  itself  will  be  di  ;- 
cussed.  This  algorithm  involves  the  parallel  processing  of  target  range  and 
bearing  data  by  several  extended  Kalman  Filters  correspond ing  to  distinct 
maneuver  characteristics  of  anticipated  target  vehicles.  At  time  of  fire 
the  filter  with  the  largest  computed  likelihood  function  is  selected  for 
lead  prediction.  Finally,  results  of  simulation  studies  in  which  actual 
target  path  data  is  used  to  generate  filter  input  data  for  hit  probability 
evaluation  is  discussed.  Comparisons  are  made  between  the  adaptive  al¬ 
gorithm  and  non-adaptive  first  order  algorithms. 

I.  INTRODUCTION.  This  paper  describes  a  multiple  model  adaptive 
Kalman  Filter  approach  to  the  problem  of  estimating  and  predicting  the 
position,  velocity  and  acceleration  states  of  tank  targets  of  varying 
maneuverability.  The  estimation  and  prediction  problem  presupposes  that  the 
range  and  angle  DATA  (measurements  corrupted  by  Gaussian  white  noise)  is 
available.  The  target  dynamics  is  described  by  a  system  equation.  Our 
solution  to  this  problem  is  an  adaptive  algorithm  implementable  in  real  tine 
with  a  microprocessor  to  compute  target  position  a  projectile  tine  of  flight 
in  the  future.  This  study  beg, ins  with  the  selection  of  the  Antitank  Miss¬ 
ile  Test  (ATMT)  Phase  II  data]  to  identify  the  filter  acceleration  models. 
It  consists  of  three  dimensional  (x,v,z)  position  data  recorded  at  approx¬ 
imately  10  samples  per  second.  Maximum  likelihood  identification  method  is 
applied  to  this  data  to  identify  a  finite  set  of  Markov  Acceleration  Models 
which  are  representative  of  a  broad  spectrum  of  vehicle  maneuvers  consider¬ 
ed  likely  to  occur  in  actual  engagements.  These  models  provide  the  requir¬ 
ed  state  variable  description  of  the  target  dvnamics  used  in  the  formulation 
of  the  multiple  model  extended  Kalman  Filter  Algorithm  for  lead  prediction. 
The  extended  Kalman  Filter  is  required  in  this  application  as  a  result 
of  nonlinearities  induced  by  target  coordinate  transformations  and  non¬ 
linear  measurement  equation. 

The  adaptive  lead  prediction  concept  is  based  on  the  simultaneous 
(parallel)  processing  of  the  discrete  extended  Kalman  Filters  corresponding 
to  the  distinct  target  models  identified  from  the  ATMT  data.  The  likeli¬ 
hood  function  associated  with  each  filter  is  computed  up  to  the  time  of 
fire  of  the  weapon,  and  the  filter  having  the  greatest  likelihood  is 
automatically  selected  for  lead  prediction. 
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In  the  present  study,  only  the  azimuth  and  range  information  of  t fie  tar¬ 
get  is  processed  in  the  filter  with  the  target  elevation  considered  con¬ 
stant.  The  performance  of  this  design  is  examined  with  a  Monte  Carlo 
simulation  and  the  sensitivity  of  the  lead  estimates  to  measurement  noise, 
level  of  target  maneuver,  range  sampling  rate,  and  time  of  flight  of  pro¬ 
jectile  are  analyzed  to  determine  the  feasibility  of  using  this  algorithm 
for  fire  control  lead  prediction  against  various  maneuvering  targets. 

II.  DATA  ANALYSIS.  The  ATMT  data  consists  of  six  tracks  pro¬ 
duced  by  a  M60A1  tank,  a  Scout  Vehicle  and  a  Twister  Vehicle  undergoing 
evasive  maneuvers.  The  M60A1  tank  is  capable  of  speeds  of  10  to  16  miles 
per  hour  and  with  a  maximum  acceleration  of  approximately  .3g.  The  Scout 
is  an  armored  reconnassiance  vehicle  capable  of  moving  at  a  speed  of 
15  to  25  miles  per  hour  and  a  maximum  acceleration  of  approximately  . 5g . 
Since  our  only  interest  is  in  modeling  the  acceleration,  the  position  data 
is  sampled  at  a  frequency  of  2  cps  and  twice  differentiated  to  obtain  the 
acceleration  estimates  which  are  then  resolved  into  along-track  and  cross¬ 
track  components.  The  power  spectral  density  of  this  data  is  omputed  by 
the  maximum  entropy  method3  which  assumes  the  data  is  generated  by  an  auto¬ 
regressive  process.  The  power  spectral  desity  S(f)  is  given  by 


S(f)  = 


h  _ 


where  is  the  standard  deviation  of  a  Gaussian  noise  process ;  ot; is  the  i-th 
coefficient  of  the  autoregressive  process;  M  is  the  number  of  coefficients, 
and  the  coef f icients  are  estimated  recursively3. 

The  number  of  the  autoregressive  coefficients  is  usually  larger  than  3 
which  Is  not  desirable  for  Kalman  Filtering.  However,  the  power  density 
spectrum  affords  enough  inf oi mat  ion  for  estimating  essential  poles  and 
zeros  of  a  simplier  model  structure.  Later  the  maximum  likelihood  identif¬ 
ication  program  is  used  to  fine  tune  the  pole  and  zero  estimates. 

The  simplified  model  determined  from  the  spectral  analysis  has  the  follow¬ 
ing  form: 

A( s)  - - — -  <?  (s) 

S2+/9,S  +/3z  t 

where  q(s)  is  the  Gaussian  noise  process;  A(c)  is  the  system  acceleration; 
t,  fi,  and  Ai  are  parameters  to  be  identified  for  the  chosen  tracks  and 
each  of  the  along-track  and  cross-track  formulations. 


III.  DISCRETE  EXTENDED  KALMAN  FILTER.  The  system  and  the  measurement 
equations  are  readily  defined  as  follows: 

in  -  ■  dt)  +  S« 

e.  -  b(xj  +  r. 


i vi&  *  ' 


where  Xk  is  the  system  state  vector  at  the  discrete  time  kdt  in  the  cartesian 
coordinate  system,  the  system  function  containing  all  information  about 
the  system  dynamics,  ^the  plant  noise  vector,  2*  the  measurement  vector, 
h(Xj<)  a  vector  containing  the  true  range  and  azimuth  angle  of  the  target 
position  at  time  kdt,  r^  the  measurement  noise  vector,  and  dt  the  time 
between  two  samples. 

The  necessary  statistics  and  conditions  are: 

C0V(S;  *  it) 

cov(t ;  ,tj)  =  Ri  *ij 
cov(Si,rj)  =o.  Vi,j 

e(x.)  =X. 

cov(X.)  =  P. 

where  £ij  is  the  Kronecker  Delta. 


Given  the  above,  the  discrete  Extended  Kalman  Filter  equations  can  be 
written  as  follows:  The  predicted  state  estimate  vector  is  given  by 

X„.i|k  =  ±  ’  Jt) 

and  the  state  error  a  priori  covariance  matrix  by 

Pk„|k  -  $  Pk  $T  +  Qk 

where  ,  i±(Z,  dt) 


$ 


42S 


x  =xk 


t  -  si  +  &<**+  £  (di)A 

/«s  d  ^ 

Xk  =  ~7r~  *k 

^  _  $  b  Xk 

Xk  -  Xk 


The  updated  state  estimate  vector  can  be  written  as  follows: 


^k+l  ^k  +  ljk  +  K  2-k  +  l 

where  2.k  +  (  3=1  2*+i  ~  h  (  X»c-n|k) 

K  -  PjchJk  Ht(HP^|kHt  +  Rk+«) 

M  __  £bte)  ‘ 

n  te 

h  (*w) 


84 


!#*-«** 


x^,  X2  represents  x,y  position  state  estimates  respectively  in  fixed  cartesian 
coordinates.  The  state  error  a  posteriori  covariance  matrix  is  given  by 

fk  +  l|k  —  K  H  Fk  +  i  jk 
\{ *  i(tk-t)Qs  ir(tk--c)  dz 

where  the  continuous  case  plant  noise  covariance  matrix,  Qs,  is  known. 

The  continuous  time  system  dynamic  equations  used  in  deriving  the  dis¬ 
crete  time  equations  are  given  by 

X|  =X5  »  *2  =  *4 

X3  —  (Xg  Aa  +  X4Ac)/V 

X4  =(X4Aa  -X3Ac)/V 

Xs  —  “7^41  X5  -/Qa.2  Xi  >  X6  =  Xs 

X7  =  -/Sci  X7  -fic2  X*  *  Xg  —X7 


Aa 

Ac 

V 


where  and  X^  are  the  corresponding  X  and  Y  components  of  the  velocity 
vector;  Aa  is  the  target  acceleration  along  the  velocity  vector;  Ac  is  the 
target  acceleration  perpendicular  to  the  velocity  vector. 

With  this  filter,  target  range  and  angle  measurements  may  be  processed  to 
generate  target  state  estimate  recursively.  Before  defining  an  adaptive- 
filter  procedure,  the  parameters  of  the  Markov  model  need  to  be  identified. 


Iv •  LIKELIHOOD  FUNCTION  &  MAXIMUM  LIKELIHOOD  IDENTIFICATION  OF 

PARAMETERS y .  liiven  a  parameter  vector  « ,  the  probability  of 
occurence  of  the  measurement  vector  sequence  Z^can  be  represented  bv  a  multi¬ 
variate  gaussian  distribution. 

?(zk  ;  2<)  =  ?(z*jzk-';«)---f  (lajz'  ;*)  Kz.ss) 

^  1 3  iak-i,  _  exp(-i  S k  zH ) 

t(-»|z  .-)  (d  St  S 

S,  =  H  +  Rk 
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where 


P(Ek«*)  =  the  likelihood  function 

n  =  number  of  elements  in  the  measurement  vector  2^. 


In  order  to  identify  the  best  parameter  vector  o<  to  give  a  maximum  p 
we  can  equivalently  minimize  the  negative  log  likelihood  function: 


_l  3^  cr'^r 
2  Zi 


Since  the  term  (27T  )n^  in  the  likelihood  function  does  not  contribute  any 
interesting  information  it  has  been  eliminated  in  forming  MUZ*;#) .  The 
Gauss-Newton  method  is  used  in  the  minimization  procedure. 


2<j  +  l 


2 Si  -  f  D"' 


M  (Zk;sS) 
b  2 


where 


1  for  this  method,  and  D, 


D  *  E 


( 


} 


the  expected  Hessian 


The  test  for  convergence  is  given  by 

(s<j+i  -«j)T  D  (— j+/  -esj)  < 


V.  PARALLEL  FILTERS  6.  ADAPTIVE  ESTIMATION.  Target  state  prediction 
for  maneuvering  ground  targets  have  never  been  a  simple  task  to  undertake. 

The  major  uncertainty  comes  from  the  target  driver's  (stochastic)  decision 

to  maneuver.  However,  it  appears  there  exists  a  maximum  level  of  maneuver 
that  the  ground  vehicles  studied  can  attain.  This  maximum  level  provides  a 
non-trivial  range  of  dynamic  motion  that  can  be  quantized  to  a  finite  number 
of  maneuver  levels.  In  this  study,  five  filters  are  incorporated  into  the 
multiple  model  filter.  Model  Ml  (Filter  1)  is  a  simple  A  states  constant 
velocity  filter.  The  remaining  A  filters  are  identified  with  various  maneuver 
levels . 

The  adaptive  estimation  is  a  straight  forward  decision  making  process. 
Measurement  in' range  and  azimuth  angle  are  processed  through  the  parallel 
filters.  The  filter  having  the  largest  likelihood  function  is  automatically 
chosen  to  provide  the  best  estimate  for  lead  prediction  and  gun  orders.  Two 
concepts  of  adaptive  prediction  are  examined.  In  concept  A  the  likelihood 
functions  account  for  the  entire  measurement  history  up  to  the  time  of  fire. 
Thus  this  adaptive  prediction  concept  is  good  against  targets  with  constant 
maneuver  level.  In  concept  B,  only  the  last  ten  samples  prior  to  the  firing 
time  are  used  to  compute  the  likelihood  functions.  This  adaptive  filter 
concept  tends  to  be  more  sensitive  to  changes  in  target  maneuver  levels. 

VI.  SIMULATION .  A  Monte  Carlo  simulation  of  100  runs  was  set  up 
to  process  a  number  of  10  second  segments  from  the  ATMT  data  representing 
various  maneuver  levels  for  the  Mb0,  twister  and  Scout  Vehicles.  Ihese 
segments  of  data  are  dltferent  from  those  used  for  the  parameter  idem i i i cat  ion 
tasks  discussed  earlier. 
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For  evaluating  the  system  performance,  the  perpendicular  miss  distance  of 
the  predicted  line  of  sight  from  the  real  target  position  is  defined  as  the 
prediction  error  Ep  in  meters.  The  firing  time  points  are  fixed  for  each 
segment  under  process.  The  performance  indicator  ph  at  each  firing  time 
point  is  defined  as  the  ratio  of  the  number  of  times  that  the  prediction 
error  Ep  is  less  than  1.15  meters  to  the  total  number  of  runs.  Actually,  they 
are  hit  probabilities  considering  the  prediction  errors  alone. 

Assuming  engagement  range  of  approximately  2000  meters,  45°  cross  range 
(across  the  range  vector),  1  s~  range  measurement  error  of  2  meters,  1  £*" 
azimuth  tracking  error  of  0.3  mils,  a  projectile  speed  of  1500  meters  per 
second  and  using  the  adaptive  prediction  concept  A,  the  hit  probability  re¬ 
sults  are  illustrated  in  Figure  1  and  summarized  in  the  following  table: 


Target 

Type 

Number  of  Case*, 
7  Firing  Points 
per  Case 

dean 

Const.  Uelocity 
Prediction 

ph 

Adaptive 
Prod ic  t ion 

M60A 1 

13 

.41 

.49 

Scout 

10 

.27 

.38 

Twister 

8 

.20 

.  26 

For  an  engagement  range  of  approximately  1158,  60°  cross  range,  1 
range  measurement  error  of  3  meters,  1  e-  azimuth  tracking  error  of  0.3  mils, 
a  projectile  speed  of  1158  meters  per  second  and  using  the  adaptive  prediction 
concept  B,  the  hit  probability  results  are  summarized  in  the  following  table: 


Target 

Number  of  Cases, 

dean 

ph 

Type 

7  Firing  Points 

Const.  Uelocity 

Adaptive 

oer  Case 

Prediction 

Prediction 

P160A1 

6 

.51 

.56 

Twister 

6 

.31 

.37 

With'  the  latter  conditions,  the  sensitivities  of  the  system  are  observed 
for  a  particular  maneuvering  segment  as  shown  in  Figure  2.  Figure  3 
illustrates  the  system  range  (hence  the  time  of  flight  of  projectile)  sen¬ 
sitivity.  Figure  4  illustrates  the  system  sensitivity  to  angular  measurement 
noise.  Figure  5  illustrates  the  system  sensitivity  to  range  measurement 
noise.  Figure  6  illustrates  the  system  sensitivity  to  range  sampling  rate. 

VII.  DISCUSSION  &  FUTURE  PLAN.  This  study  has  demonstrated  that 
maneuvering  target  acceleration  may  be  adequately  modeled  as  a  discrete 
set  of  stationary  Markov  processes  whose  parameters  can  be  identified  off 
line.  Parallel  discrete  extended  Kalman  filters  have  been  used  to  success¬ 
fully  process  range  and  angle  measurements.  The  adaptive  selection  of  the 
most  appropriate  filter  at  each  time  step,  based  on  its  largest  likelihood 
function,  has  been  accomplished  on  line.  Representatives  maneuver  patterns 
and  levels  used  in  this  study  were  taken  from  the  ATMT  data  base.  The 
results  from  the  Monte  Carlo  simulations  indicate  that  the  performance  of 
the  multiple  model  adaptive  filter  design  is  generally  comparable  to  a 
filter  which  is  tuned  to  the  target  dynamics  of  that  particular  tracking. 
Interval.  In  particular,  the  results  show  that  the  adaptive  prediction 
consistently  performed  better  than  the  constant  velocity  prediction  with  an 
Improvement  in  prediction  ranging  f  r  >m  10  to  percent. 

Since  the  range  data  is  currently  not  a  uniformly  accessable  measurement, 
the  range  sampling  rate  has  been  examined  as  an  area  of  uncertainty 


together  with  range,  angular  measurement  noise,  and  range  measurement  noise. 
The  results  indicate  that  the  system  performance  for. the  azimuth  channel  is 
heavily  dependent  of  angular  measurement  noise  and  projectile  time  of 
flight  in  terms  of  range,  and  is  not  very  sensitive  to  range  measurement 
noise  and  range  sampling  rate.  The  results  also  indicate  that  higher 
probability  of  hit  can  be  obtained  in  the  cross  range  geometry  than  in 
the  down  range  (coming  down  along  the  range  vector)  geometry. 

Implementation  of  this  filter  algorithm  in  real  time  with  a  state  of 
the  art  microprocessor  is  in  the  planning  stage.  We  have  noticed  that 
Merman's  UD  factorization^ for  the  state  error  covariance  propagation  is  a 
desirable  feature  considering  computation  accuracy  and  stability.  Several 
variations  of  the  existing  filter  algorithm  are  also  under  consideration. 
Finally,  a  complete  real  time  simulation  of  the  fire  control  system  with 
the  auto-tracker  or  human  operator  in  the  loop  and  filter  modifications 
to  improve  maneuver  detection  will  be  subjects  of  our  future  work. 
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APPLICATIONS  OF  DELAY  FEEDBACK  IN  CONTROL  SYSTEMS  DESIGN 


N.  P.  Coleman,  E.  Carroll,  D.  Lee  and  K.  Lee 
US  ARRADCOM 
Dover,  NJ  07801 


ABSTRACT :  Necessary  and  sufficient  conditions  for  exact  state  recon¬ 
struction  using  delays  are  discussed  together  with  an  example  in  which  the 
technique  is  implemented  in  real  time  using  an  8080/8085  microprocessor. 

Also,  a  frequency  domain  technique  for  synthesizing  certain  feedback  control 
laws  with  delays  is  developed  and  several  examples  discussed. 

I_. _ INTRODUCTION:  In  designing  a  control  system  using  optimal  control 

theory  or  classical  frequency  domain  techniques,  one  often  encounters  sit¬ 
uations  in  which  certain  required  signals  or  states  of  the  system  are 
unavailable  by  direct  measurement.  In  modern  control  design  this  problem  is 
usually  handled  by  implementing  some  form  of  reduced  order  or  full  order 
observer  which  provides  an  asymtotic  estimate  of  the  unmeasured  state.  In 
this  paper  a  technique  is  developed  for  exact  state  reconstruction  of  unmeasured 
system  states  using  values  of  the  measured  variables,  their  delayed  values 
and  the  control  variables  on  the  maximum  delay  interval .  Several  examples 
are  discussed  which  demonstrate  the  application  of  this  technique  on  a 
laboratory  servo  system  using  an  8080  microprocessor. 

A  second  application  of  delay  feedback  for  frequency  domain  compensation 
is  also  discussed.  A  frequency  domain  technique  is  developed  for  selecting 
appropriate  gain  and  delay  parameters  for  synthesizing  a  feedback  controller 
using  delays  in  the  output  and  several  applications  as  discussed. 

II.  REAL  TIME  STATE  RECONSTRUCTION  USING  DELAYS:  In  this  section  a 
technique  is  presented  for  exact  state  reconstruction  using  delay  feedback 
of  measured  states  of  a  control  system  and  the  values  of  the  control  input 
over  the  delay  interval.  A  real  time  application  of  this  technique  in  a 
servo  control  system  using  an  8080  microprocessor  is  also  discussed.  For 
simplicity,  consider  the  linear  time  invariant  system: 

x(t)  *  Ax(t)  +  Bu(t)  (1) 

where  x  is  an  nxl  state  vector,  u  is  an  rxl  control  vector,  A  is  an  nxn 
constant  matrix,  and  B  is  an  nxr  constant  matrix.  Let  the  observation  vector 
y(t)  be  given  by: 

y(t)  -  Hx(t) 

where  y  is  a  mxl  vector,  and  H  is  an  mxn  constant  matrix.  Let  o ^hj ^  ■  •  <h- ca 

be  time  delays. 

The  problem  is  to  reconstruct  the  state  x(t)  from  the  measurements  y(t), 
y(t-h),  •**,  y(t-hj^)  and  the  measureable  control  vector  u(s),  t-h^ssct. 

The  following  argument  due  to  D.  H.  Chyung,  Reference  (  1  )  provides  the 
basis  for  a  real  time  state  reconstruction  algorithm  discussed  in  the  examples. 
This  argument  makes  use  of  the  well  known  variation  of  parameter  expression 


for  the  time  response  x(t)  of  the  system  (1)  given  by: 


x(t)  -  e 


.  eA(t" 


1 

hi)x(t-hi)  +Ji 

t-hl 

o 

)  +J  e'AsBu(t- 


*  e^^xCt-hi)  +J  e”AsBu(t+s)ds 

-h£  i  *  1 , 2 ,  *  ■  'X. 

Multiplying  both  sides  of  equation  (2)  by  He~A^i  results  in  the  equation: 

o 

He-A^^x(t)  ■  Hx(t-hj)  +  He-A^^  fe-AsBu(t+s)ds 


y(t-hi)  +  He-  1  )  e-AsBu(t+s)ds 


i  =  1,2, ---l 


in  which  the  right  hand  side  is  completely  known.  Letting  C  denote  the  matrix 
given  by: 


A 

C  -  He 

t 


He_Ahl 


ie-AhJL  ■ 


we  can  now  write  equation  (3)  in  the  form 


Cx(t)  -  Z(t) 


where; 


u 

y(t-hi)  +  He_A^l^e_AsBu(t+s)ds 


z(t)  -  y(t-h2)  +  He~Ah2^ e~AsBu(t+s)ds 

j  -h2 

y(t-hjl)  +  He_Ah*J°e~AsBu(t+s)ds 


la  a  known  mXxl  vector  and  C  is  an  mjjxn  constant  matrix  depending  on  the 
parameters  hp  h2»  •••  hjp  If  the  matrix  C  has  rank  n,  then  equation  {k) 
can  be  written  as: 


x(t)  -  f-c]  ~lCrz{ty 


where  CT  denotes  matrix  transpose. 
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Equation  (5)  has  several  Important  implications;  First,  If  the  matrix  |cTc] 
exists,  then  the  state  x(t)  can  be  exactly  reconstructed  from  the  measurement 
y(t),  its  delayed  values  and  the  input  signal  u(t) ,  o«tsh^  ;  secondly,  the 
C  matrix  depends  only  on  the  delays  h},  ***  so  the  right  hand  matrix  cal¬ 
culation  can  be  performed  completely  off  line.  This  leaves  only  the  relatively 
straight  forward  calculation  of  x(t)  and  a  matrix  multiplication  for  on-line 
microprocessor  computation.  This  latter  comment  is  of  particular  importance 
in  real  time  control  applications  in  which  relatively  low  speed  microprocessors 
are  utilized  for  control  law  implementation.  The  following  result  establishes 
the  condition  under  which  the  matrix  C  has  rank  n. 


Result:  There  exists  a  set  of  n  delays  o^hi-<h2  •••*:hn  a,  for  any  a^o 
such  that  the  matrix  C  has  rank  n,  if  and  only  if  rank  (Q)  -  n,  where 


Q 


H 

HA 

0 

HAn~l 


Proof:  Let  a>o  and  assume  rank  (Q)  -  n.  Then  the  row  vectors  of  the  matrices 
He“Ahhf  [o,a)  contains  n  independent  vectors  since,  if  not,  there  exists  b fc  Rn 
such  that  He-Ahb“0  for  all  he  G>>H  •  Repeated  dif ferentation  with  respect  to 
h  gives  He-Ahb“HA-Ahb“HAn"^ e  “®b»0.  This  implies  that  the  non  zero  vector 
e-Ahb  is  in  the  null  space  of  the  matrix  Q  and  hence  rank  (Q)<n. 


Conversely,  assume  rank  C  ■  n,  then  rank  (Q)  *  n  since,  if  not,  there  exists 
b  +  o£Rn  such  that  Hb  ■  HAb  ••••<*  HAn“ib  «  o.  This  implies  He-Ahb  *  o  for 
all  h  and  hence  rank  C<n,  a  contradiction. 


Example:  Evaluation  of  the  state  reconstruction  technique  given  by  equation  (5) 
was  carried  out  on  an  8080  microprocessor  development  system  which  was  in  turn 
interfaced  with  a  laboratory  servo  system  as  shown  in  Figure  1.  In  this  example 
the  system  state  vector  is  given  by  pen  where  xj  is  the  motor  shaft  output 

*14] 

position  and  X2  is  the  motor  shaft  velocity.  The  measured  signal  is  x;  and  X£  is 
reconstructed  using  equation  (5).  Once  the  software  was  developed  and  debugged 
the  program  was  down-loaded  to  a  single  board  8085  microprocessor  shown  in 
Figure  2,  for  faster  execution.*  The  block  diagram  of  the  servo  system  without 
tack  feedback  is  shown  in  Figure  3. 


i* 


y 


Figure  3 


*  The  8085  configuration  shown  in  Figure  3  is  currently  being  used  to  evaluate 
digital  control  concepts  for  the  XM97  turret  system  shown  in  Figure  21. 


& 
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The  state  space  equation  is  given  by: 


—  “ 

*1 

m 

x2 

0  1  xx  0 

+ 

-2406  -4  x2  2406 

•»  -C1’0] 


The  state  transition  matrix  for  this  system  is  readily  computed  to  be 


(6) 


eAt  =  e_2t 


cos49t  +-=— sin49t 
49 


vr* 


in49t 


cos49t  -  sin49t 


-49. Isin49t 

with  the  associated  C  matrix  of  equation  (4)  being  given  by 

IT 


(7) 


C  - 


■—  — 

H 

He'A* 

- 

0 

_e2h 


e2h(cos49h  -  ~  sin49h)  49  sin49H 
49 


(8) 


with  hj  =  0  and  h2  *  h 

For  values  of  h  4  2nTT,  the  matrix  C  is  non  singular  and  we  may  compute  (CTC)-1CT 

-1  49 
C  directly  as 

1  0 

I 

■1 


49cot(49h)  -  2 


-49e-^csc49h 


(9) 


Using  equation  (5)  one  obtains  the  required  state  reconstruction  equation  for 
x2(t)  in  terms  of  the  measurements  xj(t),  x^(t-h)  and  u(s),  k-hssst. 


x2(t)  »  j49cot(49h)  -  xj(t)  -  (49e”^^csc49h)xj (t-h) 

+  2406  |e21lcos49h  -  2e^llsin49h  f-e^s(sin49s)u(t+s)ds 
l_  49  _ [  <  49 

+  2406  |  -e^*1sin49h|  Ue^scos49s  +  2e^ssin49s)u(t+s)ds 

L — ^5 — zl)  ~~vr 


(10) 


The  implementation  of  this  state  reconstruction  algorithm  was  carried  out  on  an 
8080  microprocessor  with  a  delay  value  h  “  .01  sec.  The  position  output  state 
was  sampled  at  2.2  milisecond  intervals  and  the  accuracy  of  the  A/D  and  D/A 
converters  was  12  binary  bits. 
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t Auction  with  8080/8085 
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Figure  4a  compares  the  actual  tach  output  signal  representing  the  X2(t)  state 
with  the  microprocessor  output  signal  which  attempts  to  reconstruct  xi(t)  via 
equation  (10)  using  only  the  first  two  terms  of  this  expression.  Note:  In 
this  case  equal  weighings  must  be  used  for  x^(t)  and  x^(t-h).  The  effects 
of  measurement  noise  are  readily  apparent  in  this  figure.  Figure  4b  again 
compares  measured  tach  output  with  the  microprocessor  output  signal,  however, 
in  this  case  the  full  state  reconstruction  equation  (10)  is  implemented.  This 
implementation  is  seen  to  give  a  very  accurate  state  reconstruction  which  is 
less  sensitive  to  measurement  noise. 

III.  FREQUENCY  DOMAIN  CONTROL  SYNTHESIS  USING  DELAY  FEED? \CK :  Several 
papers,  (see  Reference  (5  )  and  (6  j)  have  appeared  in  the  recent  literature 
which  address  the  problem  of  feedback  control  using  delays.  Reference  (6  ) 
develops  several  feedback  control  laws  using  delays  in  the  state  and  derivative 
of  the  state  which  are  shown  to  drive  the  full  state  of  the  system  to  zero 
and  keep  it  there.  The  constructions,  however,  have  limited  utility  in  servo 
control  applications  since  they  assume  first  that  the  control  space  has  the 
same  dimenstion  as  the  state  and  all  states  of  the  system  are  accessible  for 
on-line  measurement. 

In  this  section  we  consider  a  restricted  class  of  delay  feedback  controllers 
shown  in  Figure  5.  This  configuration  has  proved  quite  useful  in  turret  and 
servo  control  applications  in  which  G(s)  represents  the  open  loop  transfer 
function  between  the  command  input  and  the  position  output.  The  two  design 
parameters  introduced  by  delayed  feedback  are  seen  to  be  K,  the  feedback  gain, 
and  T,  the  feedback  time  delay.  The  reason  for  chosing  the  two  feedback  gains 
in  the  form  K  and  K-l  differing  by  unity  in  the  general  case,  will  be  made  clear 
below.  The  equivalent  feedback  transfer  function,  H(s),  for  the  system  in 
Figure  is: 

H(s)  *  K  -  (K-l)e-Ts  (II) 

We  may  represent  the  e"Ts  term  by  its  equivalent  Taylor  series  form  as: 

e-Ts  =  1  -  Ts  +  T2s2  -  T^s^  +  ... 

2  6 

The  frequency  band  of  primary  interest  from  a  stability  and  transient  response 
point  of  view  in  |s:  |G(s)H(s)|  f]  or  Q:  ssw^J  where  wc  denotes  the 

gain  crossover  frequency  of  the  compensated  open  loop  system.  Setting  s  =  jw 
and  assumming  IwTl^^l,  we  may  approximate  e~Ts  by  the  first  two  terms  of  its 
Taylor  series  expansion  or; 

e“Ts  *  1  -  Ts  -  1  -  jwT  (12) 

Substituting  (12)  into  equation  (11)  yeilds; 

H(s)  -  K  -  (K-l) (1-jwT)  =  1  +  j (K-l)Tw  (13) 

Since  K->  1  will  be  required,  this  corresponds  to  a  phase  lead  network  (on  a 
first  order  approximation  basis)  in  the  controller  feedback  path.  If  this  phase 
lead  term  is  properly  positioned  in  frequency,  it  will  produce  a  stabilizing 
effect  upon  the  control  systems  unit  step  or  impulse  response  characteristics. 

As  will  be  seen  in  the  examples,  the  time  delay  or  feedback  gain  can  be  adjusted 
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to  provide  any  desired  damping  in  the  system  response.  The  procedure  for  intro¬ 
ducing  a  lead  network,  effect  around  w  =  wc  using  delayed  feedback  can  now  be 
developed. 

First,  choose  wc  such  that  |c(jwc)|  ■  1 

<T 

Second,  select  the  feedback  time  delay,  T,  such  that  Twc<-Cl.  The  choice  of 
Twc  =  .  1  is  reasonable  and  is  used  in  the  examples.  For  this  choice,  the  first 
term  disregarded  in  the  series  expansion  has  magnitude  .005  at  wc  and  rapidly 
becomes  smaller  for  higher  frequencies.  Third,  select  the  feedback  gain  para¬ 
meter,  K,  such  that  the  lead  time  constant  becomes  effective  at  or  near  w  =  wc 
i.e.  (K-l)T  =  1  .  Note  under  this  condition  using  step  1  and  equation  (8),  that; 

wC 

|g(wc)H(wc)|  =  |g(wc)|  |  H(wc) 

< — yr 

~  *  and  K  *  11 

The  delayed  feedback  design  procedure  thus  is  seen  to  be  straight  forward  in 
concept.  The  effect  of  the  particular  delayed  feedback  configuration  discussed 
here  is  to  replace  the  more  standard  tach  feedback  stabilization  loop.  When 
the  delay  time  and  gain  parameters  are  properly  chosen,  system  response  charact¬ 
eristics  may  be  improved  substantially. 

Example: 

We  consider  first  a  simple  laboratory  servo  system  whose  open  loop  transfer 
function,  G(s),  is  given  by; 

G(s)  -  600 

s(l+s)  ( 1 A ) 

4 

The  -3db  crossover  frequency,  wc,  of  the  open  loop  transfer  function  G(s)  is 
56  rad/sec  and  the  delay  time,  T,  is  computed  from  step  3  and  satisfies  10T  =  1 
or  T  *  .0017  sec.  The  gain  K  is  fixed  and  satisfies  the  relation;  56 

K  -  1  -  1  -  1  -  10 

Twc  .1 

Due  to  limitations  of  the  8080  microprocessor,  the  above  design  using  a  delay  of 
1.7  ms  could  not  be  implemented.  The  smallest  delay  which  could  be  implemented 
with  the  8080  was  2.2  ms.  The  performance  of  this  design  for  a  step  input  command 
is  shown  in  Figure  6b.  Figure  6  illustrates  that  the  effective  damping  intro¬ 
duced  by  the  feedback  delay  can  be  further  increased  by  increasing  the  delay 
parameter  T.  The  desired  damping  can  also  be  achieved  by  adjusting  the  gain 
parameter  K. 

To  evaluate  the  effects  of  delay  parameters  which  were  too  small  for  imple¬ 
mentation  on  an  8080  microprocessor,  simulations  were  run  for  values  of  T  r  .8 
msec,  1,7  msec,  2.2  msec,  4.4  msec,  8.8  msec  and  17.6  msec,  using  the  servo 
transfer  function  (14).  These  results  are  shown  in  Figures  7  -  12.  Deficiencies 
in  the  linear  model  of  the  servo  system  are  readily  apparent  since  the  simulations 
indicate  more  damping  than  is  evident  in  the  test  results  of  Figure  6  and  Figure  12 
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indicates  an  instability  with  the  17.6  msec  delay  in  constrant  to  the  over 
damped  response  in  the  hardware  test  shown  in  Figure  6e. 

Example: 

In  this  example  we  illustrate  the  application  of  the  delay  feedback  control 
synthesis  technique  to  the  design  of  a  controller  for  an  XM97  helicopter  turret 
control  system  shown  in  Figure  13.  The  transfer  function  block  diagram  of  this 
system  is  shown  in  Figure  14.  The  -3db  crossover  frequency  for  the  open  loop 
system  (tach  loop  open)  was  computed  to  be  20  rad/sec  resulting  in  a  feedback 
time  delay  of  .005  sec.  The  step  response  of  the  original  XM97  design  is 
shown  in  Figure  15  and  that  of  the  delay  feedback  design  in  Figur.:  16.  The 
latter  design  exhibits  a  dramatic  improvement  with  respect  to  overshoot  and 
settling  time.  This  improvement  can  be  explained  partially  by  the  fact  that  the 
original  system  uses  motor  tachometer  feedback  for  stabilization  while  the  delay 
feedback  design  effectively  uses  actual  turret  position  and  rate  for  feedback 
stabilization.  Figures  17  -  22  also  show  the  effects  of  increasing  and  decreasing 
the  delay  feedback  parameter.  Saturation,  columb  friction  and  deadband  non- 
linearities  are  included  in  the  simulation. 

IV. _ CONCLUSION:  Applications  of  delay  feedback  for  state  construction 

and  feedback  control  are  presented  together  with  simulation  results  and  examples 
of  actual  implementations  using  Intel  8080  and  8085  microprocessors.  These 
examples  demonstrate  the  practicality  of  the  ideas  and  suggest  that  these  tech¬ 
niques  may  provide  a  useful  adjunct  to  the  more  standard  frequency  domain  and 
state  variable  techniques  for  estimation  and  control  applications. 
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